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ABSTRACT 


Dynamic  response  of  metals  and  metallic  structures 
is  reviewed  with  respect  to  current  procedures  and 
understanding  in  order  to  identify  research  needs. 
Emphasis  is  placed  on  response  in  military  environments 
which  involve  very  high  loading  rates,  such  as  those 
due  to  air  blasts  and  underwater  explosions.  Included 
in  this  review  is  discussion  of  basic  inelastic 
deformation  mechanisms,  mechanical  properties,  fracture 
analysis  and  testing,  structural  analysis  (including 
consideration  of  buckling  response)  and  scale  modeling. 
The  principal  conclusions  and  recommendations  are 
given  in  Section  II. 
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I. 


INTRODUCTION 


The  Department  of  Defense  (DoD)  has  a major  interest 
in  structural  evaluation,  design,  and  testing  of  equipment. 
Current  approaches  involve,  to  varying  degrees,  modern 
methods  of  structural  analysis,  linear  elastic  fracture 
mechanics,  and  nondestructive  inspection,  as  well  as  the 
more  traditional  methods  of  impact  testing  of  specimens  and 
scale  modeling. 

This  report  concentrates  on  military  environments  that 
involve  very  high  loading  rates,  such  as  those  due  to  air 
blasts  and  underwater  explosions.  In  these  problems,  valid 
approaches  must  incorporate  the  strain  rate  and  temperature 
sensitivity  of  inelastic  behavior  possessed  by  many  struc- 
tural alloys.  There  is  presently  a considerable  amount  of 
research  underway  to  improve  the  various  technologies  and 
their  interaction  in  dynamic  problems.  However,  the  state 
of  the  art  is  quite  limited,  especially  for  the  problems 
involving  dynamic  elastic-plastic  response  and  failure. 
Furthermore,  in  many  cases  where  theoretical  methods  exist 
there  are  insufficient  mechanical  and  fracture  property 
data  in  the  high  strain  rate  regime  to  permit  implementation 
and  verification  of  the  analyses. 

Being  cognizant  of  these  factors,  the  DoD  requested 
that  the  National  Materials  Advisory  Board  (NMAB)  appoint 
a committee  to  recommend  areas  where  current  procedures 
and  understanding  should  be  strengthened.  Special  attention 
was  directed  to  the  examination  of: 

• Materials  tests  which  reflect  more  advanced 
concepts  in  the  mechanics  of  failure, 

• Developments  which  could  successfully  merge 
material  property  testing  with  modern  methods 
of  computational  structural  mechanics, 

• Developments  which  will  lead  to  quantitative 
analysis  procedures  for  assessing  the  flaw 
and  load  severity  in  defense  engineering 
systems,  and 
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• Development  of  small-scale  testing  and  evaluation 
procedures  which  will  reduce  the  need  for  large 
model  or  full-scale  testing. 

Important  benefits  accruing  from  implementation  of  the 
Committee's  recommendations  are: 

• Reduction  in  testing  costs, 

• More  rapid  assessment  of  materials  feasibility, 
and 

• Development  of  better  materials,  processing, 
and  fabrication  to  cope  with  dynamic  loads. 

As  a practical  matter,  considering  the  very  broad 
interdisciplinary  nature  of  the  study,  this  Committee 
decided  to  restrict  its  assessment  to  metallic  structures 
and  to  deal  with  only  nominal  strain  rates  of  loading  up  to 
10®  seconds"^;  this  latter  restriction  excludes  phenomena 
such  as  those  produced  by  high-speed  projectile  penetration 
and  particle  impact.  Fatigue  also  is  excluded  from  this 
study . 

It  is  believed  that  significant  improvements  in  the 
state  of  the  art  can  be  achieved  through  a proper  under- 
standing and  application  of  the  various  interrelated  dis- 
ciplines which  are  involved:  metallurgy,  mechanical 
property  characterization,  fracture  analysis  and  testing, 
structural  analysis  (including  consideration  of  buckling 
response)  and  scale  modeling.  This  report  addresses  each 
of  these  areas  to  the  extent  believed  necessary  to  respond 
to  the  aforementioned  DoD  request  and  provide  the  basis  for 
the  Conclusions  and  Recommendations  given  in  Section  II. 

(For  the  reader's  convenience,  individual  sets  of  conclusions 
are  also  listed  at  the  end  of  the  most  lengthy  subsections; 
however,  only  the  main  ones  are  collected  in  Section  II.) 
Figure  I-l  depicts  the  way  in  which  the  several  disciplines 
interact  in  the  assessment  of  damage  tolerance  of  structures. 

The  time  variation  of  the  loading  is  obviously  an 
important  consideration  in  dynamic  problems.  Therefore, 
background  information  on  the  nature  of  some  transient 
load  environments  is  given  in  Section  III.  Current 
understanding  of  dynamic  response  of  materials  and 
approaches  employed  in  its  investigation  are  discussed  in 
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Section  IV.  Of  the  various  areas,  dynamic  fracture  mech- 
anics and  related  testing  are  covered  in  the  most  detail; 
the  committee  believes  this  emphasis  is  necessary  because 
of  the  complexity  of  the  subject  and  the  fact  that  there 
have  been  many  important  developments  in  just  the  past  few 
years  that  are  not  summarized  elsewhere. 

There  are  many  important  problems  in  which  the  failure 
mode  is  excessive  deformation  and/or  buckling.  Prediction 
of  this  type  of  failure  for  ductile  materials  subjected  to 
large  impulse  loading  can  often  be  made  quite  simply  and 
accurately,  using  the  assumption  of  rigid-plastic  behavior. 
Inasmuch  as  Jones  (1975,  1976,  1978)  and  Jones  and  Okawa 
(1976)  have  recently  reviewed  both  theoretical  and  experi- 
mental studies  in  this  area,  it  is  not  covered  in  this 
report.  They  discuss  relatively  simple  formulas  which  pro- 
vide physical  insight  into  dynamic  behavior  and  quite 
accurate  estimates  of  a variety  of  responses,  including 
failure  of  beams,  plates,  shells,  and  rings,  and  slamming 
and  ice  damage  of  ships.  The  importance  of  such  dynamic 
plastic  analyses  cannot  be  overemphasized;  they  often  are 
entirely  adequate  in  the  early  stages  of  a preliminary 
design,  provide  guidelines  for  the  selection  a meaningful 
set  of  parameters  for  an  experimental  program,  and  can 
serve  as  special  cases  for  checking  more  general  numerical 
procedures . 

The  reader  will  see  that  different  sets  of  units  are 
used  in  this  report.  Many  figures  are  reproduced  from 
publications,  and  conversion  of  data  and  units  would  be 
rather  time-consuming  and  costly.  Therefore,  we  implore 
the  indulgence  of  the  reader  and  offer  a brief  conversion 
table  between  English,  metric,  and  SI  units  in  Table  I-l. 
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TABLE  I-l  Standard  International  Values  of  Some 
U.S.  Units 


1 in.  = .0254  m 

1 ft.  = .3048  m 

1 Ibf  = 4.448  N 

1 kgf  = 9.807  N 

1 dyne  = 10  E-6  N 

1000  psi  = 6.895  MN/m^ 

1 bar  = E 6 dyne/cm^  = .1  MN/m^ 
1 atm  = .1013  MN/m^ 

1 torr  = 1 mm  Hg  = 133.3  N/m^ 

1 kgf/mm^  = 9.807  ]y[N/m2 

1000  psi  Vin.  = 1.099  MN/m^^^ 

1 Ibf/in.  = 175.1  N/m 

1 Ibm  = .4536  kg 
1 ft  Ibf  = 1.356  Nm 
1 erg  = 1 dyne  cm  = . lE-6  Nm 
1 hp  = 745 .7  Nm  /s 
1 kw  - 1000  Nm/s 

1 Ibm/in.^  = 27.68  E 3 kg/m^ 

1 gm/cm^  = 1000  kg/m^ 

1 centipoise  = .001  N s/m^ 

1 Ibf  s/ft2  = 47.88  N s/m^ 

= 344  m/s 
*^water  ~ 1470  m/ s 
^steel  “ 5000  m/s 
^light  ~ 299.8  E 6 m/s 

g = 9.807  m/s^ 
y = 66.70  E-12  Nm  ^/kg^ 
or  m^/kg-s^ 
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II.  CONCLUSIONS  AND  RECOMMENDATIONS 


These  conclusions  and  recommendations  represent  the 
composite  opinion  of  the  committee.  Every  effort  was  made 
to  be  objective  in  the  evaluation  of  the  recent  literature 
and  current  activities  in  the  field. 

As  an  aid  to  the  reader,  section  numbers  in  which  the 
relevant  background  discussion  appears  are  indicated  in 
parentheses . 

A.  Conclusions 


1.  A reasonable  understanding  of  the  processes  govern- 
ing constant  strain  rate  and  temperature  effects  on  yield 
and  stable  plastic  response  has  been  achieved.  Experimental 
techniques  and  data  are  available  for  many  structural 
materials  in  the  ranges  of  strain  rate  and  temperature  of 
engineering  interest.  However,  existing  constitutive 
equations  and  data  are  not  adequate  for  predicting  inelastic 
response  under  conditions  of  unloading  and  reloading  into 
the  plastic  state.  (Sections  IV. B. 1-3.) 

2.  One  of  the  oldest  ideas  in  fracture  control  is  the 
certification  of  materials  for  service  by  testing  under  high 
rates  of  loading.  Impact  tests  such  as  the  Charpy  test,  the 
drop— weight  tear  test,  and  the  dynamic  tear  test  have  been 
(and  continue  to  be)  widely  used.  These  tests  were  original- 
ly intended  to  give  qualitative  information  only,  such  as 
the  transition  temperature.  Recently,  however,  it  has 
become  widely  realized  that  it  may  be  possible  to  extract 
useful  quantitative  information  from  the  results  of  impact 
tests.  What  is  not  yet  generally  recognized  is  that  the 
conventional  analytical  interpretation  of  such  tests  via  a 
completely  static  point  of  view  can  be  substantially  in 
error.  Recent  research  work  using  dynamic  (i.e.,  inertia 
forces  included)  analyses  described  in  this  report  have 
clearly  shown  this  for  precracked  charpy  specimens  and  for 
the  drop-weight  tear  test.  Hence,  there  is  a strong  need 

for  additional  dynamic  analyses  to  accompany  high  strain 
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rate  testing.  These  will  likely  lead  to  improved  specimen 
design  and  test  techniques  to  evaluate  fracture  properties 
of  more  direct  use  in  engineering  structural  design 
(Section  IV. B. 4.) 

3.  The  two  main  approaches  in  dynamic  fracture  analysis 
and  testing  are:  1)  the  dynamic  propagation  of  a crack 
initiated  under  quasi-static  loading  conditions,  and  2) 
initiation  of  crack  growth  under  dynamic  loading.  Obviously, 
a full  and  proper  treatment  of  dynamic  fracture  must  combine 
dynamic  initiation,  propagation  and  arrest  of  cracks;  such 

a capability  does  not  yet  exist.  (Sections  IV.  B.4-6.) 

4.  Provided  the  loading  is  known,  or  can  be  established 
if  structural  coupling  exists,  prediction  of  dynamic 
response  without  flaw  growth  and  large  deformations  can  be 
made  for  many  engineering  structures.  (Section  IV. B. 6.) 

5.  Large  and  full-scale  models  are  currently  used  for 
structural  evaluation  despite  the  high  costs  and  time  con- 
sumption, generally  because  there  is  a lack  of  confidence 
in  small  scale  model  tests.  However,  analytical  and 
numerical  techniques  to  predict  scale  model  and  prototype 
response  have  improved  dramatically  in  recent  years.  These 
techniques  have  led  to  a significant  increase  in  understand- 
ing of  dynamic  structural  response  and  of  the  way  in  which 
flaws  and  geometric  inperfections,  respectively,  influence 
fracture  and  buckling  behavior.  Although  uncertainties 
still  exist,  the  progress  made  so  far  appears  to  provide 
the  basis  for  significant  cost  and  time  savings  through 
greater  reliance  on  analysis  and  scale  models  in  structural 
design  and  evaluation.  (Sections  IV. B. 6 and  IV. C.) 

6.  Similitude  analysis  is  vital  to  the  success  of  any 
scale  model  test,  even  though  it  may  prove  impossible  to 
match  all  geometric,  kinematic  dynamic,  and  constitutive 
similitude  parameters.  Many  materials  of  DoD  interest 
exhibit  moderate  to  high  toughness  values;  however, 
similitude  analysis  and  testing  is  not  restricted  to  linear 
elastic  behavior  and  can  be  extended  to  treat  elastic-plastic 
behavior.  (Section  IV. C.) 
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B . Recommendations 


The  two  most  critical  needs  are  concerned  with  fracture 
criteria  and  scale  modeling: 

1.  Establish  criteria  and  associated  tests  for  the 
initiation,  propagation,  and  arrest  of  cracks  in  tough 
structural  materials.  Current  dynamic  fracture  mechanics 
can  treat  linear  elastic  behavior  only.  Both  experimental 
and  mathematical  analysis  work  needs  to  be  done  to  identify 
and  establish  the  appropriate  material  parameters  needed 

for  a dynamic  elastic-plastic  fracture  mechanics  methodology 
and  to  verify  the  basic  concepts.  (Section  IV. B. 4-6.) 

2.  Develop  methods  of  structural  model  testing  which 
incorporate  fracture  mechanics  principles.  In  view  of  the 
importance  of  flaws  in  the  damage  tolerance  of  structures 
subjected  to  dynamic  loading,  model  tests  must  be  able  to 
incorporate  flaws  in  sensitive  areas.  To . c ircumvent  the 

(difficulty  arising  from  thichness  dependence  of  the  fracture 
toughness,  a research  program  addressed  to  toughness- 
scaled  materials  should  be  conducted.  This  program  should 
include  analysis  and  experiment  in  order  to  select  appro- 
priate scaling  parameters  and  to  verify  that  model  tests 
using  toughness-scaled  materials  will  give  reliable  infor- 
mation on  actual  structures  with  flaws.  (Section  IV. C.) 

Additional  recommendations  follow.  (It  should  be 
noted  that  recommendations  3 and  4 are  made  with  little 
support  given  in  the  body  of  the  report  because  of  the 
paucity  of  information  reported  in  the  open  literature.) 

3.  Further  work  is  needed  to  develop  better  analysis 
techniques  to  handle  shock  wave-structure  interaction 
problems.  These  analyses  should  be  verified  by  experiment- 
model  and  prototype.  (Section  III.) 

4.  For  close-in  explosions,  such  as  contact  mines  or 
other  near  field  problems  with  high  peak  pressure,  addi- 
tional definition  of  reflected  shock  parameters  is  needed 
(Section  III.) 
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5.  Where  gaps  exist  in  data  for  dynamic  flow  stress 
of  structural  materials  of  interest,  they  should  be  filled 
and  a catalog  of  data  prepared  by  one  of  the  appropriate 
DoD  information  centers  such  as  MCIC.  (Section  IV. B. 2) 

6.  Through  experimental  and  theoretical  investigations, 
develop  improved  constitutive  equations  for  dynamic  inelastic 
behavior  under  multiaxial  stresses  and  time  varying  strain 
rates,  (Section  IV, b. 3) 

7-  Standardization  of  test  methods  for  measuring 
stress-strain  behavior  at  high  rates  should  be  considered 
before  these  data  are  used  in  critical  design  applications. 
(Section  IV. B. 3) 

8.  Develop  improved  analysis  methods  for  prediction  of 
dynamic  buckling  and  large  deformation  response  of  shell 
structures.  This  work  should  include  both  theoretical  and 
scale-model  investigations,  (Section  iv.B.6) 
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III.  LOADING 


The  starting  point  of  any  dynamic  response  analysis 
is  a definition  of  the  geometry  of  the  structure,  the 
physical  properties  of  the  materials  from  which  the  struc- 
ture is  fabricated,  and  the  nature  of  the  loading.  The 
details  of  structural  geometry  are  generally  prescribed 
and  known  within  the  accuracy  limits  of  manufacturing 
tolerances.  For  statically  loaded  metallic  structures,  the 
set  of  material  properties  needed  to  adequately  describe 
the  elastic,  inelastic  and  fracture  behavior  is  relatively 
well  defined.  In  the  dynamic  loading  regime,  there  remain 
some  uncertainties  as  to  the  adequacy  of  our  ability  to 
describe  material  behavior.  Proper  prescription  of  dynamic 
material  properties  will  be  a major  concern  of  remaining 
sections  of  this  report.  The  remaining  prerequisite,  that 
of  load  definition,  will  be  the  subject  of  this  section. 

The  committee  has  restricted  its  attention  to  a 
regime  of  high  strain-rate  loading  which  excludes  high 
velocity  projectile  or  particle  impact.  These  latter  pro- 
blems involve  localized,  intense  deformation  which  produces 
very  high  strain  rates  and  therefore  should  be  treated 
separately.  The  present  report  is  concerned  with  dynamic 
loading  sources  which  produce  gross  structural  damage,  such 
as  impingement  by  a propagating  blast  wave,  structural 
impact  or  collision,  or  the  sudden  release  of  energy  as  in 
propellant  expansion.  These  loading  sources  generally 
produce  impulse  durations  on  the  order  of  milliseconds.  The 
response  depends,  of  course,  on  the  overall  stiffness  and 
mass  distribution  of  the  structure.  Typical  inelastic 
strain  rates  associated  with  permanent  damage  are  of  the 
order  of  1-100  sec“  , with  local  rates  at  strain  concentra- 
tion sites  (notches,  corners,  cracks,  etc.)  reaching  and 
perhaps  exceeding  1000  sec“  . Our  review  of  material  con- 
stitutive properties  will  concentrate,  therefore,  on  the 
strain-rate  range  between  "static"  and  1000  sec“^ . 

In  this  section  some  typical  dynamic  loading  sources 
of  interest  to  DoD  are  described.  This  characterization  is 
intended  only  to  illustrate  the  nature  of  the  loading 
function  and  the  relative  range  of  the  load  defining  para- 
meters. The  values  presented  represent  typical  ranges  and 
therefore  are  not  recommended  for  use  in  specific  problems. 


TABLE  III-l  Ranges  of  Pressures,  Durations,  Impulse  for  Various 
Surface  Ship  Loadings 
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Land  Mine  2x10  -1.5x10  0.1-10 


Table  III-l  presents  a list  of  dynamic  loading  sources 
of  concern  to  military  structures  or  components  with  a 
typical  range  of  peak  pressures,  load  durations  and  total 
impulses  for  each  source.  Most  problems  of  the  type  we  are 
considering  involve  a single  loading  event;  in  some  cases, 
however,  the  impact  load  may  be  repetitive.  Examples  of  the 
latter  are  ship  hull  slamming  in  rough  seas  or  repetitive 
firing  of  a weapon.  In  such  cases,  damage  is  cumulative 
and  failure  may  occur  only  after  long  time  periods  (number  of 
impacts) . Figure  III-l  plots  the  peak  overpressure  versus 
duration  for  sources  listed  in  Table  III-l. 

The  details  of  the  pressure  pulses  from  various 
sources  are  illustrated  in  subsequent  figures.  For  an  air 
blast,  either  from  conventional  or  nuclear  sources,  the 
idealized  pressure  wave  is  as  given  in  Figure  III-2.  It  has 
a sharp  rise  time  followed  by  a positive  and  then  a negative 
pressure  phase.  From  a structural  point  of  view,  the 
critical  parameters  are  the  peak  positive  pressure,  the 
duration,  and  total  impulse  of  the  positive  phase.  Figures 
III-3a,  b,  and  c give  the  magnitude  of  these  parameters  for 
different  standoff  distances  for  nuclear  and  scaled  TNT 
equivalents  (Baker,  1973) . Reflected  pressures  from  the 
exposed  surface  of  a structure  in  the  field  of  a blast  wave 
may  be  many  times  the  free  field  pressure.  The  pressure- 
time history  can  be  modified  as  well  by  the  presence  of 
surrounding  reflecting  or  refracting  surfaces,  such  as  the 
ground . 

Underwater  explosion  loading  (Cole,  1965)  is  normally 
of  two  types:  plane  waves  of  moderate  peak  pressure  with 
a slow  exponential  decay  and  spherical  waves  with  rapid 
exponential  decay,  as  illustrated  in  Figure  III-4.  The 
plane  wave  is  associated  with  a nuclear  explosion  located 
at  a considerable  distance  from  the  structure  and  the 
spherical  wave  is  associated  with  a conventional  chemical 
explosion  located  close  to  the  structure.  For  the  case  of 
the  chemical  explosion,  the  loading  is  much  different  than 
the  free  field  wave,  as  shown  in  Figure  III-5  due  to  non- 
linear interaction  between  the  structure  and  the  shock  wave 
(Barton  and  Pilkey,  1974) . This  nonlinear  interaction  is 
associated  with  the  deformation  of  the  shell.  Also,  local 
cavitation  occurs  as  the  shell  responds  to  the  shock  wave. 
Typically,  a few  milliseconds  after  shock  wave  arrival  a 
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“ la  t2  + T+  ta  + T + + T' 

TIME 

FIGURE  III-2  Idealized  Blast  Wave 


SOURCE:  Baker,  W.E,,  1973 

reloading  occurs  due  to  closing  of  the  cavitation.  Other 
secondary  loads  due  to  surface  reflected  waves  or  bottom 
reflected  waves  can  also  occur. 

In  the  case  of  plane  wave  excitation  (Huang,  1969)  due 
to  an  underwater  nuclear  explosion  at  a considerable  distance 
from  the  structure,  the  response  of  a typical  cylindrical 
shell  consists  primarily  of  rigid  body  motion,  beam  whipping, 
and  axisymmetric  (breathing  mode)  motion.  Other  motions, 
such  as  asymmetric  shell  bending  modes,  contribute  to  a 
lesser  extent. 

In  terms  of  strain  response,  the  strain  rate  for  the 
nuclear  explosion  is  an  order  of  magnitude  less  than  the 
strain  rate  for  a conventional  chemical  explosion  (cf.  Table 
- The  spherical  wave  from  the  chemical  explosion 
excites  the  high  frequency  shell  bending  modes,  resulting  in 
a high  frequency  strain  response.  The  order  of  magnitude 
difference  in  strain  rates  is  based  on  both  analytical  and 
experimental  data  (Cole,  1965)  . 

For  the  case  of  mechanical  impact  on  underwater  struc- 
tures, the  structure  is  subjected  to  a high  impulse  loading 
over  a small  area  and  for  fairly  long  durations.  For  both 
mechanical  impact  and  chemical  explosions,  impulses  on  the 
order  of  1 to  10  psi-sec  are  typical  (cf . Table  III-l) . 

In  rough  seas,  the  impact  of  the  ship  hull  (Inter- 
national Ship  Structures  Congress,  1976)  against  the  free 
water  surface  can  produce  pressures  which  results  in 
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FIGURE  III-3a 


Scaled  Peak  Overpressure  vs. 
Ground  Range 


Baker,  W.E.,  1973 
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FIGURE  III-3b  Scaled  Positive  Duration  vs. 
Ground  Range 
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POSITIVE  OVERPRESSURE  IMPULSE,  psi-ms 
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FIGURE  III-3C 


Scaled  Positive  Impulse  vs 
Ground  Range 


Baker,  W.  E.,  1973 
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FIGURE  III-4. 


Free  Field  Pressure  for  Conventional 
and  Nuclear  Explosions 
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FIGURE  III-5  Shock  Wave  Pressure  History  bn 
a Cylindrical  Shell 

permanent  hull  plating  damage.  Figures  III-6a,  b,  and  c 
show  the  nature  of  this  ship  slamming  problem  for  a typical 
large  ore  carrier.  For  high-speed  surface  effects  ships, 
data  have  been  taken  for  slamming  events  and  are  reproduced 
in  Figure  III-7.  Slamming  represents  a repetitive  impact 
which  may  lead  to  damage  accumulation. 

The  detonation  of  propellant  in  guns  and  other  projec- 
tile launch  systems  is  another  source  of  impulse  loading. 
Interior  ballistics  is  concerned  with  the  pressure-time 
loading  within  the  breech  and  the  launch  tube  (barrel) . In 
some  cases,  the  transient  pressure  field  external  to  the 
muzzle,  created  as  the  projectile  exits,  may  significantly 
affect  surrounding  structures.  A particular  example  of  the 
latter  is  that  of  heavy  weapons  externally  mounted  on  light 
weight  aircraft  such  as  helicopters.  Figure  III-8  illus- 
trates typical  conditions  for  internal  ballistics  of  a 
calibre  .30  gun  (Handbook  AMCP  706-342,  1970). 

A severe  threat  to  military  land  vehicles  is  the  buried 
mine.  Severe  shock  pressures  may  be  experienced  by  tank 
hull  bottom  plating  due  to  mine  explosions.  In  this  case, 
the  total  impulse  is  produced  both  by  the  expanding  gases 
and  by  the  added  momentum  of  the  soil  mass  carried  by  the 
explosion  (Westline,  1972).  The  pressure-time  history  can 
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FIGURE  III-6a  Frequency  Distribution  of  Impact 
Pressure  on  Bow  Flare  of  Large 
Ore  Carrier 

SOURCE:  Baker,  W.  E.,  1973 


Number  of  samples 
388 


Impact  Duration  Time,  (sec.) 


FIGURE  III  6-b  Frequency  Distribution 

of  Duration  Time  of 
Impact  Pressure  on  Bow 
Flare  of  Large  Ore 
Carrier 

SOURCE:  6th  International  Ship  Structure 

Congress,  1976 
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SOURCE: 


FIGURE  III-6c  Relationship  Between 

Magnitude  and  Duration 
of  Impact  Pressure 
Observed  on  Bow  Flare 
of  Large  Ore  Carrier 

6th  International  Ship  Structure  Congress 
1976 
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SOURCE: 


FIGURE  III— 7 Manned  Testcraft  XR— 1C  Slamming 
Event  (24  Knots,  Significant 
Wave  Height  0.8  - 0,9  m.  Head 
Seas) 

6th  International  Ship  Structure  Congress, 
1976 
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FIGURE  III-8  Interior  Ballistics  of  Caliber  .30  Gun 
SOURCE:  Handbook  AMCP  706  -342,  1970 

be  very  complex  relative  to  an  air  blast,  with  the  total 
impulse  being  a function  of  charge  weight,  stand-off 
distance  and  depth  of  burial. 

The  preceding  discussion  illustrates,  by  example,  the 
nature  of  some  typical  loading  sources.  Several  points  are 
important  to  keep  in  mind.  First,  if  the  load  duration  is 
short  in  comparison  with  the  fundamental  response  time  of 
the  structure,  then  it  may  be  considered  as  an  impulse.  The 
dynamic  structural  analysis  is  then  greatly  simplified  in 
that  the  loading  can  be  treated  simply  as  an  initial  velocity 
condition  and  load-structure  interactions  are  not  important. 
In  many  cases,  however,  the  load  duration  is  comparable  to 
the  structural  response  time  (e.g.,  the  fundamental  period 
of  vibration  of  the  structure  in  the  dominant  excited  mode) . 
In  these  cases,  the  analysis  can  become  very  complex, 
including  nonlinear  interactions  between  the  loading  and  the 
deformation.  A complete  solution  in  this  case  may  require 
solving  the  transient  field  equations  in  the  surrounding 
medium  as  well  as  in  the  structure.  Adequate  analyses  for 
many  problems  of  this  type  do  not  as  yet  exist. 
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IV.  CURRENT  UNDERSTANDING  AND  APPROACHES 


A.  Introduction 


In  order  to  design  adequately  for  damage  tolerance  one 
usually  needs  a mathematical  description  of  elastic  and 
inelastic  deformation  behavior  and  associated  structural 
analysis  techniques  for  the  expected  range  of  deformation 
rates.  (As  used  here,  damage  tolerance  is  the  ability  of  a 
structure  or  component  to  absorb  mechanical  energy  without 
failure.  In  turn,  failure  is  defined  as  a loss  in  function- 
al capability  according  to  some  specified  criterion:  e.g., 
excessive  deformation  (transient  or  residual) , structural 
instability  (collapse) , or  loss  of  residual  strength  (frac- 
ture) . The  use  of  material  screening  tests  for  dynamic 
fracture  resistance  (Charpy,  dynamic-tear,  drop-weight  tear 
test,  etc.),  normally  assures  that  fracture  in  dynamically 
loaded  structures  is  accompanied  by  large  inelastic  deforma- 
tions. Therefore,  with  tough  materials  a major  need  exists 
for  realistic  models  of  rate-dependent  plastic  deformation 
and  fracture.  The  dynamic  fracture  problem  involves 
prediction  of  initiation  of  crack  growth  under  rapid  loading 
as  well  as  dynamic  crack  propagation  and  arrest.  Section  IV 
will  include  discussions  on  the  modeling  of  these  deformation 
and  fracture  phenomena. 

The  overall  problem  considered  is  that  of  describing 
the  behavior  of  materials  and  structures  when  subjected  to 
the  dynamic  loading  of  the  type  illustrated  in  Section  III. 
This  loading  is  essentially  of  short  duration  and  high 
intensity.  We  are  interested  in  both  the  transient  and  the 
residual  or  permanent  deformation  of  the  structure.  The 
permanent  deformation  also  includes  the  creation  of  new 
surface  area  if  separation  or  fracture  of  initially 
continuous  structural  elements  has  taken  place. 

Analytical  modeling  of  the  entire  process  of  deforma- 
tion and  fracture  is  a complex  process  which  is  often  broken 
up  into  subelements;  e.g.,  "elasticity,"  "plasticity,"  and 
"fracture  mechanics,"  which  are  studied  separately.  This 
introductory  discussion  will  emphasize  the  problems  of  devel- 
oping constitutive  behavior  models  for  a continuum  which 
contains  a range  of  internal  defects  of  varying  scale. 
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For  analytical  purposes materials  are  usually  modeled 
by  taking  a finite  (but  perhaps  irreducibly  small)  element 
for  which  there  can  be  defined  a deterministic  relation 
between  the  applied  surface  tractions  and  body  forces  and  the 
resulting  deformations  or  deformation  rates.  Usually,  this 
element  is  taken  to  be  homogeneous  and  continuous.  This 
homogeneous  assumption  is  certainly  satisfactory  for  problems 
in  elasticity;  however,  for  plasticity  and  for  material  fail- 
ure it  is  the  local  nonhomoqeneities  which  govern  the  behav- 
ior. Therefore,  let  us  think  of  a material  element  construc- 
ted of  the  following  components: 

1.  Ideal  crystalline  lattice. 

2.  Point  defects  - interstitials,  vacancies,  etc. 

3.  Line  defects  - dislocations. 

4.  Areal  defects  - grain  boundaries,  twin  boundaries, 

cracks . 

5.  Volume  defects  - inclusions,  pores. 

Thus,  the  volume  element  consists  of  a continuous,  homo- 
geneous crystalline  lattice  interspersed  with  a variety  of 
localized  nonhomogeneities  or  "defects."  Elastic  behavior 
is  based  upon  the  reversible  small  distortions  of  the 
crystalline  matrix.  Strength  behavior  is  governed  by  the 
much  more  complex  phenomena  associated  with  the  nucleation, 
migration,  growth,  interaction,  and  coalescence  of  the 
multiplicity  of  defects. 

In  particular  experimental  situations.,  we  may  be  able 
to  isolate  a deformation  process  where  a single  mechanism  is 
operating  alone  or  is  dominant.  For  example,  very  small 
deformation  where  only  reversible  lattice  distortions  occur 
allows  us  to  measure  the  constants  in  the  elastic  constitu- 
tive law.  A second  example  is  a single  crystal  of  a pure 
metal  or  solid  solution  alloy  where  irreversible  deforma- 
tion results  from  dislocation  glide  only.  This  is  the 
simplest  form  of  plasticity.  A specimen  with  a single  large 
crack  is  a third  example;  if  stressed  to  a level  sufficient 
to  extend  the  crack,  but  not  sufficient  to  produce  large 
scale  plastic  deformation,  this  specimen  provides  the  data 
required  for  the  application  of  linear  elastic  fracture 
mechanics . 

Generally,  however,  we  are  faced  with  developing  a 
constitutive  model  for  a structural,  polycrystalline  alloy 
which  contains  a finite  population  of  many  of  the  defects 
listed  above.  The  response  of  this  defect  structure  to  a 
given  load  environment  is  very  complex  and  depends  on  the 
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relative  population  and  distribution  of  defect  species  and 
their  modes  of  interaction.  Analytical  description  with 
relatively  good  accuracy  is  available  for  each  of  the 
isolated  individual  mechanisms:  e.g.,  interstitial  or 
vacancy  diffusion,  dislocation  glide,  twinning,  grain 
boundary  sliding,  crack  propagation,  etc.  The  challenge 
is  to  better  understand  the  nature  of  the  interaction  of 
these  mechanisms  and  thereby  develop  constitutive  and 
fracture  relations  for  real  materials  when  a multiplicity 
of  effects  are  occurring  simultaneously.  Little  progress 
has  been  made  in  this  direction. 

The  alternative  has  been  to  rely  on  stress-deformation 
response  measurements  of  bulk  specimens;  such  measurements 
effectively  integrate  all  of  the  internal  defect  motions. 

The  constitutive  relations  derived  from  this  empirical 
process  then  relate  only  the  stimulus  to  the  response  of  a 
"black-box"  material  element  (test  specimen) ; analytical 
relationships  are  fit  to  a variety  of  specific  mechanical 
test  data  corresponding  to  typical  load  or  deformation 
histories  such  as  creep,  constant  strain  rate,  fatigue, 
impact,  etc.  This  approach,  although  it  has  proven  useful, 
suffers  from  a lack  of  predictive  capacity  for  "untested" 
situations.  Improvement  must  come  from  constitutive  model- 
ing based,  as  nearly  as  possible,  on  actual  physical 
phenomena . 

Obviously,  the  matter  of  scale  is  critical  when  defin- 
ing a basic  material  element  for  the  purposes  of  developing 
constitutive  models.  The  scale  of  many  of  the  defects  which 
play  a significant,  if  not  critical,  role  in  the  constitutive 
behavior  is  of  the  same  dimension  as  the  lattice  parameter 
(spacing  between  nearest  neighbor  atoms) . Therefore,  unless 
we  reduce  the  problem  to  a discrete  lattice  model  (individ- 
ual atoms  bound  together  by  pair-wise  potentials) , our 
minimum  material  element  must  contain  a distribution  of 
defects  (inhomogeneities) . 

Finally,  it  is  important  to  consider  the  nonhomogenei- 
ties to  be  structural  defects,  having  not  only  specific 
geometrical  configurations  but  also  to  have  associated  with 
them  local  residual  stress  fields  produced  by  distortions 
of  the  surrounding  matrix,  usually  assured  to  be  an  elastic 
continuum.  In  addition,  the  deformation-induced  stress 
field  may  be  further  intensified  by  the  presence  of  other 
local  structural  inhomogeneities  (inclusions,  porosity. 
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ect.).  The  development  and  interaction  of  these  small-scale 
inhomogenities  are  the  factors  which  control  both  plastic 
deformation  and  fracture.. 

B . Analytical  Modeling  and  Related  Materials  Property  Tests 

1.  Basic  Deformation  Mechanisms 

Let  us  consider  those  processes  in  metals  which  are 
dependent  upon  the  plastic  strain  rate.  In  considering  the 
underlying  causes  of  plasticity,  one  normally  deals  with  the 
effect  of  stress  and  temperature  on  the  initiation  and 
propagation  of  dislocations.  In  fracture,  one  is  concerned 
at  least  with  the  effect  of  the  same  variables  on  the 
initiation  and  propagation  of  cracks  in  brittle  materials 
or  of  voids  and  ductile  materials.  If  the  problem  is  one 
of  solely  elastic  deformation,  structural  analysis  is 
greatly  simplified  because  the  material  properties  are  by 
definition  independent  of  the  time-dependence  of  the  deforma- 
tion. Inelastic  deformation  in  metals  is  primarily  the 
result  of  the  motion  of  dislocations  and  their  interactions 
with  other  internal  defects  in  the  lattice.  The  rate 
controlling  mechanisms  for  this  process  are  outlined  in  this 
section. 

Table  IV-1  and  Figures  IV-1  and  IV  -2  summarize  the 
known  plastic  deformation  rate  controlling,  mechanisms  which 
apply  to  metals  over  the  entire  ranee  of  temperature  and 
stress.  The  four  basic  mechanisms  are  viscous  drag,  thermal 
activation,  athermal  stress  barriers,  and  diffusion.  Several 
specific  examples  of  each  mechanism  are  given  in  Table  IV-1. 
Figure  IV-1,  which  is  a modified  version  of  deformation  map 
as  used  by  Ashby,  (1972)  illustrates  curves  o'f  constant 
strain-rate  in  the  stress  x - temperature  (T)  plane.  Stress 
and  temperature  are  normalized  by  the  shear  modulus,  jj,,  and 
melting  temperature,  T^j,  respectively.  The  curves  are  for 
constant  values  of  dislocation  velocity,  v,  normalized  by 
the  shear  wave  velocity,  Cs-  This  parameter,  v/Cg,  can  be 
related  to  plastic  shear  strain  rate,  y , by  the  kinematic 
relation: 


y = pbv  (IV-1) 

where  P is  the  density  of  mobile  dislocations  and  b is  the 
Burger's  vector  of  the  dislocations.  In  Figure  IV-1,  the 
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TABLE  IV-1  Basic  Deformation  Mechanisms 


1.  Viscous  Drag  on  Dislocations 

• Thermoelastic 

• Phonon  scattering 

• Phonon  viscosity 

• Electron  viscosity 

2.  Thermal  Activation  of  Dislocations 
Past  Short-Range  Barriers 

• Dislocation  intersection 

• Solute  atoms 

• Peierls  mechanisms 

• Cross  slip 

3.  Athermal  Stress  Fields 

• Long-range  stress  fields 

• Incoherent  precipitates 

4.  Diffusion  Controlled  Dislocation  Motion 

• Dislocation  climb 

• Nabarro-Herring  creep 

• Coble  creep 


dimensional  scales  are  roughly  appropriate  for  aluminum 
where  the  dislocation  density  is  taken  as  10^°cm“^  . Figure 
IV-2  shows  curves  of  constant  homologous  temperature  for  the 
same  material. 

The  four  basic  mechanisms  can  be  separated  roughly  by 
stress  level.  At  the  highest  stresses,  t > Tq,  dislocations 
are  not  impeded  by  other  defects  and  the  only  drag  forces 
are  those  generated  by  the  interaction  of  the  high  velocity 
dislocations  with  the  basic  crystalline  lattice.  This 
interaction  can  result  in  several  dissipative  mechanisms  as 


31 


NORIVIALIZED  STRESS  tI\i 


TEMPERATURE  T (*K) 


FIGURE  IV-1  Deformation  Map.  Dimensional  Parameters 
Correspond  Roughly  to  Aluminum 


SOURCE;  Adapted  from  Ashhy,  1972 
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STRESS  T<  MN/m' 


NORMALIZED  Sil?ESS,  T/^ 


STRAIN  RATE-y  (s'^) 


FIGURE  IV-2  Deformation  Map  in  Terms  of  Dislocation  Velocity 
or  Strain  Rate.  Alternate  Presentation  of  Fig- 
ure IV-1 

SOURCE:  Adapted  from  Ashby,  1972 
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listed  in  Table  IV-1,  all  of  which  yield  a linear  vicous 
relation  between  dislocation  velocity  (or  strain  rate)  and 
applied  shear  stress;  viz., 


bT  = Bv 


(IV-2) 


where  B = B (T)  is  a net  drag  coefficient.  The  magnitude  of 
this  coefficient  is  such  that  only  the  viscous  force  is 
significant  at  very  high  dislocation  velocities. 

For  stresses  in  the  range  < t < , the  dominant 

rate  controlling  mechanism  is  thermal  activation  of  dis- 
locations past  short-range  obstacles , shown  schematically  in 
Figure  IV-3.  In  this  range,  the  governing  equation  is  of 
the  Arhennius  form 


y 


B ' exp 


U(T) 

kT 


(IV-3) 


where  the  coefficient  B^  of  the  exponential  is  often  taken 
as  a constant  but  must  contain  the  number  of  activation 
sites,  area  swept  out  by  a successful  activation  event,  the 
Burger's  vector,  and  the  attempt  frequency.  U = U(T)  is  the 
stress-dependent  thermal  activation  energy  and  k is  the 
Boltzmann  constant. 

Thermal  activation  is  effective  for  only  short  range 
stress  fields  with  dimensions  of  the  order- of  the  lattice 
spacing.  Movement  of  dislocations  through  internal  stress 
fields  of  longer  range  can  be  accomplished  only  by  supplying 
the  energy  by  means  of  the  applied  stress.  The  athermal 
stress,  Tg , is  a mean  stress  required  to  overcome  such  long 
range  internal  barriers. 

The  total  picture  given  above  is  useful  conceptually, 
but  it  should  be  kept  in  mind  that  within  a given  range  of 
parameters  a number  of  distinct  defect  interactions  are 
probable,  only  a few  of  which  are  listed  in  Table  IV-1. 
Therefore,  in  reality,  especially  for  commercial  alloys,  the 
situation  is  more  complex  than  shown.  Obviously,  a simple 
constitutive  equation  cannot  account  for  all  possible 
internal  defect  interactions  explicitly.  Therefore,  for  each 
alloy  we  must  have  some  guidelines  as  to  the  dominant 
controlling  mechanisms  for  a given  set  of  the  independent 
variables.  This  is  the  intent  of  the  Ashby  type  deformation 
map.  In  the  same  sense,  fracture  maps  may  be  constructed 
also  (Ashby,  1977) . 


34 


FIGURE  IV-3  Thermal  Activation  Past  Square  Array 
of  Localized  Obstacles 


a.  stress  vs.  displacement 

b.  successive  positions  of  dislocation 
in  plan  of  motion  with  regular  array 
of  obstacles 
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2 . Metallurgical  Aspects  of  Strain  Rate  Sensitivity 

As  discussed  above,  the  stress  associated  with  plastic 
flow  of  a material  is  customarily  represented  as  a function 
of  two  external  variables,  the  instantaneous  plastic  strain 
rate  and  the  temperature  (cf . Figure  IV-2) . The  degree  of 
sensitivity  to  these  variables  depends  very  much  upon  the 
material  as  well  as  its  thermo-mechanical  history.  For 
example,  the  variation  in  this  "flow  stress"  with  either 
plastic  strain  rate  or  temperature  is  much  more  pronounced 
in  low  strength  as  compared  to  high  strength  steels,  and 
the  dependence  upon  these  variables  is  much  greater  in  body- 
centered  cubic  metals  (e.g.,  iron)  than  in  face-centered 
cubic  metals  (e.g.,  aluminum).  There  are  important  practi- 
cal ramifications  to  this  dependence  of  the  flow  stress  on 
these  variables,  as  in  high  speed  manufacturing  operations 
where  the  flow  stress  at  the  strain  rate  involved  rather 
than  the  static  strength  determines  the  requisite  machine 
capacity.  In  addition,  the  fracture  process  in  steels,  one 
matter  which  will  be  of  concern  herein,  can  be  deleterious- 
ly  influenced  by  either  an  increase  in  strain  rate  or  a 
decrease  in  temperature.  This  review  will  be  primarily 
directed  at  strain-rate  sensitivity,  which  is  almost  always 
coupled  to  temperature  sensitivity,  and  will  focus  on  the 
controlling  factors  giving  rise  to  a dependence  of  the  flow 
stress  of  an  alloy  on  these  parameters.  In  addition,  the 
influence  of  strain  rate  and  temperature  on  fracture  be- 
havior of  precracked  alloys  will  also  be  considered.  In 
order  to  set  some  bounds  on  the  range  of  the  variables  con- 
sidered, the  temperatures  of  interest  will  be  well  below 
th^  creep  range,  and  an  upper  limit  in  strain  rate  of 
10  /sec  will  be  set  so  that  consideration  of  shock  phenom- 
ena which  can  involve  strain  rates  of  up  to  10^/sec  will  be 
exc luded . 

In  the  strain  rate  range  of  interest ; i . e . , below 
10^/sec,  the  principal  factors  leading  to  a strain  rate 
effect  are  temperature  and  the  time  available  for  thermal 
activation  of  dislocation  motion;  for  aluminum  see  Figure 
IV-2.  The  sensitivity  of  the  fracture  process  to  strain 
rate  and  temperature  results  directly  from  the  fact  that 
the  fracture  toughness  of  an  alloy  depends  upon  the  extent 
of  plastic  deformation  involved  in  the  rupture  process; 
normally,  the  greater  the  amount  of  localized  plastic  de- 
formation, the  tougher  the  material,  and  any  factor  which 
inhibits  plastic  deformation  can  be  expected  to  also  lower 
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the  fracture  toughness.  Further,  in  materials  that  can 
fail  by  cleavage,  below  a critical  temperature  an  increase 
in  strain  rate  will  promote  brittle  behavior  if  the  local 
stresses  at  a notch  or  crack  rise  rapidly  enough  to  satisfy 
a fracture  criterion  before  sufficient  time  has  elapsed  to 
permit  the  thermal  activation  of  enough  plastic  deformation 
to  relax  these  high  stresses.  For  such  materials,  whenever 
fracture  under  dynamic  loading  conditions  is  of  concern  it 
is  important  that  the  factors  responsible  for  strain-rate 
sensitivity  be  understood,  and  the  magnitude  of  related 
strain-rate  effects  be  established. 

Before  going  into  a more  detailed  consideration  of  the 
underlying  causes  of  strain-rate  sensitivity,  a few  examples 
will  be  given  to  illustrate  the  extent  to  which  strain  rate 
can  influence  the  flow  stress.  In  Figure  IV-4  a comparison 


FIGURE  IV-4  Stress-Strain  Curves  of  Low  Carbon 
Steel  Under  Static  and  Dynamic 
Loading.  Times  in  Microseconds  are 
Shown  on  Dynamic  Curve 

SOURCE:  Campbell  & Harding,  1960 
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is  given  of  flow  stress  of  a "low  carbon  steel  under  condi- 
tions of  static  (=“10“^/sec)  and  dynamic  (“^2  x 10^/sec) 
tensile  straining  conditions  (Campbell  and  Harding,  1961) . 
In  this  material,  plastic  deformation  initially  involves 
the  spread,  of  a Luders  band,  and  the  increase  in  the  upper 
yield  point  with  strain  rate  comes  about  because  disloca- 
tions are  initially  locked  by  interstitial  atoms  and  a 
finite  time  is  required  to  initiate  dislocation  motion  and 
the  spread  of  this  Luders  band.  The  dip  in  the  dynamic 
stress-strain  curve  is  related  to  the  spread  of  this  band. 
Beyond  the  yield  region,  the  flow  stress  initially  rises 
and  then  slowly  falls  off  and  approaches  the  static  curve 
at  higher  strains.  This  falloff  in  strain-rate  sensitivity 
with  increasing  strain,  but  without  a decrease  in  the  flow 
stress  itself,  also  occurs  in  the  case  of  iron  tested  in 
compression  as  shown  in  Figure  IV-5  (Hockett  and  Zukas, 

1973) . With  high-purity  aluminum.  Figure  IV-6,  at  least 
over  the  strain  range  tested,  this  trend  is  not  followed 
(Hauser,  et  al.,  1961). 


FIGURE  IV-5  True  Stress  Against  True  Strain 

C'^urves  for  Iron  Compressed  at  2 95  K 
at  the  Indicated  Constant  True 
Strain  Rates 

SOURCE:  Hockett  and  Zukas,  1973 
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FIGURE  IV-6  Effect  of  Strain  on  Stress  at 
Constant  Strain  Rate 

SOURCE:  Hauser  et  al,  1961 

In  Figure  IV-7,  the  dependency  of  the  lower  yield  of 
a low  carbon  steel  is  shown  as  a function  of  grain  size, 
radiation  dosage,  and  impact  velocity  (Campbell  and  Harding, 
1961) . The  stress  given  by  the  intersection  of  straight 
lines  with  the  ordinate  is  known  as  the  friction  stress  and 
it  is  seen  that  the  large  increase  in  flow  stress  with  in- 
creasing velocity  is  due  to  an  increase  in  the  thermally 
activated  friction  stress.  The  grain  size  dependency  is 
independent  of  strain  rate,  and  hence  independent  of 
thermal  activation,  since  it  is  due  to  a long-range  stress 
associated  with  grain  boundaries.  Additional  examples  of 
the  influence  of  strain  rate  and  temperature  on  the  lower 
yield  stress  of  a 0.2  weight  percent  carbon  steel  are  given 
in  Figure  IV-8  and  IV-9  (Taylor  and  Malvern,  1961) . The 
strong  dependency  of  the  lower  yield  stress  on  those  test 
variables  is  evident. 
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SOURCE : 


FIGURE  IV-7  Lower  Yield  Stress  of 

Steel:  (a)  unirradiated ; 

(b)  radiation  dose  8 x 10 
neutrons/sq.  cm;  (c)  radiation 
dose  5 X 10^®neutrons/sq.  cm. 

A = static  loading 

o = impact  velocity  198  ips 

X = impact  velocity  507  ips 

Campbell  and  Harding,  1961 


40 


o 


o 


SOURCE : 
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FIGURE  IV-8  Dynamic  Lower  Yield  Stress 
as  a Function  of  Specimen 
Strain  Rate  for  Normalized 
(Present  Investigation) 
and  Unnormalized  (Taylor) 
Mild  Steel  at  Room 
Temperature.  (0.26  steel) 

Taylor  and  Malvern,  1961 
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FIGURE  IV-9  Dynamic  Lower  Yield  Stress 
as  a Function  of  Test 
Temperature  (.2c  steel) 

SOURCE:  Taylor  and  Malvern,  1961 

In  order  to  understand  the  basis  for  strain-rate 
sensitivity,  it  is  useful  to  consider  again  the  interaction 
of  a dislocation  with  the  various  types  of  barriers  in  a 
material  which  oppose  the  motion  of  the  dislocation  under 
stress.  These  barriers  to  dislocation  motion  can  be  divided 
into  two  groups  according  to  the  distance  over  which  they 
interact  with  glide  dislocations  (Basinski,  1959,  and 
Basinski  and  Christian,  1960) . In  face-centered  cubic  metals 
such  as  copper  and  aluminum,  the  principal  obstacles  exert- 
ing short-range  forces,  such  as  forest  dislocations  or  jogs, 
can  be_  overcome  at  finite' temperatures  with  the  help  of 
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thermal  fluctuations,  even  though  the  force  on  the  glide  dis 
locations  is  smaller  than  the  force  exerted  by  the  obstacles 
This  type  of  interaction  will  lead  to  a temperature  and 
strain-rate  dependent  flow  stress.  On  the  other  hand,  the 
activation  energies  involved  when  the  dislocation  moves 
under  adverse  stress  over  a large  distance  would  be  outside 
the  range  of  thermal  fluctuations  at  ordinary  temperatures. 
We  must  expect,  therefore,  that  the  flow  stress  originating 
from  long-range  elastic  interactions  between  dislocations 
should  be  independent  of  the  strain  rate,  and  depend  upon 
the  temperature  only  through  the  elastic  moduli. 

Body-centered  cubic  metals  such  as  iron  differ  from 
face-centered  cubic  metals  in  that  there  is  an  additional 
important  type  of  barrier  known  as  the  Peierls-Nabarro 
stress  which  results  from  an  interaction  of  a dislocation 
with  the  lattice.  This  barrier  is  of  a short-range  type 
and  can  be  overcome  with  the  aid  of  thermal  fluctuations. 
Additionally,  dislocation  interaction  with  interstitial 
atoms  such  as  carbon  and  nitrogen  can  occur,  but  it  appears 
that  the  Peierls-Nabarro  stress  rather  than  an  interstitial 
effect  gives  rise  to  significant  friction  stress  and  is 
responsible  for  the  high  degree  of  temperature  and  rate 
sensitivity  in  bcc  metals. 

The  following  type  of  equation  was  proposed  by  Seeger 
(1957)  to  account  for  the  temperature  and  strain  rate  sen- 
sitivity of  the  flow  stress  in  the  presence  of  a particular 
short-range  barrier  such  as  the  cutting  of  a glide 
dislocation  by  forest  dislocations.  An  applied  shear  stress 
T,  can  aid  in  the  motion  of  the  dislocation  through  this 
type  of  barrier;  however,  a portion  of  the  applied  stress, 
Tg,  may  be  needed  to  overcome  long-range  barriers  due  to  the 
stress  fields  of  parallel  dislocations.  In  this  model,  the 
number  of  activated  sites  is  independent  of  temperature. 

The  plastic  strain  rate,Y  » is  expressed  as 


Y = B'  exp 


XJ  —V  (t— T ) 

icT 


(IV-4) 


where  B=bA'NvQ,  and  b is  the  Burgers  vector.  A'  is  the  area 
swept  out,  N is  the  number  of  activated  sites,  Vq  is  a 
frequency  factor  which  depends  upon  the  obstacle  and  the 
way  it  is  overcome,  k is  Boltzmann's  constant,  is  the 
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constant  activation  energy,  and  is  known  as  the  activa- 
tion volume.  From  Equation  IV-4,  the  following  expression 
can  be  derived  from  the  shear  stress  for  flow. 


Uq  - kTln(Na'bt^o/'y)  : 

' 'g  

V_ 


T 


where  Tq= 


Up 

kin  (Na ' huo/y) 


T^Tq 


(IV-5a) 


(lV-5b) 


The  dependency  of  the  flow  stress  upon  temperature  and 
strain  rate  based  on  this  equation  is  shown  in  Figure  IV- 10, 
which  depicts  a limited  region  of  interest  of  Figure  IV-1. 
The  linear  dependency  of  the  flow  stress  on  temperature  at 


FIGURE  IV-10  Dependence  of  the  Flow^Stress  on  Temperature 
as  Predicted  by  Equation  IV-5 


44 


a given  strain  rate  as  predicted  by  Equation  IV-5a  is  in 
general  agreement  with  observations  on  the  hexagonal  of 
trigonal  metals;  e.g.,  magnesium,  zinc,  cadmium,  and  bis- 
muth, and  on  the  face-centered  cubic  metal  aluminum 
(Seeger,  1957) . In  the  case  of  titanium,  a hexagonal  metal, 
the  variation  is  not  linear  but  the  thermal  and  athermal 
regions  are  apparent  as  shown  in  Figure  IV-11  (deMeester, 
et  al.,  1975).  It  is  clear  from  Equation  IV-5a  that  tem- 
perature and  strain  rate  effects  are  interrelatable,  and 
the  following  parametric  relationships,  derivable  from 
Equation  IV-5a  have  been  used  to  correlate  the  uniaxial 
strain  rate  c and  temperature  sensitivity  of  the  uniaxial 


FIGURE  IV-11  The  0.2  Percent  Yield  Stress  Divided  by  the 
Shear  Modulus  Versus  Temperature.  Also 
included  are  values  of  the  athermal  component 
of  the  flow  stress  divided  by  the  modulus,  as 
determined  by  various  methods 


SOURCE:  deMeester  et  al,  1975 

flow  stress,  ct,  with  a and  e equated  to  2t  and  2y/3,  res- 
pectively; 

cr  = fl  (Tin  A/e)  (IV-6a) 

or,  equivalently, 

0^  = ^2^'^m)  where  T^  s (1  - Q/ln  I)  t (IV-6b) 

^ 0 
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where  A,  a,  and  Cq  constants,  and  Tj  is  the  velocity- 

modified  temperature  of  MacGregor  and  Fisher  (1945). 

Figure  IV-12  is  an  example  of  the  use  of  Equation  IV-6a 
(Bennett  and  Sinclair,  1966)  . Figures  IV-13  and  IV- 14  are 
illustrations  of  the  use  of  the  velocity-modified  tempera- 
ture approach  (MacGregor  and  Fisher,  1945) . 


FIGURE  IV-12  Yield  Strength  for  Seven  Steels 
in  Terms  of  the  Rate-Temperature 
Parameter,  TlnA/c.  A = 10°sec 

SOURCE:  Bennett  and  Sinclair,  1966 
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FIGURE  IV-13  Stress  as  a Function  of  True  Strain 

and  Velocity-Modified  Temperature  for 
Annealed  1020  and  1045  Steels.  Strain 

C 

rates  from  5 to  500  x 10“  /sec  are 
correlated  on  each  curve 

SOURCE:  MacGregor  and  Fisher,  1945 
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FIGURE  IV-14  Stress  as  a Function  of  True  Strain 
and  Velocity -Modified  Temperature 
for  Annealed  Brass 

SOURCE:  MacGregor  and  Fisher,  1945 


Basinski  (1959)  has  used  a modified  form  of  Equation 
IV-4  in  the  analysis  of  the  temperature  and  strain-rate 
sensitivity  of  face-centered  cubic  metals.  This  analysis 
allows  for  an  increase  in  the  number  of  activated  sites 
with  an  increase  in  temperature  and  also  includes  consid- 
eration of  the  decrease  of  the  elastic  moduli  with  tempera- 
ture. Again  dislocation  intersections  are  considered  to  be 
the  important  barriers.  Basinski  and  Christian  (1960)  have 
also  analyzed  the  behavior  of  decarburized  iron  and  iron 
containing  a small  amount  of  carbon  and  nitrogen  (0.0027% 
carbon  and  0.004%  nitrogen).  The  principal  barrier  to  dis- 
location motion  in  these  materials  was  found  to  be  the 
Peierls-Nabarro  stress,  and  the  equation  for  the  strain 
rate  was  similar  to  that  employed  by  Seeger,  Equation  IV-4, 
but  with  a modified  dependency  of  the  activation  energy  on 
stress.  The  variation  of  the  flow  stress  with  temperature 
for  iron  is  shown  in  Figure  IV-15a.  Basinski 's  analysis 
which  is  confirmed  by  the  data  leads-  to  the  conclusion  that 
the  strain-rate  sensitivity  (A?)  should  varv  as 

T,  f 
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FIGURE  IV" 15a  Temperature  Dependence  of 

Reversible  Flow  Stress  in  Iron. 

Represent  Points  Taken  at 
2 Percent  Strain;  O,  Represent 
Points  Taken  at  10  Percent  Strain; 

□ Represent  Points  Taken  at  20  Per- 
cent Strain.  e = 10"^/sec 

SOURCE:  Basinski  and  Christian,  1960 

(|f)T,|  = 

where  C and  D are  independent  of  temperature.  The  strain- 
rate  sensitivity  should  therefore  rise  at  low  temperatures 
and  then  fall  off  at  higher  temperatures  as  a critical 
temperature  is  approached.  This  behavior  is  shown  in 
Figure  IV-15b.  We  note  that  the  degree  of  strain-rate 
sensitivity  in  iron  is  much  greater  than  for  aluminum  and 


49 


Temperature,  °K 

FIGURE  iV-15b  Strain  Rate  Sensitivity  of  the  Flow 
Stress  of  Polycrystalline  Iron 

SOURCE:  Basinski  and  Christian,  1960 
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other  face-centered  cubic  metals  also  studied  by  Basinski. 
Also,  the  fact  that  the  strain-rate  sensitivity  of  decarbu- 
rized  iron  is  quite  similar  to  that  of  non-decarburized  iron 
provides  further  evidence  of  the  importance  of  the  Peierls- 
Nabarro  stress  rather  than  some  interstitial  effect  in 
controlling  the  strain-rate  sensitivity  (Basinski  and 
Christian,  1960) . 

Next  let  us  consider  results  of  Davies  and  Magee  (1975) 
from  constant  strain  rate  tests  carried  out  at  room  temper- 
ature on  a variety  of  steels  and  aluminum  alloys.  These 
results  are  shown  in  Figures  IV-16  through  IV-21.  At  strain 


figure  IV-16  Ultimate  Tensile  Strength  as 
a Function  of  Strain  Rate  for 
the  -Series  of  Hot-Rolled  Steels 
Containing  0.10  to  0.82  Percent 
Carbon 

SOURCE:  Davies  and  Magee,  1975 
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SOURCE: 


SOURCE: 


FIGURE  IV-17  Ultimate  Tensile  Strength  as 
a Function  of  Strain  Rate  for 
Annealed  and  Hard  Rolled  1010 
Steel,  Showing  Reduction  in 
Strain  Rate  Sensitivity  with 
Prior  Cold  Work 


Davies  and  Magee,  1975 


FIGURE  IV-18  The  Yield  and  Ultimate 

Tensile  Strength  as  a Function 

of  Strain  Rate  for  4330  Steel 

that  was  Quenched  and  Tempered 

to  R 37 
c 

Davies  and  Magee,  1975 
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SOURCE : 


SOURCE: 


FIGURE  IV-19  Yield  and  Ultimate  Tensile 
Strength  as  a Function  of 
Strain  Rate 


Davies  and  Magee,  1975 


FIGURE  IV-20  Yield  and  Ultimate  Tensile 
Strength  as  a Function  of 
Strain  Rate  for  302  and  310 
Stainless  Steel 

Davies  and  Magee,  1975 
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FIGURE  IV-21  Yield  and  Ultimate  Tensile 
Strength  as  a Function  of 
Strain  Rate  for  the  Fully 
Hardened  Aluminum  Alloys 
7075  and  6061 

SOURCE:  Davies  and  Magee,  1975 

rates  in  excess  of  80/sec  reflection  of  elastic  waves 
("ringing")  occurred  which  interfered  with  the  determina- 
tion of  the  flow  stress  at  low  strains,  but  beyond  100  per-- 
cent  elongation  this  effect  died  out  and  therefore  the 
ultimate  tensile  strength  rather  than  the  yield  strength 
is  plotted  in  some  of  these  figures.  For  comparison,  data 
for  Ti-6A1-4V  are  also  shown  (Hatch,  1965)  as  are  the  re- 
sults of  MacGregor  and  Fisher  (1945)  from  Figure  IV-13  for 
the  flow  stress  of  1020  steel. 

From  equation  IV-5a,  at  constant  temperature  uniaxial 
stress  cr  can  be  expressed  as 

c-  = CTg  + Cl  Ine/co  (IV-8a) 


which  predicts  that  the  flow  stress  should  vary  linearly 
with  the  log  of  the  strain  rate;  this  behavior  seems  to  be 
true  for  aluminum  alloy  (cf . Figure  IV-21) , but  not  true 
for  the  steels  (cf . Figures  IV-16  through  IV-20) . 

If  both  sides  of  Equation  IV-8a  are  divided  through 
ap , then 

a C]_ 

= 1 + — In  f/c-  (IV-8b) 

rr  rr 
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or 


R = 1 + In  (IV-8c) 

where  R = cr/ag  , K = Ci/a^  . The  quantities  R and  are 
thus  indices  of  strain-rate  sensitivity.  Values  for  both  R 
and  from  Figures  IV-16  through  IV-21  are  given  in  Table 
IV-2.  It  can  be  seen  that  the  value  of  R decreases  with 
increasing  strength  level.  Note  that  the  strain-rate  sensi- 
tivity of  the  titanium  alloy  appears  to  be  greater  than  that 
of  a steel  of  comparable  strength  level. 

With  respect  to  the  influence  of  strain  rate  on  frac- 
ture behavior  of  precracked  steel  specimens,  it  is  of  inter- 
est to  note  that  the  trends  in  fracture  behavior  have  been 
related  to  the  same  parameter  T In  A/f  (Eqgation  IV-6a) 
which  has  been  used  to  correlate  the  strain  rate  and  tem- 
perature dependence  of  the  flow  stress,  Figure  IV-12.  The 
fracture  toughness  (cf.  Section  IV-4)  of  several  steels 

as  a function  of  this  parameter  is  shown  in  Figures  IV-22a,b, 
and  c (Hatch,  1965).  The  degree  of  correlation  is  quite  good 


FIGURE  iy-22a  Crack  Toughness  as  a 'Function  of 

the  Parameter  for  E Steel,  A = 10  /sec 

SOURCE;  Corten  and  Shoemaker,  1967 
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TABLE  IV-2  Dynamic  Factors  for  Several  Alloys 
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SOURCE:  Davies  and  Magee,  1975 
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FIGURE  IV-22b  Crack  Toughness  as  a Function  of  the 


Parameter  for  A302B  Steel,  A = 10  /sec 
SOURCE:  Corten  and  Shoemaker,  1967 


5: 


FIGURE  IV-22c  Crack  Toughness  as  a Function  of 

the  Parameter  for  ABS  Class  C Steel 
A = 10^/sec 

SOURCE:  Corten  and  Shoemaker,  1967 
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for  two  of  the  steels,  but  the  extent  of  the  scatter  in  the 
data  for  the  third  steel  does  not  enable  a conclusion  to  be 
drawn  as  to  the  degree  of  correlation.  (In  this  analysis 
the  strain  rates  used  were  correlated  with  a point  within 
the  plastic  zone  using  an  assumed  elastic  stress  distribu- 
tion. The  strain  rate  for  other  points  such  as  on  the 
elastic-plastic  boundary  were  also  used  without  any  signifi- 
cant shift  in  data  correlation.)  A good  correlation  implies 
that  there  is  a strong  interrelationship  between  the  strain 
rate  and  temperature  sensitivity  of  the  flow  stress  and 
fracture  toughness. 

On  a plot  of  the  fracture  toughness  or  impact  energy 
absorbed  versus  temperature,  an  increase  in  the  rate  of 
loading  often  leads  to  a shift  of  the  toughness-temperature 
curve  to  lower  temperature,  as  illustrated  in  Figure  IV-23; 


FIGURE  IV-23  Slow-Bend  and  Impact  CVN  Test 

Results  for  ABS-C  Steel, 

O’  = 39  ksi 
Y 

SOURCE:  Rolfe  and  Barsom,  1977 
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no  significant  shift  exists  for  the  high-strength  steel  of 
Figure  IV— 24.  Barsom  (1975)  has  developed  a useful  empiri- 
cal relationship  to  account  for  the  magnitude  of  this  shift 
in  °F: 


= 215  - 1.5  a for  36  < a < 140  ksi  (IV-9) 

shift  ys  y ' 

and 

T , . = 0.0  for  cr  > 140  ksi 

shift  y 

Referring  to  Figure  IV-23,  this  shift  is  the  horizontal  dis 
tance  between  the  intersection  of  the  dashed  lines  (tangent 
lines  drawn  from  the  lower  shelf  level  aid  the  transition 
region  for  the  impact  test)  to  the  onset  of  the  transition 
on  the  static  curve.  This  relationship  was  originally 
developed  for  Charpy  V-notch  specimens  but  as  Figures  IV-25 
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FIGURE  IV-24  Slow-Bend  and  Impact  CVN  Test 
Results  for  18Ni(180)  Steel 
= 180  ksi 

SOURCE:  Rolfe  and  Barsom,  1977 
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FIGURE  IV-25  Use  of  CVN  Test  Results  to 

Predict  the  Effect  of  Loading 
Rate  on  for  ABS-C  Steel. 

CTy  = 50  ksi 

SOURCE:  Rolfe  and  Barsoiti,  1977 

IV-26  (Rolfe  and  Barsom,  1977)  indicate  -it  has  also  been 
applied  to  fracture  toughness  tests.  We  again  note  that 
strain  rate  effects  are  much  more  pronounced  for  low- 
strength  steels  as  compared  to  high-strength  steels.  In  the 
ductile  region  above  the  transition  range  it  is  believed 
that  the  increase  in  flow  stress  with  strain  rate  produces 
the  increase  in  toughness  of  these  alloys  (cf.  Figure  IV-23) 
In  conclusion,  it  appears  that  a reasonable  under- 
standing of  the  processes  governing  strain  rate  effects  on 
flow  stress  has  been  achieved,  and  that  experimental  results 
are  available  for  a range  of  engineering  materials  of  inter- 
est. The  situation  with  respect  to  a basic  understanding 
of  the  dependence  of  the  fracture  process  on  strain  rate  is 
not  as  well  established  as  for  flow  behavior.  However,  it 
is  encouraging  that  in  certain  cases  the  trends  in  fracture 
behavior  with  temperature  and  strain  rate  have  been  related 
to  the  same  parameter  (T  In  A/c)  which  has  been  used  to  cor- 
relate the  flow  stress;  also,  experimental  results  are  avail 
able  to  assist  in  the  assessment  of  the  effects  of  strain 
rate  in  particular  instances. 
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3 . Dynamic  Test  Techniques  and  Constitutive  Equations 

The  previous  section  has  dealt  with  metallurgical 
related  phenomena  which  result  in  plastic  strain  rate 
sensitivity.  This  section  will  review  briefly  the  test 
techniques  which  allow  us  to  measure  quantitatively  the 
strain-rate  effect  in  bulk  specimens  and  the  continuum 
constitutive  equations  which  are  used  to  include  strain-rate 
sensitivity  in  dynamic  structural  analysis. 

An  extensive  review  and  bibliography  of  high  strain- 
rate  test  techniques  has  been  given  by  Lindholm  (1971) . The 
early  work  in  this  area,  which  began  during  and  after  World 
War  II,  initially  emphasized  measurements  of  large  amplitude 
plastic  waves  in  bars.  Plastic  wave  profiles  (strain  or 
particle  velocity)  were  measured  as  functions  of  time  or 
position  in  the  bar.  Derivation  of  the  constitutive 
properties  of  the  bar  material  from  such  measurements 
required  the  solution  of  the  transient  initial  and  boundary 
value  problem  using  an  a priori  assumed  form  of  the  consti- 
tutive relation.  This  is  a consequence  of  the  fact  that 
stress  is  not  measured  directly  and  independently  in  the 
experiment.  Adequacy  of  the  assumed  form  of  the  constitutive 
relation  must  then  be  based  upon  the  degree  of  agreement 
between  the  analysis  and  other  experimental  information  on 
wave  profiles  not  used  as  part  of  the  material  characteriza- 
tion procedure.  Two  critical  problems  arise  with  this 
method.  First,  it  has  been  found  that  the  plastic  wave 
profiles  are  not  strongly  sensitive  to  the  explicit  form 
of  the  rate-dependent  terms  in  the  constitutive  relation. 
Thus,  the  experiment  is  not  very  discriminating  in  terms  of 
the  material  property  of  interest.  Secondly,  as  the  wave 
propagates,  the  strain  rate  can  vary  widely;  however,  a 
constant  strain-rate  test  is  normally  desired  in  order  to 
minimize  strain-history  effects.  These  factors,  and  others, 
led  to  considerable  uncertainty  and  some  confusion  during  the 
50s  and  60s  regarding  the  importance  of  strain  rate  in 
dynamic  plasticity. 

A significant  advance  in  dynamic  testing  was  made  when 
Kolsky  (1949)  developed  the  Hopkinson  pressure  bar  technique 
for  compressdon  testing.  In  this  technique  a short  compres- 
sion specimen  is  sandwiched  between  two  long,  elastic  bars 
which  serve  to  transmit  the  loading  to  the  specimen  and,  with 
suitable  instrumentation,  to  monitor  the  transient  load, 
displacement  and  velocity  of  the  bar  faces  in  contact  with 
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the  specimen.  The  load  is  generated  by  a single  elastic 
pulse  and  the  measurements  rely  only  on  the  accuracy  of 
elastic  wave  theory  in  long  bars.  Since  both  load  and 
deformation  are  measured,  the  stress-strain  relation  is 
derived  directly.  The  average  strain  rate  in  this  test  is 
controlled  by  the  amplitude  to  the  elastic  loading  pulse  and 
the  length  of  the  test  specimen.  By  varying  these  parameters 
it  has  been  possible  to  use  the  technique  to  generate  stress- 
strain  data  in  |:he  range  of  strain  rates  from  approximately 
10  s“^  to  10^  s~  . Since  its  introduction,  the  Hopkinson 
pressure  bar  technique  has  been  extensively  studied, , analyz- 
ed and  modified  such  that  today  it  is  a generally  accepted 
and  utilized  method  in  dynamic  testing.  Modifications  have 
included  its  adaptation  to  tension,  torsion,  and  direct 
shear  testing  as  well  as  compression.  Despite  its  wide  use 
and  the  requirements  for  test  data. in  this  range,  there  is 
as  yet  no  move  toward  standardization  of  this  test  method. 

In  the  intermediate  strain-rate  range  between  convent- 
ional "static"  testing  and  10® s-^ , hydraulic,  ram-type, 
actuated  test  machines  are  now  available.  These  machines  are 
basically  conventional  type  universal  load  frames,  but  with 
special  hydraulic  actuator  and  valving  that  in  combination 
produce  very  high  crosshead  velocity.  Maximum  velocities 
reach  strain  rates  on  the  order  of  10^s”^with  small  test 
specimens.  However,  at  these  crosshead  velocities  consider- 
able attention  must  be  paid  to  the  effects  of  inertia  in  the 
total  system  since  the  accelerations  are  very  high.  The 
frequency  response  of  the  system  and  of  individual  components, 
such  as  the  load  cell  and  extensometer , must  be  accounted  for 
if  the  measured  forces  are  to  provide  information  on  the 
constitutive  properties.  This  is  again  an  area  in  which  some 
move  toward  standardization  of  procedures  is  needed  if  the 
data  are  to  be  widely  used  in  critical  design  applications. 

One  other  type  of  equipment  should  be  mentioned  which 
has  application  in  the  intermediate  strain-rate  range.  This 
is  the  cam-plastomer  (Hockett,  1959) . In  this  system  energy 
is  stored  in  a rotating  mass  which  drives  a cam  having  a 
controlled  profile.  By  means  of  a cam  follower  mechanism 
the  cam  is  coupled  to  the  specimen  and  compresses  it  in  one 
revolution.  An  advantage  of  this  technique  is  that  a 
logarithmic  profile  on  the  cam  and  constant  rotational  speed 
results  in  a constant  true  strain  rate  test.  At  the  higher 
test  speeds,  the  hydraulic  machine  operates  with  an  open 
loop  so  that  the  strain  rate  is  not  positively  controlled. 
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Details  of  these  and  other,  less  widely  used,  test 
techniques  are  given  by  Lindholm  (1971) . Limited  work  has 
been  achieved  under  biaxial  stress  conditions  and  varying 
strain  rate.  This  type  of  data  will  be  necessary  to 
evaluate  generalized  constitutive  equations  needed  for 
structural  analysis. 

Extension  of  mathematical  models  of  classical  plasticity 
to  rate-  or  time-dependent  flow  is  generally  based  upon  the 
concept  of  a yield  function.  In  many  cases  such  a yield 
function  may  be  assumed  of  the  form 


F 


(Jp  1/2 


-1 


c 


(IV- 10) 


where  is  the  second  invariant  of  the  stress  deviator  and 
c is  the  low  rate  or  "static"  yield  stress  in  shear.  The 

plastic  strain  rate,  cij,  is  then  assumed  to  be  proportional 
to  the  magnitude  of  F (where  F is  a measure  of  the  difference 
between  a scalar  measure  of  the  instantaneous  stress  state, 
in  Equation  (IV-10) , and  the  yield  surface, 

=c7  for  positive  values  of  F.  This  stress  difference  ^ 
is  often  called  'overstress ' . In  a general  formulation 
for  one  dimensional  problems  of  viscoplastic  wave  propa- 
gation in  the  x,  direction,  Malvern  (1951)  introduced  the 
functional  dependence 


c'll  = A-(7ii+'r(cTii,  e^i) 


'l'('^ii>  ^ii)  — 0 for  < /(^ii)  , , 

for  <7,,  >/(.,i). 


Malvern,  however,  retained  the  overstress  dependence  in  that 
the  plastic  component  of  the  strain  rate  \j;,  is  finite  only 
when  the  stress  is  greater  than  the  static  stress-strain 
c rve,  = f(^ii).  Perzyna  (1963)  extended  the  Malvern 
formulation  to  generalize  stress  states  as  follows : 


gj'y 

(J')l/2 


<I>(7^)  = 0 for  < c 

¥=  0 for  (J^)i/2  > ^ 


(lV-12) 


In  the  above,  X and  "Y  are  appropriate  elastic  and  viscous 
material  constants  and  e^j  and  cr^j  are  the  deviatoric 
strain  and  stress  tensors,  respectively.  These  continuum 
formulations  for  dynamic  plasticity  rely  on  macroscopic 
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measurements  of  plastic  flow  to  provide  the  arbitrary  con- 
stants and  the  form  of  the  material  functions  (e.g.,  <J>(P)  ) . 

Concurrently,  there  has  been  a significant  effort  to 
describe  macroscopic  flow  in  terms  of  the  controlling 
microscopic  mechanisms,  predominantly  slip  produced  by  the 
motion  of  dislocations,  as  discussed  in  the  preceding 
section.  In  this  approach  we  noted  that  the  plastic  strain 
rate  is  a product  of  the  dislocation  density  and  the  mean 
dislocation  velocity.  By  assuming  a variety  of  functional 
dependencies  of  dislocation  density  and  velocity  on  stress 
and  the  deformation  state,  qualitative  and  often  quantitative 
description  of  macroscopic  behavior  can  be  obtained. 

This  dislocation  dynamics  approach  is  difficult  to 
apply  because  the  relevant  paramenters,  dislocation  density 
and  dislocation  velocity,  are  not  easily  measured  under 
transient  conditions.  Also,  since,  they  occur  as  a product, 
their  independent  contribution  to  the  total  plastic  strain 
rate  is  difficult  to  define.  Nevertheless,  a very  signifi- 
cant increase  in  understanding  of  metal  plasticity  has 
resulted  from  this  approach.  Somewhat  recent  work,  the 
most  notably  by  S.  R.  Bodner  (1968),  combines  the  results 
of  dislocation  dynamics  into  the  framework  of  continuum 
plasticity  theory.  As  a generalization,  Bodner  proposes 
that  the  second  invariant  of  the  plastic  deformation  rate 
deviator  is  a function  of  the  corresponding  stress  invariant, 
or 

(lV-13) 

He  then  develops  incremental  forms  of  the  flow  equations. 

A significant  difference  in  this  approach  and  the  classical 
plasticity  approach  is  the  absence  of  an  explicit  yield 
function.  From  dislocation  dynamics,  plastic  flow  may 
occur  at  any  value  of  finite  stress  and  temperature.  This 
eliminates  the  need  for  partitioning  of  the  stress  space 
into  elastic  and  plastic  domains,  with  the  attendant  book- 
keeping problems  in  associated  computational  procedures. 

Yield  phenomena  can  still  be  adequately  described  by  a 
proper  choice  of  functions  representing  dislocation  multipli- 
cation or  annihilation.  This  approach  also  eliminates  the 
reliance  on  a 'static'  stress-strain  curve  or  yield  surface, 
which  has  no  real  physical  significance  other  than  a ref- 
erence. The  'overstress  hypothesis'  is  not  needed. 
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Another  current  approach  is  to  rely  heavily  on  the 
concepts  of  irreversible  thermodynamics  (e.g.,  Schapery, 

1968  and  Valanis,  1972).  This  approach  brings  appropriate 
attention  to  the  need  to  formulate  the  problem  in  terms  of 
true  state  variables.  These  theories  introduce  internal 
variables  which  characterize  the  deformed  state  of  the 
material  and  enable  the  formulation  of  constitutive  equat- 
ions that  account  for  complex  strain-history  and  temperature 
effects.  For  further  discussions  of  the  thermodynamic 
approach  as  well  as  detailed  reviews  of  other .constitutive 
equations  see  Stricklin  and  Saczalski  (1976)  and  Valanis 
(1976). 

These  approaches  have  all  been  developed  in  the  attempt 
to  describe  the  same  phenomena,  attacking  the  problem  from 
different  conceptual,  bases . At  present,  the  amount  of 
experimental  data  available  concerning  rate'-dependent 
materials  under  combined  stresses  and  time  varying  strains 
and  nonproportional  strain  histories  and  temperature  is  not 
sufficient  to  assess  the  validity  of  any  generalized  theory. 
The  number  of  independent  variables  can  become  very  large 
when  testing  in  multiaxial  stress  states.  Some  guidance 
from  theory  is  needed  in  the  design  of  critical  experiments 
which  can  differentiate  between  alternate  models  of  the 
plastic  deformation  process. 
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4 . Dynamic  Fracture  Analysis  and  Tests 


a.  Introduction 


The  discussion  given  in  this  section  of  the  report 
is  based  on  the  point  of  view  of  fracture  mechanics.  This 
requires  adoption  of  the  premise  that  a flaw  in  the  form  of 
a crack  exists  in  a given  structure  at  the  time  that  a 
specified  high  rate  of  loading  is  applied.  The  distinction 
between  fracture  mechanics  and  a "failure  mechanics"  ap- 
proach is  illustrated  for  an  idealized  situation  in  Figure 
IV-27.  In  failure  mechanics,  a dominant  flaw  does  not  exist 


FIGURE  IV-27  Distinction  Between  Failure  Mechanics 
and  Fracture  Mechanics  Viewpoints 
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and  the  material  is  considered  to  be  completely  homogeneous. 
Then  it  is  often  appropriate  to  develop  a failure  surface 
that  is  a function  of  the  principal  stresses  at  each  point 
in  a body,  not  unlike  the  Tresca  and  von  Mises  criteria  in 
the  theory  of  plasticity.  In  fracture  mechanics,  in  con- 
trast, attention  is  focused  on  the  conditions  at  the  tip  of 
the  dominant  crack- like  flaw.  The  dimensions  of  the  crack 
are  explicitly  included  in  a criterion  for  failure  through 
the  growth  of  this  flaw  which  also  depends  upon  the  "frac- 
ture toughness"  of  the  material;  e.g.,  the  parameter  K in 
Figure  IV-27.  ^ 

All  materials,  of  course,  contain  flaws.  These  may 
arise  from  inherent  defects  in  the  material  or  from  defects 
introduced  by  processing  and  manufacturing  operations.  In 
either  instance,  it  is  likely  that  a more  or  less  uniform 
distribution  of  small  defects  exists  in  every  mater  ial  at 
the  beginning  of  its  service  life.  ’ When  subjected  to  loads 
while  in  service,  these  flaws  may  become  enlarged  and  grow 
in  a stable  fashion  by  creep  or  fatigue  mechanisms.  This 
growth  may  be  aided  by  an  aggressive  environment.  Because 
of  the  presence  of  stress  concentrations  or  high  localized 
loads,  some  flaws  will  grow  more  than  others.  When  at  least 
one  well-defined  crack  exists,  a fracture  mechanics  point  of 
view  becomes  appropriate.  More  definitely,  fracture  mechan- 
ics is  the  discipline  concerned  with  failure  by  the  initia- 
tion and  subsequent  propagation  of  a flaw  in  its  most  dele- 
terious condition:  a sharp  crack. 

Note  that  the  term  fracture  mechanics  is  used  here 
in  a broad  sense.  It  includes  the  growth  initiation,  stable 
growth,  and/or  the  unstable  propagation  of  a crack  in  a body 
without  regard  to  the  constitutive  relation  obeyed  by  the 
material  and  the  environmental  and  loading  conditions  to 
which  it  is  subjected.  This  definition  should  be  distin- 
guished from  linear  elastic  fracture  mechanisms  (LEFM)  which 
is  a special  case,  albeit  the  most  widely  known,  of  the 
general  methodology;  by  definition,  LEFM  is  concerned  with 
the  fracture  of  media  which  are  linearly  elastic  except  for 
a relatively  small  zone  of  damage  or  plastic  deformation 
at  the  crack  tip. 

Just  as  for  any  fracture  problem,  the  characteriza- 
tion of  crack  propagation  initiated  by  high  rates  of  loading 
involves  three  distinct,  but  highly  interrelated  components: 
dynamic  analysis,  dynamic  material  properties,  and  a dynamic 
crack  growth  criterion.  This  report  is  primarily  aimed  at 
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delineating  the  material  properties  necessary  to. cope  with 
high  strain  rate  loadings.  However,  presentation  of  such 
data  by  themselves  is  meaningless.  Equal  prominence  must 
be  given  to  the  theoretical  basis  by  which  the  data  are 
first  extracted  from  a laboratory  experiment  and  subse- 
quently applied  to  evaluate  a structure.  Accordingly,  the 
discussion  in  this  section  contains  an  integrated  descrip- 
tion of  all  three  aspects  of  the  problem. 

While  some  background  material  is  provided,  the  dis- 
cussion primarily  describes  the  extensions  and  generaliza- 
tions of  ordinary  fracture  mechanics  treatments  that. are 
required  in  high  strain  rate  loading  problems.  This  is  in 
accord  with  the  general  approach  taken  in  this  report  which 
is  to  emphasize  the  theoretical  concepts  and  material  prop- 
erties beyond  those  of  the  corresponding  static  problem. 

An  exhaustive  review  of  the  literature  has  hot  been  attempt- 
ed. Instead,  the  intent  was  to  focus  on  a few  key  results, 
emphasizing  reports  and  journal  articles  that  have  appeared 
fairly  recently  in  the  literature,  in  an  attempt  to  provide 
an  up-to-date  perspective  on  the  subject.  The  primary  goal 
is  to  provide  a basis  for  defining  the  directions  that  sub- 
sequent research  on  the  mechanics  of  dynamic  deformation 
and  failure  under  high  strain  rates  must  take. 

b . Background  Discussion,  Definitions  and  Terminology 

of  Fracture  Mechanics 


The  words  "failure"  and  "strength"  are  often  used 
but  seldom  defined  in  a precise  manner.  Consider'  first  the 
term  failure.  There  are  many  different  ways  that  a struc- 
ture can  become  unable  to  adequately  perform  its  primary 
function.  In  each  such  instance,  failure  can  be  considered 
to  have  occurred.  As  noted  in  the  introduction  to  Section 
IV,  the  possible  failure  modes  range  fran  simple  loss  of 
structural  stiffness  due  to  gross  inelastic  deformation; 
e.g.,  yielding  to  the  complete  loss  of  load-carrying  capac- 
ity by  gross  macroscopic  deformation  and  separation;  e.g., 
fracture.  As  an  intermediate  situation,  failure  may  involve 
a reduction  in  load-carrying  capacity  due  to  localized 
damage;  e.g.,  spallation.  Thus,  failure  can  be  gradual  or 
rapid,  may  or  may  not  be  readily  detectable,  and  need  not 
be  catastrophic  in  nature. 
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The  term  strength  is  conventionally  associated  with 
the  load  level  at  which  a failure  occurs  by  some  means  in  a 
standard  test  specimen.  Clearly,  the  strength  will  depend 
upon  the  failure  mode  under  consideration  and  may  be  a 
function  of  many  different  parameters  in  the  test  program, 
such  as  specimen  size.  It  may  also  depend  upon  the  assump- 
tions in  an  analysis  technique  that  is  introduced  to  equate 
the  numbers  actually  measured  in  the  experiment  to  the 
desired  strength  parameter.  Consequently,  even  for  a given 
failure,  the  strength  parameter  may  not  be  a true  material 
property  at  all,  but  will  depend  upon  the  geometry  and  other 
details  of  the  test  specimen  used  to  measure  it.  Hence, , the 
"strength"  of  a material  determined  in  a laboratory  test  is 
often  not  directly  applicable  to  the  same  material  when  it 
is  used  differently  in  service. 

In  view  of  the  foregoing  remarks,  it  is  clear  that 
an  analysis  procedure  that  correctly  extracts  a material 
property  from  a standard  test  procedure  which  is  applicable 
to  engineering  applications  involving  different  geometries 
and  loads  is  absolutely  necessary.  Only  in  this  way  can 
reliable  estimates  of  the  service  failure  loads  be  made 
using  strengths  determined  in  small-scale  laboratory  tests. 
But,  even  though  substantial  experimental  work  has  been 
performed  to  date,  this  capability  does  not  seem  to  exist 
for  materials  subjected  to  high  strain  rate  loading.  It  is 
likely  that  the  least  progress  has  been  made  in  the  most 
critical  problem  area — when  crack-like  defects  exist  in  the 
material.  To  treat  this  subject,  it  is  appropriate  to 
briefly  review  the  subject  of  fracture  mechanics.-  Accord- 
ingly, the  basis  of  existing  fracture  mechanics  techniques 
for  quasi-static  problems  (i.e.,  inertia  due  to  straining 
is  neglected)  will  be  given  next  as  a prelude  to  a more 
thorough  examination  of  the  fracture  mechanics  of  structural 
materials  under  high  rates  of  loading. 

Consider  that  a sharp  planar  crack  exists  in  a linear 
elastic  material  subjected  to  in-plane  symmetric  loading 
remote  from  the  crack.  Consider  further  that  the  crack 
lies  along  the  x axis  and  will  extend  in  a self-similar 
manner.  Then,  in  the  close  proximity  of  the  tip  of  a crack, 
the  component  of  stress  which  acts  in  the  y direction  (i.e., 
the  direction  normal  to  the  crack  plane)  can  be  written  as 
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(r,e) 


I 0 r 0 30  1 

° 2 2 ~ J ■"•••  (IV-14) 

where  0=0  coincides  with  the  prolongation  of  the  crack, 

0 = TT  coincides  with  the  crack  face,  and  r is  the  distance 
from  the  crack  tip.*  Similarly,  the  displacement  normal  to 
the  crack  plane  is,  for  local  plane  strain  conditions. 


v(r,0)  = -i  sin  j^2(l-v)  - cos^  -^j  + ;(IV-15) 

where  la  is  the  shear  modulus  and  v is  Poisson's  ratio.  In 
equations  IV- 14  and  lV-15,  K ■ is  a parametei;  called  the 
stress-intensity  factor.**  li  depends  upon  the  applied  loads, 
the  crack  length,  the  geometry  and,  in  some  cases,  the  elas- 
tic properties  of  the  components.  In  fact,  one  major  element 
of  the  LEFM  approach  is  that  these  parameters  influence  the 
stress  field  close  to  the  crack  tip  only  via  the  stress- 
intensity  factor. 

As , illustrative  examples,  consider  that  a uniform 
tensile  loading  in  the  y direction,  Cy(“)  = CTf, , acts  on  a 
large  body.  Then  for  a single  isolated  crack  'of 

length  2a  on  the  x-axis.  For  an  infinite'  periodic  sequence 
of  cracks  on  the ^x  axis  with  midpoints  a distance  2b  apart, 

= Gn  tan  ' (■^,)  . And,  for  ,a  central  crack  in  a long 

T y ’ 2,lo  — ■ jl/  TT^ 

strip  of  width  2w,  =“  a(,Vrra  sec  {— ) . 

The  second  major  element  of  LEFM  is  the  connection  be- 
tween the  local  crack  tip  quantity — the  parameter  — and  a 
global  energy  balance  related  parameter  G.  This  arises  as 
follows.  Consider. a virtual  extension  Aa  of  the  crack. 

Then,  under  conditions  where  LEFM  applies,  G is  defined  as 


* Several  standard  formulas  are  given  in  this  background 
discussion  withoug  citing  publications;  they  can  be  found 
in  recent  books  on  fracture  (Hertzberg,  1976  and  Rolfe 
and  Barson,  1977) . 

**  Conventionally,  the  subscript  I is  used  in  order  to 
distinguish  between  this  mode  of  loading,  the  so-called 
opening  mode,  and  Kjj  and  Kjjj  which  correspond  to  in-plane 
and  out-of-plane  shear  loadings  (i.e.,  the  skew- symmetric 
and  antiplane  shear  modes,  respectively). 
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G 


1 


(IV- 16) 


dW  dU 
b da  da 

where  W denotes  the  work  done  by  the  applied  tractions,  U 
denotes  the  elastic  strain  energy  in  the  body,  and  b is 
the  length  of  the  advancing  crack  front;  the  quantity  G is 
observed  to  be  the  excess  energy  dW  - dU  divided  by  the 
area  of  crack  extension,  dba.  Therefore,  assuming  inelastic 
behavior  is  limited  essentially  to  the  crack  plane,  and 
neglecting  heat  flow  and  kinetic  energy,  G must  be  the  work 
done  by  the  traction  cr  (r,o)  on  separating  the  material  at 
the  crack  tip  (per  unit  of  new  surface) . For  a linear  elas- 
tic body  this  work  is 


(r,o)  v(Aa-r,TT)dr  (lV-17) 


^ a 

G = 2 J ~ a 

Aa-*o  Aa  ■'  2 y 

o 


Upon  substituting  from  Equation  IV- 14  and  IV-15,  performing 
the  integration^  and  taking  the  limit,  it  is  found  that 


G 


(lV-18) 


Now,  it  may  be  taken  as  a basic  principle  of  LEFM  that  the 
crack  will  extend  when,  and  only  when,  G exceeds  a critical 
value  of  energy  G . Then,  in  view  of  Equation  lV-18,  it  is 
exactly  equivalent  to  consider  that  crack  extension  takes 
place  when  exceeds  a critical  value.  Denoting  this  latter 
parameter  by  K , we  can  write 


^ 

1-v^ 


G 

c 


(IV- 19) 
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The  energy  is  a material  property  and  thus  Kj^  is  also  a 
material  property.  The  latter  quantity  is  called  the’ plane- 
strain  fracture  toughnesst 

In  essence,  LEFM  provides  a criterion  for  fracture 
by  crack  extension  by  equating  a function  of  the  applied 
loads  and  the  geometry  of  a structural  component — the 
stress-intensity  factor — to  a material  property — the  frac- 
ture toughness.  To  illustrate  the  application  of  this  idea, 
the  stress- intensity  factors  quoted  above  can  be  used  to 
express  the  fracture  stress  ^f  for  materials  with  given 
toughnesses  and  geometries.  These  are  obtained  by  simply 
replacing  Kj  by  and  taking  the  applied  stress  equal  to 

CTf.  The  results  are 

(IV-20a) 
(IV-20b) 


(IV-20C) 


^f 


= Kjj,  (na) 


-1/2 


isolated  crack 


Kjc  1 2b  tan 


TT  a 
2b 


V-1/2 


, cqlinear  cracks 


= K 


Ic 


finite-  width  strip 


In  an  ideal,  completely  elastic  body  for  which  the 
material  separates  along  a single  plane,  one  writes  G^-2y 
where  y is  now  the  so-called  surface  energy.  Hence, 
from  Equations  IV-19  and  IV-20a, 


(Tf  = 


2 Ey 
5:a  (l-V^) 


(lV-21) 


* According  to  ASTM  E399,  plane  strain  conditions  are 
achieved  near  the  tip  of  a through-the-thickness  crack  when 
the  plate  thickness  B s:  2 . 5 (Kjc/ay)^,  where  is  the  yield 

stress.  A similar  expression  has  been  used  for  fast  running 
cracks  where  is  interpreted  as  the  dynamic  yield  stress 

corresponding  to  a strain  rate  within  or  adjacent  to  the 
plastic  zone. 
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which  is  the  familiar  Griffith  equation  except  for  the 
factor  l-u^ ; this  factor  accounts  for  the  fact  that  the 
crack  tip  is  in  a state  of  plane  strain,  rather  than  plane 
stress.  In  metals,  however,  it  is  generally  true  that 
G(-;»2y  because  of  the  complex  failure  processes  that  occur 
at  the  crack  tip  when  crack  extension  initiates.  In  situa- 
tions where  the  loads  are  not  applied  symmetrically  with 
respect  to  the  crack  plane,  more  than  one  type  of  stress 
intensity  factor  prevails;  indeed,  in  some  cases  all  three 
factors,  Kj,Kj];,  and  Kjjj  may  be  needed.  The  use  of  the 
strain  energy  density  factor  S associated  with  an  element 
ahead  of  the  crack  has  been  proposed  as  a criterion  of 
fracture  for'  static  and  dynamic  problems  (Sih,  1977) . Crab 
instability  is  taken  to  coincide  with  the  failure  of  this 
element  when  the  stored  energy  reaches  a critical  value  of 
S,  denoted  as  S^,,  which  is  assumed  to  be  a characteristic 
property  of  the  material.  For  the  case  of  an  angle  crack 
of  length  2a  subjected  to  uniaxial  dynamic  loading,  the 
fracture  stress  may  be  expressed  as 


(lV-22) 


c 


a 


f 


TiaF  (/3,  t,  e^) 


in  which  F is  an  analytically  determined  function  that  de- 
pends on  time  as  well  as  the  geometrical  parameters  (cf. 
Section  IV.B.4.e).  In  Equation  IV-22,  B is  the  angle 
between  the  crack  plane  and  the  axis  of  loading,  while  0 
is  the  direction  of  crack  initiation  determined  from  the 
condition  of  minimum  strain  energy  density.  As  discussed 
in  Section  IV.B.4.e  other  criteria  have  also  been  proposed 
for  mixed  mode  problems . 


One  of  the  basic  considerations  in  LEFM  is  that 


the  amount  of  plastic  deformation  surrounding  the  crack  tip 
must  be  small  enough  that  the  perturbation  of  the  linear 
elastic  stress  and  displacement  fields  can  be  neglected 
when  calculating  energy  release  rate  G.  While  this  is  es- 
sentially true  for  brittle  materials,  it  will  not  be  a good 
assumption  for  ductile  materials.  In  the  latter  case,  how- 
ever, a plasticity  correction  to  crack  Idngth  known  as  the 
Dugdale  or  strip  yielding  model  has  been  usefully  employed. 
In  this  model,  plastic  deformation  without  strain  hardening 
is. considered  to  occur  in  a thin  strip  which  is,  in  effect. 
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an  extension  of  the  crack.  Then,  the  plastic  zone  length  I 
for  an  isolated,  through- the-thickness  crack  in  a body  under 
uniform  tension  is  given  by 


(IV-23a) 


where  cTp  is  the  applied,  tensile  stress  and  ay  is  the  yield 
stress  of  the  material.  If  l is  large  compared  to  sheet 
thickness,  plane  stress  theory  is  employed  to  compute  the 
crack  opening  displacement  at  the  tip. 


S 


6t 


8 

TT  E 


log 


sec 


(IV-23b) 


A critical  value  of  the  crack-opening  displace- 
ment fit  is  sometimes  offered  as  an  alternative  fracture  cri- 
terion. .But  when  i/a  « 1,  it  is  found  from  Equation  IV-23b 
and  IV-20a  that  f = K^/Eg  , which  shows  that  a critical 
crack-opening  displacement  criterion  is  equivalent  to  a 
criterion  on  the  stress  intensity  factor  when  the  amount  of 
plastic  deformation  is  small.  Hence,  there  is  no  basic  con- 
tradiction between  these  two  approaches. 

We  should,  add  that  a strip-like  yield  zone  is  cer- 
tainly not  universally  observed  and,  in  fact,  does  not  even 
represent  the  most  usual  case  in  metals;  it  appears  to  be 
limited,  to  thin  sheets  of  moderate  or  low  strain-hardening 
material.  When  this  is  not  the  case,  the  plastic  zone  has 
;a  two  or  three-dimensional  character,  and.  a very  involved 
,ftiathematical  treatment  may  be  required. 

An  even  more  useful  and  general  approach  to  char- 
acterizing the  critical  condition  for  fracture  initiation  in 
the  presence  of  large-scale  plasticity  is  based,  on  the  so- 
called.  J-integral  (Rice,  1968;  see  also  Rolfe  and.  Barsom, 
(1977)).  Initiation  is  assumed  to  occur  when  J reaches  a 
critical  value  for  small-scale  yielding,  and 

therefore  this  criterion  is  consistent  with  LEFM.  It  should, 
be  noted  that  the  use  of  J for  the  opening  mode  is  limited 
to  plane  strain  or  generalized,  plane  stress  problems,  and 
to  situations  where  the  deformation  theory  of  plasticity  is 
a realistic  approximation  of  the  elastic-plastic  behavior. 

It  is  of  interest  here  to  evaluate  the  J-integral  for  the 
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case  of  a body  in  which  material  separation  occurs  entirely 
within  a thin  layer  or  strip  ahead  of  the  crack  tip,  and 
outside  of  this  layer  the  material  is  either  elastic,  or  the 
deformation  theory  applies.  The  J— integral  becomes  (Rice, 
1968) 


o-  ^5  d^^ 

y d ix 


(lV-24a) 


where  Cy  is  the  reaction  of  the  separating’’ material  on  the 
intact  (elastic  or  plastic)  material  and  j is  the  separation 
between  the  edges  of  the  intact  material;  6^13  the  value  of 
6 at  the  tip.  If  Cy  = cry(6),  then  Equation  IV-24a  becomes 

J =1^^  cr  (6)d6  (lV-24b) 

_ t)  ^ 

which  is  the  work  done  in  pulling  to  failure  a ligament  of 
material  in  the  thin  layer  at  the  crack  tip;  note  that  the 
j-value  does  not  include  the  W3rk  of  plastic  deformation  of 
material  outside  of  the  thin  strip,  if  we  further  assume 
0'y=  constant  s (O'Y)then  J = CTySt,  which  implies  a 
criterion  is  equivalent  to  a crack  opening  displacement 
criterion.  This  fact  that  J equals  the  work  done  in  separ- 
ating the  material  at  the  tip,  and  not  the  work  of  large- 
scale  plastic  deformation,  may  be  the  reason  for  the  success 
of  the  criterion, 

A substantial  portion  of  current  fracture  testing 
work  is  focused  on  the  recommendations  of  the  Pressure  Ves- 
sel Research  Committee  (PVRC)  and  the  requirements  of 
Sections  III  and  XI  of  the  ASME  Boiler  and  Pressure  Vessel 
Code,  These  are  based  largely  on  the  concepts  and  method- 
ology of  LEFM.  More  specifically,  the  PVRC  recommendations 
and  the  ASME  Code  criteria  represent  the  fracture  toughness 
of  pertinent  steels  by  Kjj^,  the  so-called  "reference  value" 
of  the  critical  stress-intensity  factor  for  unstable  crack 
propagation.  The  curve  of  Kj-^  as  a function  of  the  temper- 
ature relative  to  the  drop-weight  nil  ductility  temperature 
(NDT)  test  (see  Section  IV, B, 4 , C, (2 ) (6 ) as  derived  in  the 
PVRC  recommendations  is  shown  in  Figure  lV-28,  As  can  be 
seen  from  Figure  lv-28,  Kjj^  is  simply  a lower  bound  to  the 
^v^ilskle  fracture  toughness  values  for  nuclear  pressure 
vessel  steels.  The  most  important  contributions  to  this  data 
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SOURCE:  Rolfe  and  Barsom,  1977 


base  are  the  parameters  Kj^  and  Kja/  where  Kjd  is  the  value 
of  the  stress  intensity  factor  for  crack  growth  initiation 
under  dynamic  loading,  and  Kja»  is  an  approximation  to  the 
value  at  the  arrest  of  a rapidly  propagating  crack. 

In  view  of  the  practical  importance  of  the  results 
given  in  Figure  IV-28,  the  basis  for  determining  and 

Kja  values  is  quite  important.  Both  values  are  purported 
to  be  material  properties.  As  noted  above,  Kj^  is  the 
toughness  value  measured  by  rapid  (dynamic)  loading  of  a 
specimen  with  a stationary  (nonmoving)  precrack.  K is 
the  toughness  value  measured  on  a specimen  immediateL^y  after 
the  arrest  of  a propagating  crack.  in  both  cases,  static 
analyses  have  been  used  to  obtain  the  values  of  stress  in- 
tensity factors.  In  view  of  the  role  of  these  values  in 
defining  the  Kjj^  curve,  their  origin  in  a static  analysis 
when,  in  fact,  the  experiments  in  which  they  are  determined 
are  highly  dynamic  is  a point  of  some  concern.  In  contrast, 
the  parameter  Kj^j,  the  dynamic  toughness  of  a rapidly  pro- 
pagating crack,  is  entirely  based  on  dynamic  (i.e.,  inertia 
forces  included)  analysis  procedures. 

c.  Crack  Growth  Initiation  Under  High  Loading  Rates 

(1)  Dynamic  Fracture  Concepts.  At  present,  both 
crack  initiation  and  crack  arrest  studies  are  analyzed 
predominately  in  a manner  mathematically  identical  to  a 
crack  initiation  analysis  using  LEFM.  Thus,  it  is  supposed 
that  crack  arrest  is  the  reverse,  in  time,  of  crack  growth 
initiation.  The  arrested  crack  is  treated  as  a stationary 
crack  whose  stress-intensity  factor  is  derived  from  a static 
analysis.  The  crack  arrest  toughness  Kjg  is  interpreted  in 
a way  that  is  analogous  to  Kj^-,.  Similarly,  in  analyzing 
crack  growth  initiation  under  high  loading  rates,  crack-tip 
inertia  effects  are  neglected;  this  leads  to  the  high  load- 
ing rate  crack  growth  initiation  toughness  Kj^j,  which  is 
analogous  to  Kjc- 

The  use  of  Kja  as  a material  property  is  based 
on  the  assumption  that  it  is  uniquely  related  to  the  minimum 
value  of  fracture  toughness  as  a function  of  crack  speed. 
This  assumption  is  supported  in  the  Kj^  measurements  that 
have  been  reported  by  Crosley  and  Kipling  (1975).  But, 
because  inertia  effects  are  neglected,  Hahn,  et  al.  (1975- 
76)  have  argued  that  Kj^  cannot  be  a true  material  property 
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and  will  instead  be  a function  of  the  stress-intensity 
factor  at  the  onset  of  rapid  crack  propagation  and  of  the 
test  piece  geometry. 

A more  correct  viewpoint  for  crack  arrest  is  based 
on  an  extension  of  LEFM  to  the  dynamic  case.  This  delineates 
the  role  of  a fundamental  dynamic  toughness  parameter 
there  is  evidence  that  Kj(j  = where  a denotes  the 

instantageous  crack  propagation  speed. 

The  emphasis  in  this  section  is  on  the  measurements 
and  analyses  associated  with  the  parameter.  Such  meas- 

urements are  made  with  impact  tests  of  various  kinds.  It  is 
interesting  that  impact  tests  traditionally  have  been  used 
■ only  in  a qualitative  fashion  for  the  acceptance  of  materials 
to  preclude  fracture  under  various  design  codes.  But,  be- 
cause such  design  procedures  (based  on  tensile  strength  with 
large  factors  of  safety)  do  not  always  suffice  when  flaws 
are  present  and,  in  addition,  do  not  allow  any  flexibility 
to  scale  up  the  fracture  resistance,  they  have  begun  to 
give  way  to  quantitative  procedures  based  on  LEFM.  In  sup- 
port of  a fracture  mechanics  approach,  more  precise  infor- 
mation has  been  required  from  impact  testing. 

The  applicability  of  fractube  resistance  data  gen- 
erated by  the  high  strain  rates  in  impact  tests  may  be 
questioned,  for  statically  loaded  structures.  Stress-inten- 
sity rates  from  3 x lO'^to  1 x 10®  MN*m  ®/®/s(3  x 10“^  to 
1 X 10®ksiVin./s)  , causing  fracture  to  occur  in  1 to  0.1 
ms,  are  common  in  specimens  under  conventional  impact  load- 
ing rates.  These  relatively  high  rates  may  be  seyeral  orders 
of  magnitude  higher  than  those  encountered  in  engineering 
structures.  However,  catastrophic  fractures  of  large 
welded  structures  are  related  to  the  dynamic,  high-strain- 
rate  properties  of  the  materials.  This  is  a consequence  of 
the  steep  strain  gradient  that  exists  ahead  of  a sharp  crack 
which  subjects  the  material  ahead,  of  a running  crack  to 
quick  high  strain  rates  even  at  relatively  moderate  applied, 
loading  rates . 

(2)  Dynamic  Fracture  Initiation  Testing  Methods.  In 
this  subsection,  the  more  prominent  of  the  testing  techniques 
used,  to  extract  dynamic  crack  growth  initiation  data  on 
engineering  materials  will  be  briefly  described.  Attention 
will  be  primarily  addressed,  to  Kj^  determination  methods; 
but,  because  of  the  correlations  that  exist  between  the 
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several  different  parameters  determined  in  impact  tests, 
attention  will  not  be  limited  exclusively  to  these.  Thus, 
in  the  following,  the  Charpy  test,  the  drop  weight  tea^r 
test,  and  the  dynamic  tear  test  will  be  discussed.  The 
brief  treatment  given  here  is  based  in  large  part  upon  the 
extensive  exposition  given  by  Lange  (1976) . (See  also 
Hertzberg,  1976.) 

(a)  Charpy  V-Notch  Test.  The  Charpy  test  is 
the  oldest  standardized  impact  test  in  general-  use  (ASTME23- 
33T) . The  energy  criterion  provided  the  first  correlation 
between  the  results  of  a laboratory  fracture  test  and  the 
service  performance  of  the  material  in  a structure;  this 
energy  was  correlated  with  the  initiation,  propagation,  and 
arrest  of  fracture  in  the  World  War  II  Liberty  ships  and 
T-2  tankers.  Fractures  initiated  from  small  flaws  when  the 
Cv  energy  was  less  than  14  J (10  ft-lb) , and  arrested 
cracks  were  found  only  when  the  Cv  energy  was  more  than  27  J 
(20  ft-lb) . These  correlations  have  been  used  to  justify 
the  use  of  a 20  J (15  ft-lb)  criterion  to  preclude  brittle 
fracture  in  all  types  of  welded  steel  structures.  Unfor- 
tunately, while  the  20  J (15  ft-lb)  criterion  was  found  to 
be  valid  for  ship  steels',  it  is  not  generally  useful. 

Although  there  are  several  specimen  designs,  the 
most  commonly  used  specimen  today  is  shown  in  Figure  IV-29; 
it  is  a 10-mm-square  specimen  with  a notch  2-mm  deep  having 
a 0.25-mm  root  radius.  The  relatively  dull  notch  in  the 
standard  Cv  specimen  is  one  reason  the  transition . of  a C^ 
energy-temperature  relation  is  not  indicative  of  structural 
performance.  To  improve  this  feature  of  the  C^  specimen, 
much  recent  research  has  been  conducted  with  fatigue-cracked 
C^  specimens.  Fatigue  cracking  the  C^  specimen  tends  to 
steepen  the  transition  region  and  shift  it  to  a higher 
temperature.  But,  the  plane-strain  measuring  capability  of 
a 10-mm  section  remains  at  2.0'\/iran  (0.4\/In.)  = Kj^j/ffyd- 
This  limitation  in  plane-strain  measurement  can  be  an  impor- 
tant restriction  to  the  characterization  of  the  transition 
region  of  certain  steels. 
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FIGURE  IV-29  Geometry  of  Precracked  Charpy 
Specimen 

SOURCE:  Ayres,  1976 

The  instrumented  precracked  Charpy  test  employs  a 
Charpy  V-notch  specimen  with  a fatigue  crack  grown  to  ap- 
proximately the  center  of  the  specimen.  The  tup  or  striker 
is  instrumented  so  that  the  force  on  the  tup  is  recorded 
as  a function  of  time.  A typical  load  versus  time  curve  for 
a test  in  the  brittle-ductile  transition  region  is  shown  in 
Figure  IV-30.  The  curve  is  used  to  determine  the  fracture 
toughness  Kjd  via  the  static  elastic  three-point  bend  beam 
relation  for  the  static  stress-intensity  factor  given  by 
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(b)  Drop  Weight-Nil  Ductility  Transition  Temper- 
ature (NDT)  test.  The  Drop  Weight-Nil  Duc- 
tility Transition  (NDT)  Temperature  Test  was  devised  at  the 
Naval  Research  Laboratory  in  the  early  1950's  (Lange,  1976). 
This  test,  currently  covered  by  ASTM  E208,  was  the  first 
test  to  be  called  a Drop  Weight  test.  The  test  can  be  used 
to  determine  the  temperature  at  which  a steel  sample  can 
plastically  deform  under  a dynamic  load  in  the  presence  of 
a small  crack,  or  it  can  be  used  as  a pop-in  test. 

The  pop-in  crack  of  the  DWT-NDT  specimen  is 
formed  when  a brittle  weld  bead  fractures  as  the  specimen  is 
loaded  in  a drop-weight  machine.  The  crack  from  the  weld 
bead  is  either  arrested,  or  it  continues  to  propagate  across 
to  the  top  surface  of  the  specimen.  If  the  crack  does  not 
propagate  to  one  of  the  top  corners  of  the  s.pecimen  before 
the  surface  strain  approaches  2 percent,  the  deformation  is 
stopped  by  an  arrestor  block  on  the  anvil  and  the  test 
result  is  called  a "no  break".  Specimens  are  subsequently 
tested  at  lower  temperatures  until  a "break"  occurs.  This 
defines  the  NDT.  The  DWT-NDT  test  is,  therefore,  a "go — 
no-go"  appearance  test. 

Another  important  function  of  the  DWT-NDT  test  is 
to  determine  the  plane-strain  fracture  limit  for  sections 
16  mm  (5/8-inch)  thick.  Because  an  NDT  performance  trans- 
lates to  brittle  fracture  initiating  from  a small  flaw,  this 
level  of  fracture  resistance  is  = 2.5  /mm  (0 . 5Vinch)  . 

Service  experience  has  shown  that  when  the  steel  in  a welded  ■ 
structure  has  a fracture  resistance  less  than  this  critical 
value,  catastrophic  fractures  can  be  expected. 

(c)  Drop  Weight  Tear  Test.  The  Drop-Weight  Tear 
Test  (DWTT) , which  employs  an  oversized  Charpy  specimen,  is 
used,  to  define  the  temperature  transition  region  for  ferritic 
steels  in  sections  from  3.18  to  19.1  mm  (0.125  to  0.75  inch). 

(The  samples  in  the  DWTT,  and  those  in  the  DT  test  discussed 
below,  are  broken  using  pendulum  or  dropweight  machines.) 

Its  primary  application  has  been  to  determine  the  temperature 
at  which  shear-type  fractures  occur  in  steel  for  linepipe. 

For  many  pipes,  the  propagation  of  cracks  that  can  be  initi- 
ated from  a number  of  accidental  sources  will  run  for  ex- 
tended distances  unless  the  appearance  of  the  fracture  is 
more  than  80  percent  shear.  Although  this  criterion  indi- 
cates the  performance  of  some  line  pipes,  the  test  and  its 
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fracture  appearance  criterion  do  have  general  application 
for  establishing  the  transition  region  for  many  steels  in 
thin  sections.  A correlation  between  Charpy  energies  and 
DWTT  energies  for  these  materials,  obtained  by  Wilkowski 
and  reported  by  Popelar,  et  al.  (1977) , is  shown  in 
Figure  IV-31, 

When  a conventional  steel  having  a yield  strength 
under  827  MPa  (120  ksi)  is  used  in  sections  less  than  16  mm 
(5/8-inch)  thick,  the  transition  in  fracture  resistance 
tends  to  be  restricted  to  a narrower  range  of • temperatures 
than  when  the  same  steel  is  rolled  to  a thicker  gage.  There- 
fore, for  many  structural  applications,  a fracture  appear- 
ance criterion  is  all  that  is  needed  to  predict  a brittle 
or  ductile  performance. 

Fracture  resistance  rises  so  sharply  in  thin 
sections  that  a 50  percent  shear  temperature  criterion  is 
a very  readily  determined  and  unambiguous  parameter  for 
predicting  service  performance  when  a full-section  specimen 
is  used. 

The  DWTT  is  not  a test  for  generalized  use.  It 
cannot  be  used  to  establish  fracture  resistance  quantita- 
tively at  upper  shelf  temperatures  or  to  elevate  transitions 
for  steels  that  do  not  undergo  a sharp  transition  in  frac- 
ture appearance.  Energy  measurements  from  the  DWTT  are  not 
meaningful  when  the  small  notch  causes  the  ligament  to  de- 
form plastically  prior  to  fracture  even  when  the  fracture  is 
brittle.  Although  the  test  has  its  limitations,  it  is  use- 
ful for  determining  the  temperature  transition  region  for 
plain  carbon  and  low-alloy  steels  in  sections  thinner  than 
19  mm  (3/4- inch). 

(d)  Dynamic  Tear  Test.  The  Dynamic  Tear  Test 
(DT)  also  uses  an  oversized  Charpy  specimen;  the  notch  is 
deeper  than  in-  the  DWTT  and  is  embedded  within  a titanium 
embrittled  electron  beam  weld.  It  was  develooed  to  fill 
the  need  for  a practical  test  that  could  precisely ' measure 
fracture  resistance  over  a broad  range,  including  that  en- 
countered in  the  transition  region  of  the  ferritic  steels. 

For  practical  reasons,  the  specimen  is  not  instrumented,  and 
the  total  energy  used  to  fracture  a specimen  is  the  criterion 
of  fracture  resistance.  This  empirical  value  of  fracture 
resistance,  DT  energy,  is  translated  into  structural  pa- 
rameters by  correlation  with  an  analytical  parameter  or 
with  specific  structural  performance. 
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figure  IV-31  Correlation  Between  Absorbed  Energy  in  the  Drop-Weight  Tea 
Test  and  the  Charpy  Upper  Shelf  Energy  for  Pipeline  Steels 
SOURCE:  From  Wilkowski  as  reported  by  Popelar  et  al,  1977 


The  primary  intent  in  developing  the  DT  test  was 
not  to  develop  an  inexpensive  or  test,  but  to  pro- 

vide a sensitive  and  reliable  fracture-resistance  criterion 
for  the  elastic-plastic  and  plastic  regimes.  Although  the 
DT  energy  criterion  is  empirical,  it  has  correlated  very 
well  with  values  for  various  steels,  titanium,  and 

aluminum  alloys.  The  correlations  cover  a sufficiently 
broad  range  of  alloys  in  each  base  metal  to  justify  the 
development  of  generalized  diagrams  for  structural  analysis. 
The  diagrams  are  called  Ratio  Analysis  Diagrams  where,  by 
a simple  graphical  procedure,  a DT  energy  value  can  be  trans- 
lated into  a Kj^/Cy  value  or  a critical  flaw  size. 

When  practical,  the  DT  test  uses  a full-thickness 
specimen  so  that  the  crack  tip  constraint  in  the  specimen  is 
the  same  as  that  in  the  structural  element  of  interest. 
Obviously,  this  is  not  very  practical  when  thick  sections 
are  involved.  So,  extensive  studies  have  been  conducted  to 
establish  procedures  to  compensate  for  the  effect  of  speci- 
men size.  These  studies  have  included  tests  on  12-inch- 
thick  sections,  and  have  led  to  methods  for  extrapolating 
DT  energy  values  obtained  with  a subsize  specimen  to  full- 
section  structural  performance. 

A correlation  between  and  DT  energy  for  the 

conventional  structural  steels  remains  in  the  developmental 
stage.  These  materials  are  strain-rate  sensitive,  and  this 
characteristic  has  always  complicated  attempts  to  measure 
their  fracture  resistance  properties  in  terms  of  Kjd*  If 
conditions  are  such  that  fracture  initiates  by  clevage, 
then  the  fracture  resistance  is  low;  if  the  fracture  init- 
iates by  microvoid  coalescence,  then  the  fracture  resistance 
is  high.  Therefore,  the  fracture  mode  must  be  consistent 
for  the  two  test  methods  or  correlation  cannot  be  developed. 

In  addition,  the  strain  rates  associated  with  impact  testing 
are  sufficient  to  cause  complex  bending  waves  in  the  specimen, 
which  makes  the  analysis  of  tests  quite  difficult.  Until 

more  direct  determinations  are  made,  correlations  between 

Kjd  snd  DT  energy  will  be  of  limited  significance. 

A correlation  between  and  DT  energy  was  devel- 

oped with  high-strength  steels  that  were  not  sensitive  to 
strain  rate.  The  initial  relationship,  which  was  developed 
with  wrought  products,  was  later  substantiated  by  data  using 
high-strength  cast  steels.  The  data  from  the  study  on  cast 
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steels  are  shown  plotted  on  the  original  relationship  devel- 
oped with  wrought  steels  in  Figure  IV-32.  The  correlation 
extends  to  a value  of  110  MPa  -^/in  (100  ksi  Vinch)  , where 

the  16-mm  {5/8-inch)  DT  specimen  becomes  substantially 
plastic  and  the  DT  energy  increases  at  a more  rapid  rate 
than  does,  the  plane-strain  parameter  for  fracture^ toughness . 

(e)  Precracked  Tension  Bar.  A test  method  has 
been  developed  by  Costin  et  al.  (1977)  for  accurately  estab- 
lishing the  dynamic  fracture  initiation  properties  of  struc- 
tural metals  at  extremely  high  loading  rates.  (Although  this 
method  is  new,  it  is  included  here  because  of  its  potential 
usefulness.).  The  apparatus,  which  is  shown  in  Figure  IV-33, 
is  an  adaptation  of  the  Koisky  pressure  bar  ’ (split-Hopkinson 
bar)  in  which  a one-inch  diameter  round  bar  specimen  with  a 
pre-fatigued  circumferential  notch  is  loaded  to  failure  by 
the  rapidly  rising  tensile  pulse  resulting  from  an  explosive 
detonation.  Using  the  standard  Koisky  technique,  the  average 
stress  at  the  fracture  site  is  measured  as  a function  of 
time.  Crack  opening  displacement  is  measured  by  optical 
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16  mm  DYNAMIC  TEAR  ENERGY  (J) 


FIGURE  IV-32  Correlation  Between  Kjc  and  5/8 

. DT  Energy  for  High-Strength  Wrought 
and  Cast  Steels 


SOURCE:  Lange,  1976 
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FIGURE  IV-33  Schematic  Diagram  of  Apparatus 
SOURCE:  Costin  et  al,  1977 
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means,  as  a function  of  time,  thus  yielding  a complete  load- 
displacement  record  for  each  test.  From  the  data  the  crit- 
ical value  of  the  crack  tip  stress  intensity  factor,  Kj^^, 
at  loading  rates,  Kj,  in  excess  of  lO^MPav^T/s  (10®ksi  vln/s) 
may  be  obtained.  This  is  nearly  two  orders  of  magnitude 
faster  than  has  been  achieved  by  other  standard  techniques. 

Analysis  of  the  apparatus  indicates  that  when  the 
pulse  length  is  long  compared  to  the  notch  width,  the  loading 
of  the  notched  section  may  be  viewed  as  quasi-static.  Thus, 
determination  of  fracture  parameters  may  be  accomplished  by 
using  the  local  load— displacement  record  in  conjunction  with 
existing  formulas.  Costin  ef  al.  (1977)  employ  the  standard 
static  relation 

P ' I 

= —5—  '/ttR  F (IV-27) 

I ttk 

when  the  pla^jiic  zone  is  very  small;  in  Equation  IV-27, 

P = load,  R = radius  of  remaining  circular  ligament, 

F = function  of  R/D,  where  D is  the  outside  diameter  of  the 
bar.  They  use  the  J- integral  when  the  yield  zone  is  not 
small  cottipared  to  R. 

(3)  Effect  of  Specimen  Inertia  on 'Dynamic 
Fracture  Initiation.  Dynamic  analysis  of  the  initiation  of 
crack  growth  under  impact  loadings  generally  requires  the 
use  of  numerical  methods.  The  finite  element  method,  dis- 
cussed in  Section  B-.6,  is  currently  being  used  by  several 
investigators.  The  power  and  flexibility  of  the  finite- 
element  method  can  be  economically  and  effectively  increased 
by  the  development  of  special  crack-tip  elements  with  an 
embedded  characteristic  crack-tip  stress  singularity.  Such 
elements  are  now  popularly  called  cracked  finite  elements. 

Results  obtained  with  a Itiodel  of  an  instrumented 
precracked  Charpy  specimen  by  Anderson,  et  al.  (1975)  are 
shown  in  Figure  IV-34.  Most  investigators  have  assumed 
failure  of  the  specimen  occurs  when  the  peak  load  develops. 
The  peak  load  is  then  used  in  a static  analysis  of  a three- 
point  bending  specimen  to  estimate  the  value  of  K^.  at  which 
the  crack  begins  to  grow;  as  an  example,  this  value  is  shown 
as  Kj  (static)  in  Figure  IV-34.  In  some  cases  the  static 
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Charpy  Specimen 


value  turns  out  to  be  twice  the  true  critical  stress-inten- 
sity factor  Kjj,  for  static  loading,  and  this  leads  to  the 
erroneous  conclusion  that  the  material  has  a dynamic  tough- 
ness about  twice  as  great  as  the  static  toughness.  Figure 
IV-34  shows  that  this  is  very  likely  not  the  case.  At  the 
time  that  the  load  peaks  (less  3|i  seconds  for  signal  travel 
time) . K is  about  180  x 10®  dynes/mm® which  is  not  twice 
the  static  toughness  for  mild,  steel,  but  it  is  nominally 
about  equal  to  it.  This  rather  dramatically  points  out  the 
significance  of  inertia  effects  in  fracture  specimens  when 
the  loads  are  applied  suddenly.  The  initial  dive  to  nega- 
tive values  of  Kj  just  after  impact,  which  is  exhibited,  in 
Figure  IV-34,  is  an  indication  of  just  how  sensitive  the 
dynamic  cracked,  element  is  to  wave  phenomena.  Until  the 
compressive  wave  passes  through  the  thickness  and.  is  re- 
flected off  the  back  side,  the  tendency  of  the  crack  will 
be  to  close. 

(4)  Material  Property  Data  for  Dynamic  Crack 
Growth  Initiation.  Krafft,  et  al.  (1970)  calculated,  the 
^id  values  shown  in  Figure  IV-35  from  unnotched,  tensile 
properties  measured,  at  high  strain  rates  using  a theoretical 
relation  between  and  yield,  stress.  Since  then,  direct 

measurements  of  have  been  reported,  for  a number  of  steels 

for  loading  rates,  ]^,  up  to  about  10® to  10®MPav(n^s.  These 
rates  are  the  upper  limits  of  the  experiments  and.  the  anal- 
yses that  have  been  used  in  the  Krafft,  et  al.  studies.  It 
is  of  interest  to  note  that  a different  test  method  has  been 
suggested  by  Shockey  and  Curran  (1973)  which  might  be  able 
to  provide  data  at  of  10®  to  10®  MPavr^s. 

^Id.  ^^ta  for  several  low  to  medium  strength  steels 
have  been  obtained  by  Shoemaker  and.  Rolfe  (1971)  . K,.  values 

3-n  A533B  steel  (HSST  Plate  02)  have  been  measured  by 
Shabbits  (1970)  at  a series  of  temperatures  and  loading 
rates.  Crosley  and.  Ripling  (1975)  also  obtained  values 

on  this  same  material  during  the  course  of  their  crack  ar- 
rest testing.  In  addition,  Kj^  values  are  also  available 
for  this  material  from  the  work  of  Shabbits  (1970)  and 
from  J-integral  and  equivalent  energy  tests.  All  of  these 
data  for  this  heat  of  A533B  steel  are  summarized  in  Figure 
IV— 36.  It  indicates  a rather  small  rate  effect  at  low  tem- 
peratures, but  a fairly  large  effect  at  higher  temperatures. 
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Kj^and  K|Q(ksi\/lrr  ) 


FIGURE  IV-35  Versus  Loading  Rate,  K,  and  Versus 

Crack  Velocity  for  a C-Mn  Steel  Derived 
by  Eftis  and  Krafft  (1965) . (Note 
positioning  of  horizontal  scale  is  based 
only  on  a highly  approximate  model  relating 
K and  crack  velocities.) 

The  trend  observed,  from  these  investigations  is 
a generally  inverse  relationship  between  sensitivity  to 
loading  rate  effects  and.  room  temperature  yield,  strength  in 
ferritic  steels;  i.e.,  as  noted,  in  Section  B.2,  the  effects 
are  most  pronounced,  in  low  strength  steels.  Barsom  (1975), 
for  example,  concludes  that  the  sensitivity  essentially  dis- 
appearp  in  steels  with  room  temperature  yield  strength  above 
140  ksi. 

The  fracture  toughness  of  nuclear  pressure  vessel 
steels  is  observed,  to  be  loading  rate  sensitive  in  the  re— 
gion  where  valid  LEFM  toughness  values  have  been  measured. 
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FIGURE  lV-36  Influence  of  Loading  Rate,  K,  on  the  Initiation  Toughness,  K-r^. 

of  A533B  Steel.  K-r^  is  also  shown  ^ ' Id' 


Experimentaily,  this  has  meant  that,  below  the  upper  shelf 
temperature  range,  the  toughness,  Kj^^,  of  a rapidly  loaded 
specimen  containing  a stationary  crack  is  less  than  the 
toughness,  K , exhibited  by  the  same  specimen  loaded  stat- 
ically. While  the  value  of  appears  to  initially  de- 

crease with  increasing  loading  rate,  the  presence  of  a 
minimum  in  the  versus  S:  curve  has  not  yet  been  firmly 

established.  The  rising  branch  of  the  curve  for  high 
loading  rates  (cf.  Figure  IV-35)  is  based,  on  calculated, 
values  of  for  propagating  cracks,  using  an  assumed  re- 

lation between  crack  velocity  a and.  loading  rate,  K,  The 
possibility  does  exist  that  achieves,  and  remains  at, 

a minimum  value  for  high  loading  rates. 

There  is  some  evidence,  based  on  a comparison  of 
: Charpy  V-notch  slow  bend,  and  impact  energy  curves,  that  the 
dynamic  and  static  fracture  toughness  curves’,  for  loading 
rate  sensitive  materials,  may  cross  in  the  upper  shelf  tem- 
perature range.  (An  example  of  this  behavior,  for  total 
energy  absorbed,  with  A302B  steel,  is  shown  in  Figure  IV-37.) 
If  such  is  the  case,  then  the  propagation  resistance  of  load- 
ing rate  sensitive  materials  may  increase  monotonically  with 
crack  speed  in  the  upper  shelf  temperature  range,  thus  init- 
ially favoring  slow  stable  crack  growth  instead,  of  rapid 
crack  propagation,  at  these  higher  temperatures. 


FIGURE  IV-37  Slow-Bend  and  Impact  Charpy  V-Notch 
Test  Results  for  A302B  Steel 

SOURCE:  Rolfe  and  Barsom,  1977 
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Since  an  experimental  technique  is  not  available 
for  detecting  the  onset  of  crack  growth  initiation  in  most 
dynamic  tests,  it  is  generally  assumed,  that  it  occurs  at 
maximum  load.  Thus,  to  calculate  energy-based  parameters, 
the  energy  to  maximum  load  is  used.  It  is  recognized  that, 
if  the  crack  initiates  before  maximum  load,  the  elastic- 
plastic  values  generated,  in  this  way  will  exceed  the  actual 
Kjd  values  for  the  material  and,  therefore,  they  must  be 
accepted,  cautiously. 

d.  Rapid  Unstable  Crack  Propagation  and.  Crack  Arrest 

(1)  Dynamic  LEFM  Concepts.  The  problem  of  de- 
termining the  crack  length-time  response  of  a rapidly  prop- 
agating crack  is  of  great  concern  in  several  different  kinds 
of  engineering  structures.  These  structures  have  in  common 
the  feature  that  unchecked,  unstable  crack  growth  would  have 
catastrophic  consequences.  They  include  aircraft,  nuclear 
pressure  vessels,  bridges,  gas  transmission  pipelines,  and. 
ship  hulls,  to  cite  only  a few  of  those  for  which  research 
has  been  specifically  addressed. 

Although  many  actual  problems  may  require  more 
complicated  elastic-plastic  treatments,  the  process  of  rapid 
unstable  crack  propagation  and.  arrest  in  structures  cur- 
rently can  be  discussed  only  in  terms  of  linear  elastic 
fracture  mechanics  concepts  and  parameters.  A dynamic  ex- 
tension of  LEFM  recognizes  four  energy  categories:  elastic 
strain  energy,  kinetic  energy,  work  done  by  applied  forces, 
and  the  energy  dissipated  by  crack  tip  flow  and  fracture 
processes.  The  first  three  of  these  depend  on  elastic 
properties,  the  crack  length,  the  applied  loads,  and  the 
geometry  of  the  body  containing  the  crack.  The  net  change 
in  these  three  components,  per  unit  area  of  crack  extension, 
is  called  the  dynamic  energy-release  rate,  or,  equivalently, 
the  driving  force  for  crack  extension.  Giving  this  the 
symbol  G,  then 
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where  U is  the  strain  energy,  T the  kinetic  energy,  W the 
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work  done  on  the  structure  by  external  loads,  a is  the  crack 
length,  and  b is  the  plate  thickness. 

Two  differences  exist  between  the  evaluation  of  G 
in  the  static  case,  (cf . Equation  IV-16) , and  for  a fast 
propagating  or  arresting  crack.  First,  a kinetic  energy 
term  is  present.  Second,  U,  T,  and.  W must  be  evaluated  from 
fully  dynamic  analyses;  i.e,  inertia  forces  explicitly  in- 
cluded in  the  equations  of  motion  for  the  structure. 

The  energy  dissipated  per  unit  area  of  fracture 
is  called  the  fracture  energy  R.  This  expresses  the  resist- 
ance to  cracking.  The  fracture  energy  is  a material  property 
essentially  independent  of  the  geometry  and  applied  loads. 

It  is,  in  fact,  a basic  postulate  of  LEFM  that  all  inelastic 
irreversible  energy  dissipation  processes  that  accompany 
crack  extension  can  be  included  in  a single  material  prop- 
erty that  is  possibly  a function  of  the  instantaneous  crack 
speed,  but  is  independent  of  the  crack  length,  the  applied, 
loads,  and.  the  external  geometry  of  the  body.  The  extent 
to  which  this  is  true  really  determines  the  applicability 
of  LEFM. 

A crack  extension  criterion  follows  from  the  prin- 
ciple of  energy  conservation.  That  is,  the  energy-release 
rate  must  be  balanced  by  the  fracture  energy  for  the  crack 
to  propagate.  This  statement  means  that  the  initiation  of 
growth  of  a stationary  crack  and  the  continued  propagation 
of  a moving  crack  are  only  possible  when  G = R.  Conversely, 
initiation  of  crack  growth  is  precluded — and,  for  a propa- 
gating crack,  arrest  must  take  place — when  G < R for  all 
values  of  R. 

Note  that  the  condition  where  G exceeds  R is  not 
possible  because  it  would  violate  the  energy  balance  prin- 
ciple. The  stationary  crack  inequality  does  not  violate  the 
energy  balance  because  the  basic  energy  balance  is  actually 
Equation  IV-28  multiplied  by  a.  Note  also  that  in  LEFM  the 
fracture  energy  for  the  extension  of  a stationary  crack  due 
to  slowly  applied,  loads,  designated,  as  G^,  corresponds  to  a 
critical  value  of  the  static  energy-release  rate. 

The  arrest  of  a rapidly  propagating  crack  in  a 
structure  under  load  can  be  considered  on  several  different 
levels  of  complexity.  Starting  from  the  simplest  (and.  least 
accurate)  and  continuing  with  more  complicated  and  more  ac- 
curate approaches,  the  various  types  can  be  classified  as 
either  a; 

• Static  or  quasi-static  analysis 
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• Infinite  medium  dynamic  analysis 

• Fully  dynamic  analysis 

The  primary  distinction  that  differentiates  between  static 
and  dynamic  approaches  is  that  specimen  inertia  terms  and, 
the  contribution  of  kinetic  energy  to  the  crack-driving 
force  equation  are  excluded  in  the  former.  Physically,  this 
means  that  static  or  quasi-static  theories  are  limited  to 
situations  where  (1)  the  loading  is  static  or  applied  slowly, 
(2)  the  crack  propagates  slowly  and.  it  (3)  changes  speed 
only  gradually.  As  the  extensive  work  done  at  Battelle 
(Hahn,  et  al.,  1975-1976)  and  other  facilities  has  shown, 
the  arrest  of  rapid  crack  propagation  tends  to  occur  rather 
abruptly.  This  alone  indicates  that  statically  based  treat- 
ments must  be  applied  to  crack  arrest  with  due  caution. 
Quantitative  results  reinforcing  this  idea  have  also  been 
produced . 

The  distinction  between  the  two  dynamic  analysis 
procedures  lies  in  the  particular  specialization  that  is 
involved  for  simplification.  By  considering  the  structure 
to  be  an  infinite  elastic  medium,  the  effect  of  stress  waves 
reflected  back  to  the  propagating  crack  tip  from  the  bound- 
aries of  the  structure  and/or  from  internal  load  points  and. 
discontinuities  (e.g.,  welded-on  stiffeners)  are  neglected. 
These  effects  are  taken  into  account  in  a fully  dynamic  an- 
alysis, albeit  for  practical  reasons  at  the  expense  of 
specializing  the  structural  geometry  under  consideration. 

(2)  Dynamic  Crack  Propagation  Analysis  Methods. 

(a)  Basis  of  Dynamic  Crack  Propagation 
Analyses . There  is  currently  no  universally  accepted  theo- 
retically-based. design  approach  to  ensure  crack  arrest.  A 
static  approach  (or,  what  amounts  to  the  same  thing,  the 
"arrest  toughness",  approach)  is  not  generally  valid. 

This  statement  is  based,  upon  a body  of  experimental  results 
together  with  a rigorous  energy-based  method,  of  anlaysis 
that  has  shown  that  crack  arrest  is  a dynamic  process  that 
must  be  treated,  within  the  context  of  a dynamic  fracture 
mechanics  theory.  This  work  includes  the  generalizations 
of  the  static  approach  discussed  above. 

One  purpose  of  this  subsection  of  the  report  is  to 
demonstrate  that  statically  based  analyses  can  dangerously 
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overestimate  the  capacity  of  a structure  to  arrest  a rapidly 
propagating  crack.  It  will  be  further  shown  that  even  anal- 
yses taking  full  account  of  the  essential  aspects  of  dynamic 
fracture  mechanics  can  still  be  inadequate  for  predicting 
crack  arrest  when  other  vital  features  of  the  problem  are 
neglected.  In  particular,  the  otherwise  admirable  analyses 
of  Freund  (1972,  a,  b)  cannot  cope  with  stress  waves  re- 
flected back  to  the  propagating  crack  tip  from  the  specimen 
boundaries.  This  feature  is  not  only  important  for  experi- 
mental work  carried,  out  using  small  laboratory  size  test 
specimens,  but  might  also  be  important  in  analyzing  other 
mechanical  crack  arrest  devices. 

It  is  a fact  that  the  arrest  point  determined  from 
either  a static,  a quasi-static,  or  an  infinite  medium 
dynamic  treatment  will  always  be  closely  related  and,  in 
some  cases,  will  be  exactly  the  same.  Consequently,  for  the 
purpose  of  this  report,  it  will  suffice  to  describe  the  most 
accurate  of  these  with  a view  toward  contrasting  its  crack 
arrest  predictions  with  those  of  a fully  dynamic  calculation. 

(b)  Crack  Propagation  in  an  Infinite  Medium. 

A very  elegant  analysis  of  the  propagation  of  a semi-infinite 
crack  in  an  infinite  medium  is  that  given  by  Freund.  (1972, 
a,bf  see  also  Freund,  (1976)  ) . Using  a Laplace  transform  in 
conjunction  with  the  Wiener-Hopf  Technique,  Freund  has  solved 
the  equations  of  motion  for  a half-plane  crack  propagating 
in  an  unbounded  medium  for  a fairly  unrestricted,  class  of 
crack  motion.  A key  result  of  the  analysis  relates  the 
dynamic  stress  intensity  factor  K,  which  is  a function  of 
instantaneous  crack  length  a and.  speed,  a to  the  product  of 
the  static  intensity  factor  Kg  and.  a universal  function 
k(6')  of  instantaneous  crack  speed  relative  to  the  elastic 
wave  speeds  according  to 


K(a,a)  = k(a)  Kg  (a) 


(IV-29) 


The  function  k decreases  monotonically  from  unity  at  zero 
crack  speed  to  zero  at  the  Rayleigh  wave  speed.  The  quan- 
tity Kg  is  a "static"  factor  in  that  it  corresponds  to  zero 
crack  speed;  however,  K„  may  result  from  dynamic  loading  of 

o 

the  type  described  by  Freund  (1976) . 
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A second  key  result  obtained  by  Freund  is  one 
that  relates  the  dynamic  energy-release  rate  to  the  opening- 
mode  dynamic  stress-intensity  factor.  For  plane  strain  con- 
ditions, this  is 


G (a  , a) 


( a,  *a) 


(IV-30) 


where  B is  also  a geometry  independent  universal  function  of 
the  crack  speed  relative  to  the  elastic  wave  speeds.  In 
contrast  to  k,  it  is  a monotonically  increasing  function 
which  is  unity  at  zero  speed,  and  becomes  unbounded  at  the 
Rayleigh  speed. 

It  is  important  to  recognize  that  Equation 
IV-30  is  valid  for  arbitrary  geometries  and  time-variation  of 
crack  speed.  (Freund,  1976)  . As  a result,  it  is  possible  to 
use  the  idea  of  dynamic  fracture  toughness  inter- 

changeably with  an  intrinsic  material  energy  dissipation  rate, 
R(a),  (also  called  "specific  fracture  energy"  or  simply 
"fracture  energy") . That  is,  let  us  define  K (a)  through  the 
equation  ^ 


R(a)  = 


1-v^ 

E 


(a ) 


(lV-31) 


Then,  by  equating  R and  G there  results 


BK^  ( a,  k)  = ( k) 


(IV-32) 


which  serves  as  the  crack  propagation  criterion;  i.e.,  as  an 
equation  for  predicting  a and  a.  It  is  also  found,  from 
Equations  IV-2  9 and.  IV-32  that 


Kg  (a)  = Kj(k)/g(k) 


(IV-33) 


where  g^Bk^  is  also  a universal  function  of  instantaneous 
crack  speed.  The  function  g = g ( a)  can  be  interpreted  as 
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the  ratio  of  the  dynamic  to  the  static  energy  release  rates. 

In  order  to  apply  Equation  IV-33  to  investi- 
gate crack  arrest,  an  explicit  relation  for  the  function 
g(^)  is  needed.  In  Freund's  analysis,  a numerical  integra- 
tion was  used  to  determine  this  function.  His  prediction 
reveals  that  the  function  g = g(^)is  more  than  adequately 
expressed  by  the  simple  relation 


g{a)  = 


(IV- 34) 


where  denotes  the  Rayleigh  wave  speed.  Then,  by  substi- 
tuting Equation  IV-34  into  Equation  IV-33,  the  equation  of 
motion  for  the  crack  tip  becomes 


Kg  (a  ) = Kj^(a) 


1 


(lV-35) 


The  next  step  is  to  introduce  a relation  for  Kg  for  the 
geometry  of  interest.  Given  Kj^(a),  Equation  IV-35  can  then 
be  solved  iteratively  for  the  crack  speed  as  a function  of 
crack  length.  Bu  numerically  integrating  this  result,  the 
crack  length  can  be  obtained  as  a function  of  time  for  com- 
parison with  results  .computed  with  a fully  dynamic  approach. 
This  comparison  is  made  below. 

A propagating  crack  will  arrest  when  ^nd  only 
when)  Equation  IV-33  can  no  longer  be  satisfied.  Suppose 
the  right-hand  side  of  Equation  IV-33  possesses  one  minimum 
point  and  that  it  occurs  at  the  speed,  ajyj  (which  may  be  zero)  . 
Then,  the  arrest  occurs  at  the  crack  length  ® , as  calcu- 
lated from 


K|(ar)  = • {IV-36) 


Equation  IV-35  shows  that  if  1,  crack  arrest  occurs 

jyl 

at  the  speed  for  which  is  a minimum;  viz 
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min 


(IV-37) 


(c)  An  Analysis  Procedure  for  the  DCB  Test 
Specimen . A laboratory  test  specimen  that  has  been  used 
effectively  by  Hahn,  et  al.,  (1975-1976)  is  the  double  canti- 
lever beam  (DCB)  specimen  (cf . Figure  IV-40) . A model  for 
predicting  dynamic  crack  propagation  in  the  DCB  specimen 
has  been  given  by  Kanninen  (1974) . The  starting  point  for 
the  derivation  is  the  equations  of  the  theory  of  elasticity 
with  inertia  terms  included.  Because  the  peculiar  "beam- 
like" geometry  of  this  specimen  can  be  exploited,  only  four 
equations  need  to  be  explicitly  considered.  These  are  the 
two  equations  for  motion  along  the  length  of  the  beam  (the 
X direction)  and  normal  to  the  crack  plane  (the  z direction) 
and  two  Hooke's  law  equations.  The  two  equations  of  motion 
are  given  by 


and 


Bt 


xz 

Bz 


Bt 


xz 

Bx 


Bt 


+ 


_yz 


By 


P 


B^u 

z 

Bt^ 


(lV-38) 


(lV-39) 


The  two  constitutive  or  Hooke's  law  equations  that  enter 
into  the  analysis  are  given  by 


and. 


Bu 


E- 


Bx 


Bu  Bu 

X z 

Bz  ^ Bx 


T 

XZ 

G 


(IV-40) 


(lV-41) 


In  the  above  equations,  E,  G,  v,  and.  p denote  Young's 
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modulus,  the  shear  modulus,  Poisson's  ratio,  and  the  density, 
respectively;  u and  u are  displacement  components;  a » a / 


T , T , 
xy  xz 


Ty2  itress  components;  and  t denotes^tim^. 


The  simplification  that  can  be  introduced,  to 
make  the  mathematical  analysis  more  manageable  is  made  by 
introducing  cross-sectionally  averaged  dependent  variables 
into  the  analysis.  If  A = A(x)  is  the  area  of  the  DCB 
specimen  cross  section  at  any  axial  position  x,  then  these 
new  variables  can  be  obtained,  formally  as  follows.  The 
deflection  w = w(x)  is 


(IV-42) 


The  rotation  i|f  = ij;  (x)  is 


■ - \U 


zu  dydz 

X 


(IV-43) 


// 


By  operating  on  Equations  IV-38  through  IV-41  by^./  dydz 
and  using  the  above  definitions,  the  equations  of  motion  for 
a rectangular  DCB  specimen  in  terms  of  the  cross-sectionally 
averaged,  variables  for  the  case  of  fixed  displacement  loading 
are  found  to  be 


5^w ^ H(x-a)  w = (IV-44) 

3x^  ’ 5x  h®  Cq  Bt^ 


and. 


Bx^ 


n Cq 


lli 

Bt® 


(lV-45) 


where  h is  the  half-height  of  the  specimen,  b is  the  specimen 
thickness,  and  Cq  = E/p  is  the  elastic  bar  wave  speed.  The 
function  H denotes  the  Heaviside  step  function  which  delin- 
eates the  position  of  the  crack  tip  (i.e.,  x=a ) . 
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Expressions  for  the  strain  energy  and  the 
kinetic  energy  of  the  system  can  be  obtained  in  terms  of 
the  variables  introduced  above.  Omitting  the  details,  the 
resulting  expression  for  the  strain  energy  U is 


(lV-46) 


while  the  kinetic  energy  T is 


(lV-47) 


where  L is  the  overall  length  of  the  specimen.  The  most 
important  use  of  the  strain  and  kinetic  energy  expressions 
is  in  determining  the  crack-driving  force.  This  is  done 
through  the  definition  of  the  dynamic  energy-release  rate 
G in  terms  of  an  energy  balance  for  the  system  as  given  by 
Equation  IV-28.  By  substituting  Equation  IV-46  and  IV-47 
into  IV-28,  it  is  found  that  G can  be  interpreted  in  terms 
of  "crack  tip"  values.  That  is 


G 


2E 

h 


,^4 

x=  a ( t ) 


(lV-48) 


where,  as  indicated.,  the  bracketed  quantity  is  evaluated,  at 
the  axial  position  representing  the  current  crack  tip. 

As  will  be  explained  below,  tests  in  the  DCB 
specimen  are  conducted,  by  blunting  the  initial  crack  tip. 
This  allows  enough  elastic  strain  energy  to  be  stored  in  the 
specimen  that  the  crack  (which  propagates  as  a sharp  crack) 
can  propagate  at  a high  speed.  The  measure  of  the  bluntness 
ks  Kg,  the  "apparent"  stress-intensity  factor  at  the  initia- 
tion of  crack  growth.  An  example  set  of  computational 
results  is  shown  in  Figures  IV-38  and  IV-39  for  the  case 
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FIGURE  IV-38  Comparison  of  Crack  Arrest  Points 

Predicted  by  a Fully  Dynamic  Analysis 
with  that  of  a Quasi-Dynamic  Analysis 
for  a Standard  Rectangular  DCB  Specimen 

SOURCE:  Kanninen  et  al,  1977 

where  Kg  = 2Kj^  and  = K^^., 

It  can  be  clearly  seen  in  Figure  IV-38  that 
the  crack  arrest  point  given  by  the  solution  for  an  infin- 
ite medium  is  considerably  underestimated.  Ihat  is,  while 
the  fully  dynamic  theory  coincides  with  the  infinite  medium 
solution  until  the  time  of  the  first  stress  wave  reflection 
(approximately  30  fi  sec  in  this  example)  , sizable  differ- 
ences can  be  seen  thereafter. 

Consideration  of  the  results  shown  in  Figure 
IV-39  reveals  that  a statically  calculated  arrest  point 
(i.e.,  a - a^ = 37  mm)  is  reached  at  about  the  same  point 
that  kinetic  energy  reaches  a maximum  in  the  fully  dynamic 
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FIGURE  IV-39  Distribution  of  Energy  During  Rapid 

Crack  Propagation  in  a DCB  Test 

Specimen  for  K /K  2.0  and  K = K 

q IC  D IC  • 

SOURCE:  Kanninen  et  al,  1977 

calculation.  That  the  average  rate  of  change  of  the  .kinetic 
energy  is  greater  (negatively)  than  the  strain  energy  after 
the  maximum  has  been  reached,  further  shows  that  the  kinetic 
energy  provides  the  greater  contribution  to  the  crack-driving 
force.  Equation  IV-28. 
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In  a recent  analytical  study  by  Freund  (1977) , 
the  arms  of  the  DCB  specimen  are  assumed  to  deform  as  shear 
beams.  This  idealization  yields  a mathematically  tractable 
model  which  predicts  qualitatively  correct  behavior  and. 
lends  itself  to  a simple  physical  interpretation  of  the 
specimen  response.  It  is  found  that  crack  arrest  occurs 
when  the  tip  is  overtaken  by  a stress  wave  which  is  set  up 
• at  fracture  initiation  and  reflects  off  the  top  end  of  the 
DCB  (cf . Figure  IV-40) . 


FIGURE  IV-40  DCB  Test  Piece  Configuration: 

Side  Grooves  are  not  Shown 

SOURCE:  Hahn  et  al,  1975-76 
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(3)  Dynamic  Crack  Propagation  Testing  Methods. 
Certain  problems  exist  in  the  measurement  of  the  high  levels 
of  fracture  toughness  exhibited  by  many  useful  engineering 
materials.  The  reason  is  that  the  most  highly  developed 
methods  for  measuring  fracture  toughness  are  restricted  to 
fractures  in  which  plasticity  is  limited.  Specifically, 
only  fractures  that  occur  under  plane  strain  conditions  or 
under  plane  stress  conditions  in  which  the  plastic  zone  size 
is  small  relative  to  specimen  dimensions  and,  crack  length 
(see  ASTM-E-399-74)  can  be  approached  in  this  way. 

As  may  be  clear  from  the  preceding  discussion,  it 
is  necessary  to  distinguish  among  several  different  fracture 
toughness  parameters.  It  is  generally  agreed  that  .the  most 
convenient  to  use  in  material  testing  work  are  the  following: 

Criterion  for  onset  of  crack  extension:  K = K 

c 

Criterion  for  continuing  propagation:  K = 

Criterion  for  crack  arrest:  K < K 

D,min 

The  plane  stress  values  of  tough  materials  display  a modest 
thickness  dependence,  , where  b is  the  thickness  and 

0.25  .<  n < 1.0.  When  plasticity  is  limited,  fractures  can 
be  analyzed  by  the  methods  of  LEFM  to  obtain  plane  strain 
fracture  toughness  parameters  such  as  K , K , K , and 
^la'  their  plane  stress  counterparts^with  the^?  removed 

from  the  subscript.*  However,  the  limited  plasticity  at 
fracture  that  is  necessary  for  successful  application  of 
LEFM  methods  is  the  very  antithesis  of  the  desired  behavior 
of  tough  materials  where  large  plastic  zones  and  significant 
shear  lips  are  essential  to  proper  performance.  Accordingly, 
problems  arise  in  attempting  to  use  LEFM  methods  for  meas- 
uring fracture  toughness  parameters  of  tough  materials. 


* The  subscript  I (i.e.,  g k ) is  used  in 

this  subsection  to  distinguish  energy  and  toughness  values 
measured  when  the  crack  tip  plastic  flow  is  predominantly 
plane  strain  as  opposed,  to  so-called  "plane  stress"  values 
which  reflect  significant  amounts  of  through- the-thickness 
deformation. 
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In  the  following  paragraphs  several  methods  for 
measuring  or  approximating  values  for  tough  steels  by 
the  methods  of  fracture  rtechanics  are  described. 

(a)  Approximating  K Values  with  K , K , or 
~IC  . When  coincides  with  K ^the^S  - 

value  at  zero  velocity)  , static  measurements  of*^Kc  can  Serve 
as  a conservative  estimate  of  K The  K -values  can,  in  turn, 
be  related,  to  the  plane  strain  R^^via  a relation  of  the  type 


K = CK 

Ic 


(IV-49) 


where  1 < C <2.  This  means  that  K -values  could  also 
serve  as  a lower  bound  measure  of  K^'^or  as  a way  of  esti- 
mating provided  the  factor  C is  Snown.  In  practice,  the 
plate  size  and  thickness  requirements  for  measuring  K - and 
Kj^-values  for  materials  with  ^ > 0.2m ^ are  prohibitive. 

CTy 

For  example,  for^a  steel  with  a yield  strength  a = 275  MNm”^ 
and  = 300  MNm"'^/  , the  width  of  a center  cracked  panel 
adequate  to  measure  is  about  3m,  and  the  thickness  re- 
quired to  measure  is  about  Im.  More  recent  J — tech- 

niques offer  the  possibility  of  reducing  the  thickness 
requirement  by  an  order  of  magnitude.  Consequently,  j - 
measurements  may  offer  one  practical  route  to  the  evaluation 
of  - or  R values  of  high  toughness  materials.  However, 
the  path  independence  of  the  J-integral  does  not  hold,  for 
dynamic  situations.  This  means  that  the  conventional  inter- 
piretations  ina.y  b©  invo^lid.  in  thos©  cases , 

(b)  Approximating  K from  Crack-Openina 
Displacement.  Robinson  and  Tetelman  (1976)  have  shown  that 
^Ic  calculated  from  measurement  of  the  crack-tip 

opening  displacement  6 (COD)  at  the  onset  of  unstable  frac- 
ture in  relatively  small  specimens,  using  a relation  that 
can  be  obtained,  from  Equations  IV-20a  and.  IV-2  3b: 
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Methods  for  measuring  COD  are  described  in  British  Standards 
DD19:1972.  Use  of  such  COD  techniques  would  then  permit 
to  be  approximated  by  this  relation. 

The  possibility  of  applying  COD  methods  to 
dynamic  tests  also  exists.  Actual  COD  measurements  are  dif- 
ficult but,  as  Robinson  and  Tetelman  have  shown,  COD  values 
can  be  approximated  reasonably  well  from  measurements  of  the 
notch  root  contraction  on  the  fractured  test  piece. 

(c)  Direct  Measurement  of  or  with  the 

Duplex-Double-Cantilever  Beam  Tests.  Although  several  kinds 
•of  specimen  geometries  can  be  used  for  obtaining  the  data 
necessary  to  determine  K -velocity  relations,  the  wedge- 
loaded  DCB  specimen  has  been  predomimately  utilized.  (The 
.rectangular  DCB  specimen  and  wedge-loading  arrangement  are 
illustrated  in  Figure  IV-40.)  In  this  specimen,  the  frac- 
ture is  initiated  with  the  aid.  of  a blunt  slot  rather  than 
a sharp  crack.  The  blunted  notch  permits  the  specimen  to 
sustain  an  "apparent"  stress  intensity  factor,  K^,  which 

is  much  larger  than  K . Consequently,  as  soon  as  a sharp 
cj^ack  emerges  from  the  blunt  notch,  the  crack  immediately 
becomes  unstable  and  propagates  rapidly.  The  stress  inten- 
sity at  the  onset  of  fast  fracture  can  be  systematically 
altered  by  varying  the  slot  root  radius.  Typically,  root 
radii  of  0.02  to  0.08  inch  have  been  used.  The  K-values  at 
initiation  and  arrest  are  measured  with  a clip  gage  fixed  to 
the  end  of  the  test  piece.  The  velocity  during  propagation 
is  measured,  with  conducting  strips  deposited  over  a thin 
insulating  layer  of  epoxy. 

The  specimen  is  slowly  loaded  in  an  ordinary 
testing  machine  by  forcing  a split  wedge  between  the  pins. 
Since  the  wedge  loading  is  inherently  stiff,  crack  propa- 
gation proceeds  with  essentially  constant  displacement  at 
the  load  point.  Under  quasi-static  conditions,  the  crack- 
driving  force  would,  decrease  as  the  crack  grows.  This 
behavior  ultimately  causes  a high-speed  crack  to  arrest 
within  the  confines  of  the  specimen,  provided  the  test  piece 
is  long  enough. 

Wedge  loading  has  two  other  virtues.  One  is 
. that,  since  little  energy  is  exchanged  between  the  DCB 
specimen  and,  the  testing  machine  during  the  propagation 
, event,  the  results  are  relatively  insensitive  to  the  char- 
acter of  the  testing  machine.  The  second  is  that  the  fric- 
tion between  the  wedge  and  the  pins  introduces  a modest  com- 
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pressive  stress  parallel  to  the  crack  plane.  This  tends  to 
stabilize  the  crack  path.  Hence,  the  side  grooves  ordinar- 
ily required  to  keep  the  crack  from  turning  can  be  omitted 
for  those  materials  in  which  the  crack  has  some  tendency  to 
branch. 

Although  the  experimental  method  outlined 
above  has  proven  satisfactory  for  testing  high  strength 
materials,  several  modifications  have  been  found  necessary 
for  testing  the  lower  strength  grades  of  steels.  One  mod- 
ification is  the  use  of  a so-called  "duplex"  specimen  where 
a high  strength  material  is  used,  for  the  starter  section  to 
avoid  excessive  plastic  deformation  in  the  loading  regions 
of  the  specimen.  A starter  section  of  high  strength  4340 
steel  welded,  to  the  test  material  has  been  successfully 
employed  for  fabricating  duplex  specimens.  Another  problem 
has  been  a tendency  for  crack  branching  in  the  lower  strength 
steels.  This  has  necessitated  side  grooving  to  maintain  a 
straight  crack  plane.  Currently,  a 60-percent  side  groove 
depth  duplex  specimen  is  being  used  for  tests  on  A533b  steel. 

Some  typical  crack  length  time  results  are 
shown  in  Figure  IV-41.  Of  interest  is  the  fact  that  these 
are  ostensibly  linear,  which  gives  considerable  credence  to 
the  use  of  the  mathematical  model  described  in  Section 
B.4.d(2) (c) (cf.  Figure  IV-38) . This  linearity  also  facili- 
tates the  interpretation  of  the  results  since  a uniform 
crack  speed  value  can  be  associated  with  the  results  of  each 
test.  Thus,  using  the  results  of  the  dynamic  analysis  of 
the  DCB,  there  are  three  methods  of  extracting  Kjj^j  values 
from  the  two  independent  quantities  measured,  on  the  test 
specimen:  the  crack  velocity  and  the  crack  length  at  arrest. 

These  are 

(1) 


(2) 


(3) 


Obtaining  by  matching  the  initial 

average  velocity  with  calculated 
velocities  for  various  K , 
specimen  parameters.  ^ 

Obtaining  K - by  matching  the  arrest 

crack  length  with  calculated  values. 

Evaluating  an  average  value  of  K (K  ) 

W-V,  a.'  ID  ID'^ 

over  the  entire  propagation  event 
through  the  approximation 
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(lV-51) 


"'id  = ■ 

which  is  based  primarily  on  the 
assumption  that  no  kinetic  energy 
is  left  in  the  specimen  at  arrest. 


TIME  (/jsec) 

FIGURE  IV-41  Assorted  Crack  Length-Time  Measurements 

from  Wedge-Loaded,  Rectangular  DCB  Specimen 
Tests 

SOURCE:  Hahn  et  al,  1975-76 

For  Method  2 (and  also  for  Method  1 to  some  degree)  an  iter- 
ative determination  is  necessary  since  the  calculated  values 
depend  on  Kjj^-  velocity  relations  which  are  not  known  before- 
hand for  the  material  being  tested.  In  the  case  of  duplex 
specimens,  some  modifications  of  each  of  the  three  methods 
are  required,  but  the  basic  principles  remain  the  same. 

There  are  several  reasons  why  valid  Kjj^  data 
can  be  obtained  from  relatively  small  duplex  DCB  specimens 
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of  moderately  tough  steels.  First,,  the  high— strength/low— 
toughness  starter  section  reduces  the  plane  strain  thickness 
requirements  drastically.  Second,  the  grooves  along  the 
fracture  path  develop  constraints  similar  to  those  associa- 
ted with  increased  plate  thickness.  Third,  the  propagating 
crack  produces  a very  high  strain  rate  at  the  crack  tip. 

This  causes  the  effective  yield  strength  to  be  raised 
(perhaps  twice  the  static  yield  strength) , hence  reducing 
the  plastic  zone  size. 

In  principle,  the  duplex  DCB  test  can  be  also 
used  to  obtain  valid  Kjj  data  by  eliminating  the  side  grooves. 
However,  for  the  specimen  dimensions  currently  employed  to 
measure  Kjj),  removal  of  the  side  grooves  would  introduce 
several  problems,  particularly  for  high  toughness  steels. 

The  amount  of  strain  energy  that  can  be  stored  in  the  arms 
of  the  test  piece  is  not  sufficient  to  drive  the  crack  into 
the  tough  test  section  for  any  appreciably  distance,  thus 
making  analysis  of  the  results  difficult.  With  existing 
specimen  dimensions  and  procedures,  the  estimated  upper  limit 
of  toughness  measurements  is  about  250  MNm  3/2 _ Also,  the 
plastic  zone  radius  may  approach  or  exceed  the  arm  height  of 
the  specimen.  Finally,  cracks  propagating  into  the  tough 
test  section  frequently  branch  in  the  absence  of  side 
grooves;  this  may  preclude  analysis  of  the  results. 

(4)  Material  Property  Data  tor  Dynamic  Crack 
Propagation.  The  treatment  of  dynamic  crack  propagation 
and  arrest  is  complicated  by  the  variation  of  the  fracture 
energy  with  crack  velocity  and  plate  thickness.  A schematic 
of  velocity  dependence  is  shown  in  Figure  IV-42a.  Eftis  and 
Krafft  (1965)  have  deduced  values  from  wide  plate,  ship 

steel  experiments.  Their  results  reflect  low  energy  cleav- 
age fractures  below  the  nil  ductility  temperature  (NDT) . 

These  results  indicate  that  first  decreases  with  increas- 

ing velocity,  becomes  a minimum  at  a finite  velocity,  and. 
then  increases  dramatically  for  crack  velocities  in  excess 
of  600  ms.  In  contrast,  recent  results  for  low  energy 
fibrous  fractures  in  AISI  4340  steel  are  reproduced  in 
Figure  IV-42b.  Here  the  fracture  energy  increases  mono- 
tonically  with  crack  velocity. 
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Fracture  Toughness 


with  a minimum  showing  K 


min 


Fraction  of  Bar  Wave  Speed,  a/C 

o 

0 0.05  0.10  0.15 


(b)  Results  for  flat  fibrous  fracture 
of  4340  steel  after  Hahn,  et  al. 

FIGURE  IV-42  Examples  of  the  Crack  Velocity- 
Dependence  of  the  Propagating 
Crack  Fracture  Toughness 

SOURCE;  Hahn  et  al,  1975-76 
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Figure  IV-43  shows  measured  values  of  K obtained 
from  duplex  DCB  specimens  for  several  grades  of  s5ip  steel 
tested,  near  the  nil  ductility  temperature.  The  behavior  of 
these  steels  is  interesting  from  several  standpoints.  The 
toughness  of  each  of  the  four  steels  is  seen  to  be  strongly 
dependent  on  crack  velocity,  being  greatest  at  small  veloc- 
ities. Nonetheless,  for  a propagating  crack  exceeds  K , 

the  energy  associated,  with  crack  initiation  by  impact. 

There  is  some  thought  that  the  minimum  in  the  K -versus- 
velocity  curve  (if  a minimum  exists)  may  approximate  the 
values,  but  this  has  not  yet  been  demonstrated. 

values  have  been  determined  for  a limited 
number  of  low  to  intermediate  strength  grades  of  steel.  In 
the  following  discussion  of  these  data,  and  the  accompanying 
figures,  rather  than  is  used  to  denote  the  measured 

fast  fracture  toughness  since  it  is  not  certain  that  the 
measured  values  are  fully  plane-strain  values. 


Crack  Velocity,  m/s 


FIGURE  IV-43  Relation  Between  Crack 

Velocity  and  Dynamic  Toughness 
for  Ship  Steels  Tested  Near  NDT 


SOURCE;  Kanninen  et  al,  1977 
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Figure  IV-44  shows  results  obtained  for  a 9 per- 
cent Ni  steel  and  an  A517F  steel  at  -196  C.  It  can  be  seen 
that  the  K values  of  the  9 percent  Ni  steel  are  nearly 
independent  of  crack  velocity  while  the  A517F  seems  to  go 


FIGURE  IV-44  Fast  Fracture  Toughness  Data  for  a 
9 Percent  Ni  and  an  A517-F  Steel 

SOURCE:  Battel le  Columbus  Laboratories,  unpublished  data 
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through  a slight  minimum  and  then  an  increase.  At  this 
temperature,  the  fracture  modes  are  ductile  dimple  for  the 
9 percent  Ni  steel  and  cleavage  for  the  A517F  steel.  Data 
for  low  carbon  steels  in  the  temperature  range  from  -58  F 
to  - 10  F are  shown  in  Figure  IV-45 . 
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FIGURE  IV-45  Fast  Fracture  Toughness  Data  for  Low  Carbon 
Steels  at  Low  Temperatures 

SOURCE:  Battelle  Colombus  Laboratories,  unpublished  data 


Figure  IV— 46  shows  K values  as  a function  of 
temprature  (relative  to  the  i^il  ductility  temperatures, 
NDT)  and  several  crack  velocities  for  three  ship-plate 
steels  (ABS  grade  C,  E,  and  EH) . These  are  basically  C-Mn 
steels  with  room  temperature  yield  strengths  in  the  40  to 
55  ksi  range.  At  temperatures  around  and  above  the  NDT,  a 
rather  large  and.  continuous  decrease  in  Kj^  with  crack  veloc- 
ity IS  indicated  for  these  steels.  Finally,  Figure  IV-47 
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shows  some  preliminary  values  for  A533B  steel  measured  in 
a current  research  program.  The  investigators  caution  that 
these  are  tentative  results  which  may  have  to  be  corrected 
pending  further  analysis. 


FIGURE  IV-46  Fast  Fracture  Toughness  Data  for 

Several  Grades  of  Ship-Plate  Steels 

SOURCE:  Battelle  Columbus  Laboratories,  unpublished  data 
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FIGURE  IV-47  Preliminary  Data  for  A533B  Steel 
SOURCE:  Battelle  Columbus  Laboratories,  unpublished  data 

Another  avenue  proposed  for  evaluating  large  frac- 
ture toughness  values  is  to  measure  the  total  energy  absorb- 
ed in  the  fracture  of  a notched  bend  specimen  in  the  dynamic 
tear  (DT)  test.  The  energy  to  fracture  the  specimen  is 
provided  by  a pendulum  whose  velocity  just  prior  to  impact 
is  approximately  5 to  10  m/s.  DT,  the  total  energy  absorbed 
in  the  process  of  breaking  the  specimen,  is  termed  the 
dynamic  tear  energy.  It  is  observed  directly  by  noting  the 
height  of  the  pendulum  swing  after  fracture.  This  quantity 
divided  by  A,  the  cross  sectional  area  of  the  test  piece. 
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is  related  to  the  fracture  energy  R assuming  a constant 
crack  velocity.  That  is 


R 


DT 

A 


(lV-52) 


where  /3  = 1 when  the  energy  losses  in  the  impact  test  remote 
from  the  crack  tip  are  zero  and  )3  > 1 when  significant  energy 
losses  occur.  The  dynamic  fracture  toughness  can  then  be 
expressed  as 


= 


(lV-53) 


if  plane  stress  is  assumed.  For  steels,  the  DT  versus 
plots  indicate  that  (3  10  for  the  1-inch  DT  test  and 

for  the  5/8-inch  DT  test  if  we  assume  = K . Limited  DT 

versus  K data  suggest  that  |8  is  only  about  1.4,  but  addi- 
tional data  are  required,  to  confirm  this. 

The  experimental  observation  that  /3  is  greater 
than  1.0  confirms  the  assumption  generally  made  about  the 
DT  test:  that  the  energy  losses  remote  from  the  crack  tip 
in  an  impact  test  are  of  a significant  magnitude  and,  in 
some  cases,  can  overwhelm  the  actual  fracture  propagation 
energy.  Included  in  these  losses  are  the  crack  initiation 
energy,  the  energy  associated  with  plastic  deformation  at 
the  loading  points  and  at  the  specimen  boundaries  as  the 
crack  approaches  the  far  side  of  the  test  bar,  and  the  kin- 
etic energy  of  the  fractured  specimen.  The  experimental  data 
are  too  meager  to  permit  estimation  of  reliable  (8  values. 


e . Fracture  Criteria  for  Moving  Cracks  and  Complex 
Loading. 


(1)  Path  Independent  Integral.  Parameters  such 
as  or  must  be  clearly  distinguished,  from  the  specific 

surface  energy  quantity* , T (t) , frequently  used  in  the 


*This  surface  energy  is  equal  to  one-half  the  fracture 
energy  R(a)  for  a propagating  crack  when  the  conditions  of 
linear  elastic  fracture  mechanics  are  met  (cf.  Equation 
IV-31)  . 
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analytical  modeling  of  moving  cracks  in  which  r(t)  describes 
the  energy  dissipated,  in  the  creation  of  free  surface  as  a 
function  of  time  and  is  not  necessarily  associated  with  any 
point  of  crack  instability.  For  a propagating  crack,  an 
energy  inequality  may  be  written  in  the  foritf 


/ 


n 

T . u , d s 
1 1 


pu . u , dA 
11 


R 


f ^ 

R 


(t)  *a  r{lV-54) 


if  temperature  changes  due  to  straining  are  neglected. 

(except  for  those  in  a thin  layer  next  to  the  new  surface). 
The  above  equation  stateg  that  the  rate  at  which  work  is 
being  done  by  tractions  T.  acting  across  the  contour  C 
surrounding  the  crack  tip  (Figure  IV-48)  is ’greater  than  or 
equal  to  the  rate  of  increase  of  stored,  and.  kinetic  ener- 
gies inside  C plus  the  rate  at  which  energy  is  dissipated 
by  the  moving  crack.  In  Equation  IV-54,  a is  the  crack 
velocity,  U is  the  Helmholtz  free  energy  per  unit  mass,  p 
is  the  mass  density,  and  the  dot  represents  the  differen- 
tiation with  respect  to  time.  The  components  of  the  dis- 
placement vector  are  u, . 

To  reiterate,^  the  specific  energy  r(t)  in  Equation 
IV-54  should  be  clearly  distinguished  from  the  fracture 
initiation  quantity,  2 = Jjc'  used  for  analyzing  stationary 
cracks  that  become  unstable  upon  a small  Segment  of  crack  ex- 
tension. While  can  be  determined  from  a Kj^  or  type 

of  test  and  used  as  a fracture  toughness  quantity,  r(t)  simply 
represents  the  energy  dissipated  in  the  continuous  creation 
of  new  crack  surface.  A reliable  method  of  determining  r(t) 
is  lacking  when  large-scale  olastic  deformation  exists. 

The  form  of  Equation  IV-54  may  be  cast  into  a con- 
venient form  for  the,  evaluation  of  r(t)  by  imposing  the 
following  assumptions  on  the  physical  system; 

(a)  The  crack  is  restricted,  to  propagate 
along  the  x-axis  (Figure  IV-50)  in  a steady  state  manner; 
i.e.,  the  stress  and  displacement  distributions  are 
assumed,  independent  of  time  when  referred,  to  a coordinate 
system  that  moves  with  the  crack  tip.  The  tip  speed  a 
is  assumed  constant. 


* ■ The  summation  convention  is  followed  whereby  repeated 
indices  are  to  be  summed  over  their  range  (i,j =1,2,3). 
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(b)  The  material  is  assumed  to  be  purely 
elastic  with  a strain  energy  function 


FIGURE  IV-48  Contour  Around  a Moving  Crack  Tip 

The  quantities  • and  e.  . are  components  of  the  stress  and 
strain  tensors  anS,  may  be^  related  nonlinearly.  The  material, 
however,  must  deform  in  a reversible  manner. 

Under  these  assumptions,  Sih  (1970)  has  reduced 
Equation  IV-54  to  the  form 

1 ^u.  n 

R(i)  =J  (W+  - pa^  ^ " ^i  ^ * (IV-56) 

As  a direct  consequence  of  assumption  (b) , the  value  of  the 
integral  in  Equation  IV-54  is  independent  of  the  path  of 
integration.  In  other  words,  the  result  obtained  from  any 
two  arbitrary  contours  surrounding  the  moving  crack  tip  will 
be  the  same. 

For  stationary  cracks  (^0),  Equation  IV-56  reduces 

to  the  J-integral  by  Rice  (1^68) . Furthermore,  if  the 

material  is  linear  with  W = — a. .e. .,  it  has  been  shown  by 

2 1]  1] 

Sih  (1970)  that  Equation  IV-56  further  simplifies  to 
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•(lV-5  7) 
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It  should  be  clearly  pointed  out  that  the  contour 
integrals  in  Equations  IV-56  and  IV-57  are  valid  only  if 
assumption  (a)  is  observed.  In  situations  where  the  speed 
changes  or  the  crack  changes  its  direction  of  propagation 
for  the  dynamic  crack  problem  containing  both  stress  inten- 
sity factors  k^  and  k^,  no  path  independent  integrals  have 
been  found;*  refer  to  Eshelby  (1974) . 

(2)  Strain  Energy  Density  Criterion.  A fracture 
criterion  that  is  free  from  the  restrictions  mentioned  in 
subsection  e.(l)  is  based  on  the  strain  energy  density 
factor  S which  has  been  suggested  as  a parameter  for  pre- 
dicting unstable  crack  initiation  and  growth  under  static 
(Sih,  1973)  and  dynamic  (Sih,  1977)  loads  of  a general 
nature,  not  necessarily  symmetrical  with  reference  to  the 
crack  plane.  This  factor  represents  the  product  of  the 
radical  distance,  say  r^ , measured  from  the  crack  tip, 
multiplied  by  the  strain  energy  density  function,  dW/dV, 
commonly  employed  in  the  theory  of  elasticity;  i.e.,  S = 
r. (dW/dV) . The  minimum  value  of  S,  S . , determines  the 

direction  of  crack  initiation  and  the  critical  value  of 

S . , S , determines  incipient  fracture.  The  value  S is 

min  c -C 

assumed  to  be  characteristic  of  the  material  and  taken'as 

the  fracture  toughness.  It  should  be  added  that  experi- 
mental verification  of  this'  criterion  is  currently  very 
limited . 

For  the  nontrivial  problem  of  an  inclined  crack  in 
Figure  IV-49,  the  critical  dynamic  stress  of  cr^  that  the 
plate  could  sustain  can  be  obtained  from  as 


VaF(0,t,e^  ) (lV-58) 

where  the  function  F(/3,t,d  ) is  given  by 

O , 

* The  lower  case  letters  k]_  and  K2  are  used  in  this  sub- 
section for  stress  intensity  factors.  For  instance,  kj 

differs  from  by  a factor  of  VtT” ; i.e.,  K]_  = Vrrk^^. 
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F{/3,t,  00  ) 


(aiif®  sin^/3+2ai2fif2sin)8cos/34-a22f2^os^/3  ) sin^ 

(IV-59) 


The  angle /3  denotes  the  crack  position  with  reference  to  the 
applied  stress  and  the  fracture  angle,  Og , corresponds  to 
the  location  of  S , through  the  coefficients  aj_j  (i,j,=  1,2) 
depending  on  Og an9^?he  elastic  constants  (say,  shear  modulus 
H and.  Poisson's  ratio  i>)  . To  be  emphasized,  is  that  the  same 
S value  is  assumed  for  both  static  and.  dynamic  loads.  The 
time  dependent  functions  f and.  f^  (Sih,  1977)  in  Equation 
IV-59  are  associated  with  the  stress  intensity  factors 


cr 


FIGURE  IV-49  Time  Dependent  Loading  on  an  Inclined  Crack 


k^  = f^(t)aA/~a  sin^^,  k^  = f^  (t)  sin/3cosy8  (IV-60) 


and,  are  unity  for  static  loading.  It  is  important  to  know 
that  under  dynamic  loading,  k^  and  k^  do  not  necessarily 


124 


attain  their  maximum  values  at  the  same  time.  The  way  in 
which  the  combination  of  and  k^  affects  the  critical 
stress,  a^»  can  be  found  from  the  strain  energy  density 
factor : 


S = a^asin^/3[a^^f^sin^/3+2a^2f]_^2®^^^‘^°®^'^^22^2^°®^^  ^ (IV-61) 

First,  the  fracture  angle  in  Equation  IV-58  is  found  from 
the  condition  5S/B0  = 0 such  that  S is  a minimum.  Hence,  a 
plot  of  16juS  . /a^a  versus  c t/a  is  established  as  illus- 
trated in  Figure  IV-50  for  steel  with  v = 0.29.  In  the 


FIGURE  IV-50  Variations  of  Strain  Energy 
Density  Factor  With  Time 

SOURCE:  . (Sih  1977) 
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S- theory,  the  crack  is  assumed  to  extend  in  the  direction 
when  Smin  reaches  the  critical  value  of  a given 
material.  Depending  on  the  angle  /3  and  the  amplitude  of 
the  dynamic  load,  equal  to  Sc  on  any  part  of 

the  curve  in  Figure  IV-50,  or  never  reach  Sc»  Since  the 
objective  of  any  dynamic  analysis  is  to  predict  the  applied 
failure  stress,  it  is  convenient  to  compare  the  ratio  of  the 
dynamic  failure  stress,  Cf,  to  that  of  the  static  failure 
stress,  CTs,  which^may  be  obtained  from  Equation  IV-61  as 
l^^Sc  = 4(1-2  v)CTga.  Keeping  the  fracture  toughness  value, 
S(^,  constant,  the  ratio  Of/oQ  may  be  written  as 


— = V^a-2v) 


N 


(lV-62) 


The  quantity  M stands  for  16^S^^^/a  V^which  may  be  obtained 


from  Figure  IV-50  with  in  Equation  IV-62  corresponding 


to  S 


min 


reaching  S 


A plot  of  CTg  versus  c t/a  is  given  in 


Figure  IV-51.  Each  curve  passes  through  a minimum  in  the 


C2t/a 


FIGURE  IV-51  Dynamic  Load  Amplification 
SOURCE:  (Sih  1977) 
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neighborhood  of  C2t/a  = 2.  For  /S  = 15  ° , this  minimum  is  at 
the  smallest  value  of  C2t/a,  corresponding  to  the  smallest 
value  of  cr^  to  cause  fracture. 

Examples  on  application  of  the  S-Criterion  to 
the  running  crack  problem  can  be  found  in  Sih  (1977-b) . 

(3)  Maximum  Normal  Stress.  The  maximum  normal 
stress  criterion  assumes  that  fracture  occurs  when  the  maxi- 
mum principal  stress  local  to  the  craclc  tip  reaches  a crit- 
ical value.  The  plane  along  which  crack  growth  takes  place 
is  assumed  to  be  normal  to  this  stress.  This  criterion  is 
widely  used  for  quasi-static  problems  and  has  received 
considerable  experimental  verification.  In  a two-dimensional 
situation,  crack  initiates  in  a radial  direction  and  hence 
the  stress  component  under  consideration  is  a circumferential 
direction  given  by 

r 30  _ 0.  ^^2  r . 30  ^ . 0. 

or^  = [cos  — 4-3  cos  -rj  - [sin  — 4-  sin  —J  (lV-63) 

'‘4,^  ' 2 . 2 2 

where  r and  0 are  local  polar  coordinates  shown  in  Figure 
j'y_52 . The  direction  along  which  the  crack  initiates  can 


FIGURE  IV-52  Stress  Components  Near  a 
Stationary  Crack  Tip 
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be  found  from  the  condition  = 0 such  that  a0  is  a 

maximum.  Equation  IV— 63  has  been  applied  to  study  the  mixed 
mode  fracture  under  static  loads  by  Erdogan  and  Sih  (1963) 
and  is  potentially  applicable  to  stationary  cracks  subjected 
to  dynamic  loads. 

The  presence  of  biaxiality  or  triaxiality,  on 
an  element  ahead  of  the  crack,  however,  causes  limitations 
on  the  maximum  principal  stress  criterion.  For  instance, 
the  stress  state  ahead  of  a crack  under  symmetric  loading 
is  hydrostatic  in  that  both  principal  stresses  are  equal  and, 
hence  there  is  no  maximum.  Generally  speaking,  this  cri- 
terion is  effective  only  when  the  magnitude  of  one  of  the 
normal  stress  components  is  much  larger  than  the  others. 

Conceptual  difficulties  also  arise  when 
applying  the  maximum  normal  stress  criterion  to  the  running 
crack  in  Figure  IV-53.  Sih  (1968)  has  shown  that  this  ratio 
is 


(I+Sp  (2s°-H-s|)-4s^s^ 
4s^S2-(1+s=) 


(lV-64) 


is  always  greater  than  unity.  In  Equation  IV-64 


^ 6 = Crack  Velocity 


FIGURE  IV-53  Running  Crack 
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1,2 


(lV-65) 


s . 
D 


[1 


(Vc.) 


j - 


depends  on  the  crack  speed  a and  the  wave  speeds. 


c^  = J{\+2^.)/p,  = »//Li/p 


(lV-66) 


The  constants  X and  p are  the  Lam^  coefficients  and  p is  the 
mass  density.  Equation  IV-64  with  implies  that  crack 

runs  in  a direction  parallel  to  the  maximum  stress  rather 
than  normal  to  it,  a condition  that  contradicts  the  original 
assumption . 

(4)  Maximum  Energy  Release  Rate.  The  amount  of 
energy  released  by  a crack  traveling  in  an  elastic  material 
that  undergoes  no  plastic  deformation  can  be  computed  from 
Equation  IV-56  since  R = G j.  The  subscript  1 means  that 
the  relation  applies  only  for  Mode  I;  i.e.,  symmetrical 
loading.  Depending  on  the  boundary-value  problem,  the 
dynamic  G may  be  related  to  its  corresponding  static  value 
G^g  (Sih,'^1970)  as 


'IS 


1-v 


- (l+s=)=] 


(lV-67) 


for  plane  strain;  see  also  Freund  (1976) . The  quantities 

s . have  already  been  defined  in  Equation  IV-65  and 

F^(s2,S2)  depends  on  the  geometric  and  loading  conditions. 

In  the  limit  as  the  crack  velocity  goes  to  zero,  a-*0,  G 

reduces  to  G,„/  the  static  value: 

Xo 


n(l-v^)k^ 


IS 


(lV-68) 


lim  G 

• ^ JL 

a-O 
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For  cracks  that  turn  or  propagate  without  symmetry, 
the  maximum  energy  release  rate  concept  can  be  used  to 
predict  both  the  condition  for  propagation  and  the  direction 
of  growth;  however,  it  is  difficult  to  apply  in  a rigorous 
fashion.  Analytically,  the  corner  singularity  in  addition 
to  that  at  the  crack  tip  poses  a uniqueness  problem  on  the 
value  of  G in  the  limit  as  the  branched  or  angled  portion 
of  the  main  crack  tends  to  zero.  Also,  the  numerical  eval- 
uation of  G for  a bent  crack  tends  to  smooth  out  stress 
peaks  and  singularities  and  also  presents  difficulties  that 
cannot  be  easily  resolved.  Nevertheless,  work  is  continuing 
on  refining  this  physically  appealing  criterion  (Palinswamy 
and  Knauss,  1977) . 

f . Conclusions 

One  of  the  oldest  ideas  in  fracture  control  is  the  cer- 
tification of  materials  for  service  by  testing  under  high 
rates  of  loading.  As  discussed  above,  impact  tests  such  as 
the  Charpy  test,  the  drop-weight  tear  test,  and  the  dynamic 
tear  test  have  been  (and  continue  to  be)  widely  used.  These 
tests  were  originally  intended  to  give  qualitative  informa- 
tion only;  e.g.,  a temperature  transition.  Recently,  how- 
ever, it  has  become  widely  realized  that  useful  quantitative 
information  can  also  be  extracted  from  the  results  of  these 
tests.  What  is  not  yet  generally  recognized  is  that  the 
conventional  analytical  interpretation  of  such  tests  via  a 
completely  static  point  of  view  can  be  substantially  in 
error.  The  recent  research  work  using  dynamic  (i.e.,  iner- 
tia forces  included)  analyses  described  in  this  report  have 
clearly  shown  this  for  precracked  Charpy  specimens  and  for 
the  drop-weigh  tear  test. 

The  certification  of  materials  for  service  is  presently 
evolving  from  a loose  qualitative  approach  into  a much  more 
quantitative  procedure,  as  illustrated  by  the  Air  Force's 
fracture  tolerance  guidelines  adopted  for  the  B-1  bomber. 
Other  Department  of  Defense  concerns  with  high  loading  rate 
applications  are  cited  elsewhere  in  this  report.  Hence, 
there  is  a strong  need  for  dynamic  analyses  to  accompany 
high  strain  rate  testing  in  order  that  material  properties 
appropriate  for  engineering  structural  design  can  be  cor- 
rectly extracted  from  such  tests. 

The  work  reviewed  in  this  section  reveals  that  there 
are  two  main  approaches  in  dynamic  fracture  analysis  and 
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testing.  Some  investigators  have  considered  the  dynamic 
propagation  of  a crack  initiated  under  quasi-static  loading 
conditions.  Other  investigators  have  confined  themselves 
to  considering  the  conditions  prior  to  the  initiation  of 
growth  under  dynamic  loading.  Obviously,  both  aspects  are 
involved  in  a proper  analysis  of  the  problem.  In  our  judg- 
ment, the  key  contribution  that  needs  to  be  made  in  this 
field  is  the  development  of  an  analysis  capability  that 
considers  both  dynamic  crack  initiation  and  propagation. 

Only  with  such  a tool  can  meaningful,  dynamic  fracture  data 
be  extracted  from  high  loading  rate  experiments. 

It  is  not  difficult  to  conclude  that  the  existing 
theories  are  not  sufficient  for  a general  treatment  of 
material  response  to  high  strain  rate  loading,  particularly 
in  handling  the  effects  of  specimen  size  and  the  range  of 
operating  temperatures.  From  the  fracture  mechanics  point 
of  view,  many  of  the  proposed  fracture  toughness  parameters 
are  still  sensitive  to  the  operating  conditions,  and  undergo 
sharp  variations,  normally  referred  to  as  the  brittle-ductile 
transition  behavior.  Figure  IV-54.  In  principle,  a material 
parameter  should  remain  invariant  with  respect  to  changes  in 
specimen  size,  and  it  is  very  desirable  if  this  invariance 
exists  for  temperature  and  loading  rate  changes  as  well. 


Transition 


Size 


Transition  Transition 


Temperature  , Loading  Rate 


FIGURE  IV-54  Transition  Behavior 
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If  a material  exhibits  a distinct  point  of  instability 
on  the  load-deformation  curve  regardless  of  whether  it  is 
linear  (Figure  IV-55)  or  nonlinear  (Figure  IV-56) , the  frac- 
ture toughness  value  (if  interpreted  as  (AW/Aa)c  insta- 
bility, where  AW^  is  the  local  energy  needed  to  cause  insta- 
bility, should  be  the  same  for  both  cases  provided  that  the 
material  behavior  within  which  Aa  is  created,  remains  un- 
changed,. It  is  believed  important  to  adopt  a local  view  of 
the  fracture  process  in  the  failure  analysis  of  structural 
components  where  predictions  on  allowable  load  and/or  net 
section  size  must  be  made;  such  a view  is  used  in  Section 
B.5.C.  in  which  three-dimensional  crack  growth  predictions 
are  discussed.  As  noted  previously,  the  critical  value  of 
J may  be  interpreted  as  a local  fracture  initiaticn  cri- 
terion (cf . Equation  IV-24)  , and.  therefore  the  great 
success  achieved  with  this  criterion  supports  this  point  of 
view.  However,  ductile  fracture  for  static  and  dynamic  con- 
ditions is  a subject  which  is  still  little  understood,  and. 
requires  extensive  future  research, 

I 

Elastic 


FIGURE  IV-55  Linear  (Brittle)  Behavior 
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Elastic-Plastic 


FIGURE  IV— 56  Nonlinear  (Ductile)  Behavior 


5 • Three-Dimensional  Elastic-Plastic  Fracture 
Considerations . 


a . Introduction 

With  the  advent  of  fracture  mechanics,  it  is  now 
possible  to  design  structures  using  the  concept  of  fracture 
toughness  in  addition  to  the  traditional  practice  of  using 
strength  and  structural  buckling  instability  criteria.  One 
of  the  major  shortcomings  of  the  strength  and  toughness 
approach  lies  in  the  inability  of  the  analyst  to  translate 
small  specimen  laboratory  data  to  the  design  of  large  size 
structures  made  of  tough  materials.  This  lack  of  under- 
standing of  fracture  size  effect  cannot  be  resolved  by 
testing  alone. 
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A realistic  modeling  of  ductile  fracture*  is  a pre- 
requisite for  resolving  the  effect  of  size  on  the  brittle- 
ductile  transition.**  The  basic  features  of  ductile  fracture 
can  be  observed  from  macrofractography  and  involve  three 
distinct  stages  (Sih,  1976)  which  should  be  included  in  an 
elastic-plastic  crack  analysis.  Fracture  first  originates 
in  a region  subjected  to  nearly  hydrostatic  tension.  It 
then  changes  from  slow  or  stable  growth  to  rapid  or  unstable 
crack  propagation.  Final  separation  of  the  material  occurs 
when  the  crack  turns  away  from  the  normal  to  the  free  sur- 
face of  the  solid.  This  slanted  fracture  surface  is 
referred  to  usually  as  the  shear  lip  or  cup-and-cone . The 
transition  from  slow  to  rapid  fracture  occurs  so  quickly 
that  it  is  difficult  to  observe  the  load-deformation  curve. 
The  last  ligament  of  material  separation,  however,  does 
coincide  with  a sharp  drop  in  load,.  Because  the  fracturing 
process  is  seldom  symmetrical,  the  shear  lip  or  cup-and- 
cone  appears  only  on  one  half  of  the  broken  specimen. 

b.  Ductile  Fracture 

Preliminary  analysis  of  a tensile  bar  with  a crack 
shows  that  plastic  zones  are  formed  on  both  sides  of  the 
macrocrack  as  it  moves  through  the  material.  Figure  IV-57a. 
Initial  instability  corresponds  to  essentially  brittle 
fracture  where  the  fracture  path  is  almost  flat.  If  the 
crack  is  allowed  to  approach  the  boundary,  the  bar  begins 


* Elastic-plastic  stress  analyses  of  two-dimensional  crack 
models  are  found  to  be  inadequate  as  they  do  not  account 
for  necking  and  the  slanted  fracture  terminating  at  the  free 
surface.  This  is  a basic  feature  of  ductile  fracture  and. 
is  a three-dimensional  phenomenon. 

**  It  should  be  noted,  however,  that  data  obtained  for 
various  structural  steels  demonstrate  that  a true  plane- 
strain  fracture  toughness  transition  exists  which  is  inde- 
pendent of  specimen  thickness.  Fractographic  analyses  show 
that  this  transition  temperature  is  associated  with  the  on- 
set of  change  in  the  microscopic  fracture  initiation  mode 
at  the  crack  tip;  the  mode  changes  continuously  from  cleav- 
age at  the  low  temperature  end  to  ductile  tear  at  the  upper 
end  (Rolfe  and.  Barsom,  1977)  . 
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to  neck  until  the  plastic  zones  intersect  the  free  boundary. 
An  additional  island  of  yielded  material  has  been  found. 

(Sih,  1977-a) , leaving  a forked  region  of  elastic  material. 
Should  the  bar  fail  unstably  at  this  point,  a possible  path 
of  fracture  is  through  the  forked  elastic  region  as  shown 
in  Figure  IV-57b.  For  the  more  ductile  material,  the  bar 
may  neck  further  until  the  entire  ligament  between  the  crack 
and  outer  specimen  surface  is  yielded.  Figure  IV-57c. 

Further  loading  suggests  a further  distinction  of 
zones  of  distortion  from  dilatation.  This  is  illustrated  in 
Figure  IV— 57d  where  the  zones  of  distortion  are  marked  as 
second  degree  yielding.  A possible  path  of  fracture  is 
again  suggested  to  take  place  on  a slant  in  which  the 
material  is  distorted  to  a lesser  degree  than  that  in  the 
cross-hatched  region. 

To  reiterate.  Figures  IV-57a  through  d illustrate 
the  various  degrees  of  yielding  that  correspond  to  materials 
with  different  degrees  of  ductility  and  different  loads  of 
instability  on  the  load  deformation  curve.  The  mechanics 
problem  of  ductile  fracture  is  to  predict  the  fracture  path 
and  load  in  terms  of  the  continuum  material  parameters  such 
as  yield  strength,  strain  hardening  exponent,  etc.,  as  well 
as  the  fracture,  toughness  associated  with  the  local  behavior 
of  the  material  at  the  crack  tip. 

Similar  phenomena  are  observed  in  the  fracture  of 
moderately  thick  plate  specimens  with  a through  crack  such 
that  the  combinations  of  material  properties  and  plate 
geometry  would  not  satify  the  ASTM  requirements  (ASTM 
Special  Technical  Publication  No.  410,  1966)  for  plane 
strain  where  the  smallest  geometric  length  parameter  must 
be  greater  than  or  equal  to  2.5  (K^^/ay)^.  Here,  as  before, 

Kjc  is  the  valid  critical  Mode  I stress  intensity  factor  and 
is  the  uniaxial  yield  strength  of  the  material.  Figure 

IV-58  shows  that  the  critical  stress  intensity  factor  K^* 
becomes  or  geometry  independent  only  when  the  plate  is 

sufficiently  thick.  The  fracture  surface  appearance  as 


*Strictly  speaking  should  not  be  referred  to  as  the  plane 
stress  fracture  toughness  value  because  the  mode  of  fracture 
for  this  specimen  is  in  no  way  related  to  the  plane  stress 
plasticity  crack  solution. 
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a.  Relatively 

Brittle  Fiehavior 


b.  Relatively 

Ductile  Behavior 
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d.  Very  Ductile 
Beluivio!’ 


FIGURE  IV-57  Brittle  and  Ductile  Fracture 
Models 
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PLATE  THICKNESS 


FIGURE  IV-58  Variation  of  K with  Plate  Thickness 

c 


shown  in  Figure  IV-59a  is  almost  all  flat;  i.e.,  the  through 
crack  grows  in  its  initial  plane  normal  to  the  applied  load. 
As  the  plate  thickness  is  decreased,  the  crack  will  initially 
tunnel  and  then  deviate  from  its  own  plane  near  the  plate 
surfaces  forming  "shear  lips."  The  crack  growth  is  stable 
at  first,  advancing  from  its  initial  configuration  in  a 
planar  fashion  marked  by  the  dotted  area  in  Figure  IV-59b 
with  a curved  front.  Upon  reaching  the  point  of  instability, 
the  curved  crack  breaks  through  the  plate.  This  can  be 
traced  on  the  loading  history  curve.  Again,  local  conditions 
will  determine  whether  the  shear  lip  will  form  on  one  crack 
surface  or  the  other.  Further  reduction  in  plate  thickness 
leads  to  the  decrease  of  flat  surface.  Figure  IV-59c,  and  a 
slanted  fracture  pattern  is  finally  observed. 

c . Three-Dimensional  Slow  Crack  Growth  Predictions 

In  order  to  predict  the  aforementioned  ductile 
fracture  phenomenon,  it  is  necessary  to  perform  a three- 
dimensional  elastic-plastic  stress  analysis  and  to  have  a 
fracture  criterion  that  is  amenable  to  analyzing  crack 
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a.  Almost  All 
Flat  Fracture 


FIGURE  IV-59  Appearance  of  Fracture  Surfaces 

tunneling  and  mixed  mode  crack  extension  of  the  type 
described  earlier.  It  is  desirable  to  analyze  the  stable 
growth  process  of  a tunne 1 crack  in  a plate  of  ductile 
material.  The  influence  of  crack  front  curvature  can  be 
examined  with  the.  aid  of  a three-dimensional  elastic  finite 
element  analysis  in  conjunction  with  a local  growth  criterion 
to  predict  increments  of  crack  growth  and  corresponding 
changes  in  crack  front  shape.  In  general,  both  the  direction 
and  magnitude  of  local  crack  growth  from  the  current  crack 
front  are  expected  to  vary  with  position  along  that  front. 

As  a possible  crack  growth  criterion,  the  strain-energy- 
density  theory  may  be  used  (Sih,  1973) . This  criterion 
assumes  that  the  path  of  growth  from  each  point  along  the 
crack  edge  will  follow  the  minimum  strain  energy  density 
path  emanating  from  that  point.  Further,  growth  at  points 
along  the  current  crack  front  will  initiate  when  the  strain 
energy  density  at  a "core"  distance  r along  the  minimum 
path  reaches  a prescribed  value  (dW/d^) . 
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Hence,  growth  of  the  crack  front  is  approximated 
by  discrete  increments  of  growth.  The  amount  of  growth  at 
a point  along  the  crack  front  in  an  increment  is  taken  to 
be  the  distance  along  the  minimum  strain  energy  density 
path,  (dW/dV^^  , to  the  point  where  (dW/dV^.^  reaches  a 
preset  value, ^^dW/dV^ , which  depends  on  the^material . (This 
value  of  (dW/dV^,  assumed  constant,  is  to  be  determined 
experimentally  for  a particular  material.)  Figure  IV-60 
gives  a descriptive  plot  of  strain  energy  density 
against  the  radial  distance  r in  a plane  normal  to  the 

applied  load  from  the  associated  points  A,B, ,F  on  the 

crack  front.  The  new  crack  profile  is  determined  from  the 
various  values  of  r and  z/h  on  the  graph  where  intersected 
by  (dW/dV^.  The  selected  value  of  (dW/dV)  will  obviously 
affect  the  magnitude  of  growth  and  shape  of  the  subsequent 
crack  fronts.  In  Figure  IV-60,  z is  the  thickness  coordin- 
ate measured  from  the  mid-plane  such  that  z/h  = 0.5  locates 
the  plate  surface.  This  process  can  be  repeated  and  the 
crack  front  geometry  with  each  of  growth  can  be 
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FIGURE  IV-60 


Crack  Growth  Predictions  Based  on 
Strain  Energy  Density 
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predicted  as  displayed  in  Figure  IV-61.  Crack  tunneling 
will  be  a stable  process  as  long  as  there  is  sufficient 
constraint  from  the  necking  of  the  plate  surfaces.  As  this 
process  continues,  there  will  be  a point  on  the  load- 
deformation  diagram  at  which  the  crack  begins  to  grow 
rapidly;  that  is,  when  the  specimen  loses  its  structural 
integrity.  It  is  expected  that  once  the  phenomenon  of 
ductile  fracture  is  analyzed  and  the  shape  and  size  of 
these  tunneling  cracks  are  calculated,  for  a variety  of 
situations,  simple  engineering  formulas  for  design  applica- 
tions can  be  established. 

6 . Methods  of  Dynamic  Failure  Analysis 

a . Introduction 

For  the  types  of  loading  described  in  Section  III, 
dynamic  (rather  than  static)  stress  analyses  are  normally 
required  in  assessing  structural  integrity.  Similarly,  as 
emphasized  earlier  in  Section  IV, B,  dynamic  analysis  of  test 


FIGURE  IV  -61  Predicted  Increments  of 
Crack  Growth 


140 


specimens  must  be  made  in  most  cases  if  basic  material  and 
fracture  parameters  are  to  be  properly  extracted  from  test 
data.  In  this  section,  a brief  review  will  be  made  of 
dynamic  analysis  techniques.  The  area  of  numerical  methods 
of  fracture  analysis  will  be  emphasized  because  significant 
advances  have  been  made  very  recently,  and  it  is  believed 
these  methods,  when  combined  with  general  structural  analysis 
computer  codes,  will  provide  the  engineer  with  powerful, 
cost-effective  tools  that  will  aid  greatly  the  design  and 
evaluation  of  complex  structures. 

Dynamic  buckling  represents  another  mode  of  failure; 
it  may  lead  to  the  total  collapse  of  a structure  and/or 
fracture,  or  at  least  excessive  deformations.  Furthermore, 
it  is  often  true  that  dynamic  buckling  loads  are  less  than 
the  corresponding  quasi-static  critical  loads,  depending  on 
the  loading  duration.  In  view  of  the  importance  of  the 
problem  of  dynamic  buckling  of  shells,  with  or  without 
stiffeners,  this  area  will  also  be  discussed.  Buckling  is 
covered  first,  and  it  is  followed  by  a more  extensive  review 
of  dynamic  fracture. 

r 

b . Dynamic  Deformation  and  Buckling  of  Shells 
and  Rings  . 

Numerical  methods  exist  that  are  suitable  for  the 
analysis  of  a great  many  dynamic  structural  shell  problems, 
with  or  without  buckling  action;  e.g.,  see  Pilkey  et  al. 
(1974),  Almroth  et  al.  (1976),  Stricklin  et  al.  (1974),  and 
Yeung  and  Welch,  (1977) . There  are,  however,  certain  types 
of  problems  for  which  existing  techniques  are  inadequate. 

In  particular,  very  little  progress  has  been  made  with 
respect  to  (1)  the  prediction  of  large  deflections  and 
rotations  (as  opposed  to  small  or  moderate  values) , (2) 

interactions  with  a surrounding  medium,  and  (3)  inclusion 
of  realistic  inelastic  constitutive  equations  (especially 
if  unloading  and  reloading  effects  are  needed  in  the 
prediction  of  failure) . Furthermore,  the  problem  of 
asymetric  dynamic  buckling  of  shells  is  so  complicated  that 
the  criteria  for  buckling  are  not  even  reasonably  well  under- 
stood; Akkas,  1976,  however,  has  made  some  recent  progress 
in  this  area  for  spherical  shells. 

It  should  be  pointed  out  that  there  has  not  always 
been  complete  agreement  on  a useful  definition  of  dynamic 
buckling,  whether  symmetric  or  asymmetric.  Currently,  the 
most  generally  accepted  criterion  is  based  on  that  proposed 
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by  Budiansky  and  Roth  (1962) : the  peak  in  time  history  of 
a deflection  is  predicted  as  a function  of  the  amplitude  of 
the  applied  pressure  or  load;  the  amplitude  for  which  this 
peak  value  undergoes  a rapid  increase  is  the  buckling  load 
or  pressure.  Unfortunately,  to  establish  a dynamic  buckling 
load  according  to  this  criterion  requires  a greater  effort 
than  that  involved  in  the  determination  of  all  details  of 
the  response  to  a given  loading  history;  the  critical  ampli- 
tude can  only  be  determined  after  a number  of  transient 
response  analyses  have  been  performed.  Nevertheless,  the 
concept  of  a dynamic  buckling  load  or  pressure  has  the 
advantage  of  making  possible  the  construction  of  design 
charts  for  at  least  relatively  simple  structures. 

If  one  simplifies  the  material  response  to  a rigid- 
plastic  material  with  linear  strain  hardening,  it  is  often 
possible  to  predict  analytically  the  dynamic  buckling  of 
rings  and  shells  (Jones  and  Okawa,  1976) . Incorporation  of 
rate-dependence  of  material  behavior  through  the  simple 
method  of  using  a yield  stress  evaluated  at  the  initial 
strain  rate  is  found  to  be  quite  accurate  in  some  problems 
(Perrone,  1965)  . Agreement  between  analytical  predictions 
and  experimental  results  for  average  permanent  radial 
deformation  of  rings  due  to  impulsive  external  pressure 
loading  is  quite  good;  similarly,  experimental  verification 
exists  for  the  theoretical  threshold  value  of  impulse  which 
produces  significant  wrinkles  (Jones  and  Okawa,  1976)  . 

c • A Review  of  Some  Recent  Literature  on  Dynamic 
Fracture  Analysis. 

Much  of  the  work  in  the  area  of  dynamic  fracture 
analysis  has  been  done  on  simple  geometries,  often  two- 
dimensional.  The  use  of  such  analytical  studies  is  three- 
fold; (1)  analytical  solution  using  simple  shapes  is 
relatively  expedient  (indeed,  often  at  the  time,  treatment 
of  more  complex  shapes  is  not  feasible) ; (2)  comparison  of 

the  dynamic  analysis  results  with  static  results  for  some 
geometries  provides  a simple  means  of  understanding  the 
dynamic  effects;  and  finally,  (3)  results  of  the  analysis 
can  be  checked  against  experimental  results,  thus  not  only 
providing  a means  of  validating  the  analysis  procedure  but 
also  helping  one  to  understand  the  experiment  more  thoroughly. 

The  methods  of  dynamic  fracture  anlaysis  can  be 
broadly  classified  into  three  categories: 


142 


Classical,  finite  difference,  finite  element 


Classical  Methods  - Classical  Methods  (i.e.,  analytical 
methods  of  solution)  have  been  in  vogue  for  the  longest 
period  of  time.  These  methods  have  been  successfully 
applied  to  many  simple  geometries.  For  example,  Sih  and 
Embley  (1972)  applied  the  integral  transforms  method 
coupled  with  the  technique  of  Cagniard  to  the  solution  of 
impact  response  of  a finite  crack.  The  Laplace  transform 
technique  was  applied,  by  Freund,  and  Rice  (1974)  to  the 
elastodynamic  crack  tip  stress  field..  More  recently, 
influence  of  specimen  boundaries  on  the  dynamic  stress 
intensity  factor  was  studied  by  Chen  and  Sih  (1974)  by 
means  of  Laplace  and  Fourier  transforms. 

Problems  of  running  cracks  have  also  been  analyzed, 
by  the  classical  methods.  Rapid  extension  of  a penny-shaped 
crack  under  torsion  was  studied  by  Kennedy  and  Achenbach 
(1972),  using  the  Green's  function  technique.  Achenbach 
(1975)  analytically  studied,  the  elastodynamic  bifurcation 
of  running  crack  in  antiplane  strain.  Chen  and  Sih  (1975) 
studied  the  problem  of  scattering  of  plane  waves  by  a propa- 
gating crack  using  the  Fourier  transform.  The  effect  of 
crack  speed  and  other  related  phenomena  in  the  dynamic 
ductile  fracture  of  a central  crack  were  analyzed  by  Tsai 
(1976)  with  the  help  of  the  integral  transform  method. 

It  must  be  stated  here  that  most  of  the  applications 
above  are  phenomena  oriented.  To  study  a problem  or  gain 
understanding  of  certain  phenomena,  the  researcher  develops 
his  field  equations  and  solves  those  equations  by  a method, 
most  attractive  to  him  for  various  reasons.  That  the  results 
obtained  from  various  such  stiidies  are  useful  cannot  be 
disputed.  On  the  other  hand,  these  methods  cannot  solve  a 
general  geometric  shape  and  are  otherwise  limited  in  scope. 
Therefore,  whereas  individual  researchers  will  continue  to 
use  these  methods,  at  the  present  time  the  classical  methods 
do  not  qualify  as  a general  purpose  engineering  analysis 
tool. 

Finite  Difference  Method  - The  finite  difference  method,  is 
a powerful  numerical  tool  for  solving  a complex  boundary 
value  problem.  The  method  itself  is  widely  illustrated,  in 
the  literature  and.  therefore  will  not  be  discussed  here. 

We  present  instead  some  of  the  recent  applications  of  the 
method  to  the  fracture  problems. 
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A two-dimensional  dynamic  fracture  experiment  was 
simulated  using  the  finite  difference  method  by  Stockl  and 
Auer  (1976) . They  solved  two  problems.  In  the  first  prob- 
lem, they  studied  the  constant  velocity  propagation  of  a 
brittle  tensile  crack  in  finite  plane-strain  model  under 
uniaxial  tension.  In  the  second  problem,  the  model  was 
modified  to  allow  for  a specified  variable  crack  propagation 
velocity  in  a nonhomogeneous  prestress  field.  The  authors 
advocate  the  use  of  a transition  zone  at  the  crack  tip 
between  the  stress-free  surfaces  and  the  uncracked  state. 
They  point  out  that  if  a transition  zone  is  not  used,  the 
"cracking  node"  is  released  suddenly  and  the  energy  stored 
at  the  crack  tip  is  converted  into  large  oscillations.  In 
solving  the  first  problem  through  250  time  steps  with  a 
110  X 110  grid,  it  took  3.5  min  on  UNIVAC  1108  computer. 

Shmuely  and  Peretz  (1976)  modeled  the  static  and 
dynamic  DCB  tests  using  the  finite  difference  method.  The 
authors  use  the  stress  normal  to  the  crack  plane  at  one- 
half  the  finite  difference  grid  length  ahead  of  the  crack 
tip  as  the  criterion  for  predicting  speed.  When  this  stress 
reaches  a specified  constant  value,  the  crack  length  is 
extended  one  grid  length.  They  found  that  for  any  one  value 
0^  , the  crack  propagation  speed  remains  approximately 
constant.  However,  it  should  be  noted  that  the  criterion 
employed  is  purely  an  assumed  one,  and  therefore  is  not 
necessarily  the  proper  one. 

HEMP  and  HEMP  3D  (Chen  and  Wilkins,  Aug  1976  and 
Nov  1976)  are  the  well-known  finite  difference  programs. 
These  programs  were  originally  designed  and  developed  for 
solving  dynamic  problems  in  continuum  mechanics  - and  have 
been  extended  to  analyze  the  stress  field  around  cracks  in 
a two  or  three-dimensional  solid  with  finite  geometry 
subjected  to  dynamic  loads.  These  codes  employ  an  explicit 
time  integration  of  the  Lagrangian  formulation  of  the 
equations  of  continuum  mechanics  in  which  the  calculation 
grid  moves  with  the  material.  A special  extrapolation 
scheme  to  simulate  the  singularity  near  the  crack  tip  has 
been  introduced  in  the  program.  There  are  two  methods 
employed  for  modeling  the  separation  of  surfaces  during 
crack  propagation.  In  one  method  the  grid  is  actually  split 
with  the  help  of  sliding  interface  logic;  in  the  other  case 
the  effects  of  the  free  surface  boundary  conditions  are 
introduced  without  actually  splitting  the  grid.  A smoothing 
procedure  is  used  to  avoid  spurious  oscillations  during  a 
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sudden  split  in  the  grid.  The  relevant  constitutive  proper- 
ties can  be  rate-9ependent . The  material  model  can  be 
elastic,  plastic,  etc.  A discussion  on  crack  propagation 
criteria  is  givfen  in  Chen  and  Wilkins  (1977). 

Figure.  IV-62  shows  the  Lagrange  grid,  used  in  cal- 
culations of  thb  effect  of  a .semielliptical  surface  crack 
and  an  imbedded!  (internal)  elliptical  crack  for  a specimen 
which  is  subjected  to  a constant,  suddenly  applied  symmet- 
rical tension  Ibading.  Figures  IV-63  and  lV-64  illustrate 


Free  surface  for  sur- 
face crack  and  plane 
of  symmetry  for  im- 
bedded crack 
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FIGURE  IV-6,:?  Lagrange  Grid  Used  in  Calculations 
of  a Sjemielliptic  Surface  Crack  and 
..an^  Injbedded  Elliptical  Crack. 


SOURCE:  Chen  and  Wilkins.,,  November  1976 
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SOURCE: 


SOURCE: 


FIGURE  IV-63  Normalized  Stress  Intensity 

Factor  versus  Time  for  Imbedded 
Elliptical  Crack 

Chen  and  Wilkins,  November  1976 


FIGURE  lV-64  Normalized  Stress  Intensity 

Factor  versus  Time  for  Semielliptic 
Surface  Crack 

Chen  amd  Wilkins,  November  1976 
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the  variation  in  the  stress  intensity  factor  with  time  for 
these  two  cases.  The  geometry  of  axial  and  circumferential 
cracks  in  the  wall  of  a cylindrical  pipe  is  shown  in  Figure 
IV-65;  stress  intensity  factors  due  to  a constant,  suddenly 
applied  pressure  for  these  two  types  of  cracks  are  shown  in 
Figures  IV-66  and  IV-67  for  several  time  steps.  More 
recently,  Wilkins  .(1977)  extended  the  HEMP  codes  to  running 
cracks  using  an  accumulated  damage  model.  In  spite  of  the 
power  of  these  codes,  severe  limitations  to  our  prediction 
capability  still  exist  because  of  uncertainties  in  the 
material  characterization  and  the  crack  initiation  and  growth 
criteria. 


FIGURE  IV-65  Geometry  of  the  Axia]  and  Circumferential 
Cracks  in  the  Wall  of  a Section  of 
Cylindrical  Pipe 

SOURCE:  Chen  and  Wilkins,  November  1976 
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3k ^/4P 


FIGURE  IV-66  Normalized  Stress  Intensity  Factor  3k^/4P 
versus  R at  Several  Steps  for  the  Axial 
Crack 

SOURCE:  Chen  and  Wilkins,  November  1976 
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FIGURE  IV-67  Normalized  Stress  Intensity  Factor  3k^/4P 
versus  R at  Several  Time  Steps  for  the 
Circumferential  Crack 

SOURCE:  Chen  and  Wilkins,  November  1976 


Finite  Element  Method  - Development  of  the  finite  difference 
method  of  solution  preceded  the  development  of  the  finite 
element  method  (FEM) . However,  the  FEM  is  used  in  most 
stress  analysis  applications  today.  Whereas  both  methods  are 
general  purpose  and  applicable  to  a wide  spectrum  of  prob- 
lems, the  finite  element  method  is  much  more  versatile  and 
efficient  in  terms  of  adapting  to  almost  any  shape  of  the 
medium  and  in  applying  the  boundary  conditions.  In  this 
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section  we  shall  review  some  of  the  recent  applications  of 
the  FEM  in  dynamic  fracture  mechanics.  The  actual  develop- 
ment of  the  method  as  applied  to  fracture  problems  will  be 
discussed  in  the  next  section. 

Aberson  et  al.  (1976)  applied  the  FEM  to  dynamic 
fracture  problems.  They  used  a computer  program  for  two- 
dimensional  elastic  structures.  It  consists  of  rectangular 
elements,  one  for  opening  mode  problems  and  one  for  mixed 
mode  problems.  The  displacements  within  the  elements  are 
represented  by  a finite  series  of  Williams'  eigen  functions. 
Consistent,  lumped  mass  matrices  are  included  in  the  program. 
Figure  IV-68  shows  the  finite  element  model  of  a plane-strain 
rectangular  strip  which  is  centrally  cracked.  The  variation 
in  stress  intensity  factor  with  time  under  a suddenly  applied 
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FIGURE  IV-68  The  Finite  Element  Model 
SOURCE;  Aberson  et  al,  1976 


150 


tensioi  is  shown  in  Figure  IV-69.  Figure  IV-70  is  a finite 
element  model  of  a dynamic  tear  test  specimen  of  Homalite- 
100,  a photoelastic  material.  (The  specimen  support  and 
loading  configuration  are  essentially  those  of  a Charpy  test 
except  that  a translating  steel  hammer  is  used  in  place  of 
a pendulum  for  the  impact  source.)  The  node  corresponding  to 
the  impact  point  is  assigned  a mass  equal  to  one-half  of  the 
mass  M of  the  hammer.  This  node  is  then  given  an  initial 
velocity  Vq  equal  to  the  impact  velocity;  all  other  nodes 
have  zero  initial  velocity.  Since  the  specimen  is  thin, 
plane-stress  elements  were  used;  however,  it  should  be 
recognized  that  plane  stress  close  to  the  crack  tip  does  not 


TIME,  /i‘SEC. 

FIGURE  IV-69  Stress- Intens ity  Factors  Predicted  Using 
Different  Mass  Distributions 

SOURCE:  Aberson  et  al,  1976 
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M • 1,47  Kg 


FIGURE  IV-70  Finite  Element  Model  of  a Dynamic-Tear 
Test  Specimen 

SOURCE:  Aberson  et  al,  1976 


exist.  The  numerical  results  are  given  in  Figure  IV-71  for 
a crack  of  constant  length  of  0.005  in. 

Aberson  et  al.  (1976)  have  also  reported  a parallel 
study  on  the  running  cracks  using  the  finite  element  method. 
In  this  case  they  used  simple  constant  strain  triangles. 

The  finite  element  mesh  is  particularly  refined  in  the  vicin- 
ity of  the  crack  path.  A constant  crack  propagation  speed 
is  imposed.  As  the  crack  reaches  a node  its  constraint  is 
released.  The  results  are  preliminary  in  nature  and  need 
further  refinement. 

Ayres  (1976)  has  performed  a two-dimensional 
elastic-plastic  finite  element  analysis  of  the  precracked 
Charpy  V notch  specimen  with  a crack  of  constant  length 
using  the  MARC  (1976)  finite  element  program.  Figure  IV-72 
shows  a finite  element  model  for  this  specimen.  The  mesh 
at  the  crack  tip  is  constructed  with  8 node  isoparametric 
quadrilaterals  with  midside  nodes  one  quarter  of  the  way 
from  the  crack  tip  to  the  element  corner.  The  nodes  of 
each  element  at  the  crack  tip  are  distinct  and  are  free  to 
move  independently  in  the  direction  of  the  crack.  This  type 
of  modeling  provides  a 1/r  strain  singularity  appropriate 
for  plastic  behavior  (cf.  Section  d) . 
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FIGURE  lV-71  Stress  Intensity  Factor 
and  Hammer  Acceleration 
for  a Dynamic-Tear- Test 


SOURCE:  Aberson  et  al,  1976 


Figure  IV-73  illustrates  the  comparison  between 
the  computed  and  measured  applied  load-time  curve  for  the 
case  when  the  specimen  is  loaded  at  a constant  velocity. 
Although  not  perfect,"  comparison  seems  to  show  reasonable 
agreement.  In  the  same  analytical  study,  Ayres  also  showed 
that  the  dynamic  crack  opening  displacement  near  the  crack 
tip  (.Ola)  is  much  greater  than  the  statically  computed 
value . 

Applicability  of  the  MARC  type  finite  element 
analysis  is  apparent  from  Ayres'  numerical  studies.  In 
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FIGURE  IV-72  Finite  Element  Model  of  Precracked 
CharpY  Specimen 

SOURCE:  Ayres,  1976 


FIGURE  IV-73  Load-Time  Curve  Computed  from  Elastic- 
Dynamic  Analysis  with  Constant  Velocity 
of  Loaded  Point 


SOURCE:  Ayres,  1976 
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principle,  similar  three-dimensional  elasto-plastic , dynamic 
analyses  can  be  performed  using  the  same  program.  Consti- 
tutive properties  for  elasto-plastic  behavior  can  be  input; 
rate  dependence  can  be  incorporated  only  through  a specified 
constant  strain  rate  . Furthermore,  a crack  propagation 
criterion  is  not  used  in  MARC  (1976) . 

d . Discussion  of  the  Finite  Element  Method.. 

As  indicated  earlier  in  this  section,  the  FEM  is  the 
most  versatile  tool  of  structural  analysis  that  is  available 
today.  As  the  field,  of  analytical  fracture  mechanics  grows 
and  the  practical  application  of  analysis  to  important  design 
problems  spreads,  greater  development  and  use  of  this  method 
is  inevitable.  The  basics  of  the  finite  element  method  are 
amply  described  in  textbooks  (e.g.,  Zienkiewicz,  1971). 
However,  for  reference  purposes,  it  will  be  desirable  to 
describe  them  very  briefly.  Important  aspects  of  the  crack 
tip  finite  element  will  be  reviewed  in  more  detail. 

For  linear  media,  the  equation  of  motion  to  be  solved 
is 

MU  + CU  + KU  = F(t)  (lV-69) 

where 

M = mass  matrix 
U = displacement  vector 
C = damping  matrix 
K = stiffness  matrix 
F = force  vector 
(‘)=  time  derivative 

Let  us  consider  a continuum  finite  element  (2  or  3 dimen- 
sional) with  M degrees  of  freedom.  The  displacement  at  any 
point  withih  the  element  is  given  by 

.-rill  m 

Ui  = N u^  (lV-70) 


where 

m 

N = shape  function  (or  interpolation  function) 

at  node  m 
m , ^ , 

u.  = value  of  displacement  u.  at  node  m 
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The  repeated  sub-  or  superscripts  imply  summation.  The  sub 
script  'i'  is  for  spatial  directions  in  the  Cartesian  space 
In  2D,  i = 1,2;  whereas  in  3D,  i = 1,2,3.  The  strain  is 
given  by 


where 


„m  m 
«r . . = Z . , u 
ij  ijk  k 


(lV-71) 


.m 


1 3 


ijk 


(lV-72) 


in  which  6^^  represents  a Kronecker  delta. 


The  linear  elastic  stress-strain  relationship  is 


'^ij  ^ijkj^^kL 


(IV-73) 


The  nodal  forces  can  be  obtained  from  the  principle  of 
virtual  work  for  an  element;  there  results 


(IV-74) 


where  the  integral  is  over  the  volume  of  the  element. 
Equations  IV-69  to  IV-74  yield,  an  element  nodal  force- 
displacement  relationship. 


_m  , mn  n 
P . = K. . u . 
1 13  3 


(IV-75) 


in  which  is  the  stiffness  matrix, 

13 


Kinn 

13 


-iti 


5N  BN' 


n 


Bx,  Bx  ^ 
k Z 
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(lV-76) 
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The  so-called  consistent  mass  matrix  for  the  element  is 
obtained  from 


mn  i*  ^ m n ^ 

M = j pN  N dV 


(IV-77) 


The  lumped  mass  matrix  is  derived  by  lumping  the  off-diagonal 
terms;  viz.. 


m 


M 


N dV 


{lV-78) 


The  overall  structural  mass  and  stiffness  matrices  are 
obtained  by  appropriately  combining  the  element  mass  and. 
stiffness  matrices.  There  are  various  methods  of  specify- 
ing the  damping  matrix,  which  are  described  in  the  litera- 
ture. The  following  observations  are  made  on  the  above 
discussion: 

(1)  If  all  the  finite  elements,  including  those 
at  the  crack  tip,  are  displacement  type, 
then  the  above  procedure  directly  applies 
to  fracture  mechanics  problems.  If  the 
stress  intensity  factor  is  included  as  an 
unknown,  then  the  vector  U can  be  treated 
as  a generalized  quantity  which  consists 

of  the  unknown  displacement  and  the  stress 
intensity  factors. 

(2)  In  a nonlinear  material  problem.  Equations  IV-70, 
71,  72,  74,  77  and  78  still  hold.  Equation 
IV-73  holds  in  an  incremental  linear 
representation  provided  E.  is  the 

tangent  material  property  tensor.  Then 
Equation  IV-76  yields  a tangent  stiffness 
matrix.  In  this  case.  Equation  IV-69 
becomes 

MAii  + CAu  + KAu  = AF 

In  the  case  of  a large  displacement/large 
strain  problem,  K must  be  further  modified 
to  include  nonlinear  geometric  effects 
(Zienkiewicz , 1971). 
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(3)  In  case  of  a crack  propagation  problem,  the  net 
effect  is  change  in  the  connectivity  of  the 
finite  element  model.  Appropriate  logic  to 
reorganize  Equation  IV-69  or  Equation  IV-79 
can  be  easily  implemented  in  a computer 
program. 

Crack  Tip  Elements 

It  is  well  known  that  a crack  tip  exhibits  a singular- 
ity. For  a linear  elastic  material,  the  stress  and.  strain 
singularity  is  of  order  1/r  , and  for  the  perfectly  plastic 
material,  the  stress  singularity  is  of  order  1/r 
(Hutchinson,  1968) . This  means  that  theoretically  strains 
are  infinite  at  the  crack  tip,  and  that  they  vary  rapidly 
in  the  vicinity  of  the  crack  tip.  Some  of  the  earlier 
attempts  in  finite  element  application  had  been  with  the 
help  of  conventional  finite  elements  without  the  required 
singularity.  The  elements  were  very  small  in  the  vicinity 
of  the  crack  tip.  However,  it  was  found  that  such  a pro- 
cedure resulted  in  an  unsatisfactory  stress  field  (Hilton 
and  Sih,  1973) . In  addition,  this  procedure  is  very  in- 
efficient computationally  because  of  having  an  unnecessarily 
high  number  of  degrees;  of  freedoms.  The  answer  to  the  problem 
was  development  of  finite  elements  which  would  include  the 
appropriate  singularity.  Use  of  such  elements  enables  one 
to  predict  accurately  stress  or  strain  intensity  factors 
which,  in  turn,  may  be  employed  in  crack  initiation  and 
growth  criteria. 

A number  of  special  crack  tip  finite  elements  have  been 
presented  in  the  literature.  They  can  be  classified  into 
two  groups.  In  the  first  group  the  crack  tip  is  imbedded 
in  the  interior  of  the  element;  in  the  second  group  the 
crack  tip  is  enclosed  by  a number  of  special  elements,  all 
of  which  have  one  of  their  nodal  points  at  the  crack  tip. 

Elements  of  the  first  type  were  first  developed  by 
Wilson  (1968) , and  by  Hilton  and  Hutchinson  (1971) . These 
elements  are  circular  in  shape.  Wilson's  element  had  pre- 
scribed displacements  representing  the  mode  III  (antiplane) 
elastic  singularity.  Hilton  and  Hutchinson's  element  dis- 
placements were  those  associated  with  the  singularity  for 
an  elastic-plastic  power  hardening  material  under  mode  I 
plane  stress  and  mode  III  antiplane  conditions.  Byskov 
(1970)  and  Wilson  (1971)  later  presented  improved  elements 
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by  adding  more  terms  to  the  displacement  function  of  the 
element.  The  displacement  coupling  between  the  circular 
element  boundary  and  its  neighboring  triangular  elements 
is  maintained  by  imposing  a displacement  constraint  on 
the  triangular  elements:  The  displacements  of  nodal  points 
of  the  neighboring  elements  which  lie  upon  the  perimeter  of 
the  circular  element  are  set  equal  to  those  of  the  circular 
element.  A minor  incompatibility  still  persists  - the  edge 
of  the  Wilson  element  is  circular  and  the  displacements  vary 
along  it  linearly  with  the  angle  d,  whereas  the  edge  of  the 
neighboring  element  is  straight,  and  the  displacement  in 
this  element  varies  linearly  along  the  edge.  However,  it  is 
easy  to  see  that  as  the  number  of  nodes  along  the  circle 
increases,  the  incompatibility  diminishes.  Holston  (1976) 
has  extended  the  Wilson  element  to  the  mixed  mode  I and  II 
problem  by  including  additional  terms  in  the  displacement 
expansion . 

Tracey  (1971)  and  Wilson  (1971)  developed  triangular 
elements  which  enclose  the  crack  tip  in  a fan-type  arrange- 
ment (cf.  Figure  IV-74) . They  both  include  a term  propor- 
tional to  in  the  displacement  expansion,  which  gives  a 

1/V^ singularity  in  the  stress  field  appropriate  for  the 
mode  I elastic  singularity.  The  main  difference  between 
the  Tracey  and  the  Wilson  elements  is  in  the  manner  in 


FIGURE  IV-74  Crack  Tip  Enclosed  by 
Triangular  Crack  Tip 
Elements 
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which  the  interpolation  is  performed  in  the  circumferential, 
direction.  Wilson  uses  interpolation  functions  which  are 
linear  in  the  polar  angle  6.  This  approach  causes  a minor 
incompatibility  depending  upon  the  number  of  nodes  on  the 
circle.  Tracey  uses  an  interpolation  function  which  is 
linear  along  the  straight  line  joining  the  two  circumfer- 
ential nodes. 

Wilson  has  shown  that  his  higher  order  circular 
element  (1971)  is  superior  in  performance  to  other  candidates 
discussed  above  (1973) . 

During  the  last  two  years,  a new  breed  of  crack  tip 
elements  has  been  proposed.  Basically,  these  elements  are 
the  regular  isoparametric  elements  (Ergatoudis,  1968)  which, 
with  a special  arrangement  of  nodes,  yields  the  desired 
singularity.  The  first  such  attempt  was  published  by 
Henshell  and  Shaw  (1975) . They  showed  that  a 1/,/T”  stress 
singularity  develops  in  the  8 node  isoparametric  finite 
element  if  the  midside  nodes  of  two  adjacent  edges  are  moved 
to  quarter  points,  as  indicated  in  Figure  IV-75.  Soon  after 
the  paper  by  Henshell  and  Shaw,  Barsoum  (1976a)  published 


FIGURE  IV-75  Eight  Node 

Isoparametric 

Element 
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a similar  paper.  He  showed  that  by  collapsing  one  of  the 
sides  of  the  8 node  isoparametric  elements  and  moving  the 
midside  nodes  of  the  edges  adjacent  to  the  collapsed  edge 
to  the  quarter  point,  a l/^/r~ stress  singularity  is  again 
obtained  (cf . Figure  IV-76) . The  same  approach  is  also 
applicable  to  the  3D  element.  He  later  extended  his  method 
to  plate  bending  and  general  shell  elements  (Barsoum,  1976) . 
Recently,  Barsoum  (1977)  has  suggested  that  if  3 common 
nodes  in  his  degenerate  8 node  element  are  left  unconstrained, 
a 1/r  strain  singularity  results,  which  is  appropriate  for 
mode  I plastic  strain  field.  The  obvious  advantage  of  using 
Barsoum' s approach  is  that  many  computer  programs  already 
have  isoparametric  element.  Therefore,  no  new  program 
development  is  required  in  their  use  in  the  fracture  problem. 

Hellen  (1977)  in  his  discussion  of  papers  by  Henshell 
and  Shaw  (1976)  and  Barsoum  (1976)  proposes  to  use  a special 
isoparametric  element  which  uses  a slightly  different  dis- 
placement field  than  the  one  aimed  at  by  the  referenced 
authors.  Hellen 's  element  gives  results  which  are  almost 
identical  to  those  of  Barsoum. 

For  most  cases  the  degenerate  isoparametric  elements 
discussed  above  appear  to  be  very  attractive.  However, 

neither  these  nor  other  elements  provide  all  the  crack  and 
corner  singularities  that  may  be  encountered  in  practice. 
Therefore  the  work  on  new  elements  continues. 


FIGURE  IV-76  Barsoum' s (1976a)  Eight  Node  Isoparametric 
Element  for  Crack  Tip  Analysis 
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e . Conclusions  on  Fracture  Analysis 

(1)  Several  problems  have  been  solved  using  the 
classical  (analytical)  methods  of  analysis.  The  problems 
solved  are  for  bodies  of  simple  shapes  and  are  otherwise 
limited,  in  scope.  As  noted  in  the  text,  solution  of  such 
problems  is  useful  in  studying  physical  phenomena.  On  the 
other  hand,  these  methods  cannot  be  applied  to  most  realistic 
engineering  problems. 

(2)  It  is  possible  to  solve  a quite  general  prob- 
lem using  the  finite  difference  method.  Application  of  this 
method  has  been  demonstrated  by  the  use  of  HEMP  and  HEMP  3D 
codes.  The  main  problem  with  running  cracks  is  in  the 
definition  of  the  propagation  criterion. 

(3)  There  is  extensive  literature  available  in 
the  application  of  the  finite  element  method  to  simple 
geometry  problems.  Programs  like  MARC  can  be  applied  to 
complicated  dynamic  problems,  after  sufficient  confidence 
has  been  achieved. 

(4)  In  the  case  of  both  the  finite  difference  and 
the  finite  element  methods  of  analysis,  the  material  can  be 
either  elastic  or  elasto-plastic , or  practically  any  other 
type.  Conceptually,  strain  rate  effects  can  be  incorporated. 
However,  the  applications  made  so  far  are  limited  to  the 

use  of  material  properties  for  an  assumed  constant  strain 
rate . 


(5)  Special  crack  tip  finite  elements  have  been 
developed  for  singularities  in  elastic  and  perfectly  plastic 
materials.  By  suitably  locating  the  nodes,  one  can  use  a 
planar  and  three-dimensional  isoparametric  element  for  both 
materials.  Similar  formulation  is  possible  for  the  plate 
bending  and  general  shell  elements.  All  of  these  elements 
are  based  on  asymptotic  expansions  to  represent  the  crack 
tip  singularity.  These  expansions  have  been  established 
so  far  for  only  elastic  and  uniform  plastic  behavior  near 
the  crack  tip,  whereas  in  reality  there  is  a nonuniform 
elastic-plastic  zone. 
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c. 


Scale  Modeling 


One  of  the  three  major  purposes  of  this  study  is  to 
recommend  areas  where  current  practices  should  be  strength- 
ened with  regard  to  small  scale  tests.  It  is  expected  that 
improvements  in  scale  modeling  techniques  will  reduce  the 
need  for  expensive  and  time-consuming  full-scale  testing. 

Modeling  is  a procedure  in  which  similitude  analysis 
is  used  to  design  tests  conducted  to  predict  the  structural 
or  material  behavior  of  a prototype  structure  subjected  to 
various  loadings.  This  procedure  may  be  used  for  a number 
of  purposes:  (1)  evaluate  the  structural  and  material 

limitations  of  the  proposed  prototype;  (2)  experimentally 
evaluate  and  guide  analysis  procedures;  (3)  obtain  response 
information  that  no  other  engineering  techniques  can  provide; 
and  (4)  develop  empirical  damage  rules  to  improve  designs 
to  resist  dynamic  loading. 

Large  and  full-scale  models  are  frequently  employed 
simply  because  of  a lack  of  confidence  in  smaller-scale 
models,  despite  the  cost  and  time  consumption.  This  lack 
is  due  to  a belief  that  small-scale  models  will  not  duplicate 
the  many  critical  structural  and  material  conditions  that 
interact  under  dynamic  loading.  However,  through  careful 
application  of  scaling  principles,  models  can  be  very  effec- 
tive for  the  purposes  listed  above. 

Basically,  there  are  four  types  of  similarities* 

• Geometric  similarity  - The  geometric  dimensions 
of  the  prototype  are  scaled  by  a constant  scale 
factor . 

• Kinematic  similarity  - The  pure  motion,  without 
reference  to  masses  or  forces,  is  duplicated. 

• Dynamic  similarity  - The  responses  such  as  dis- 
placements, velocities,  accelerations ■ forces  , 
strain  rates,  etc.  are  scaled.. 

• Constitutive  similarity  - All  relevant  mater  ial 
response  parameters,  such  as  .fracture  toughness, 
are  scaled. 

* See  Baker  et  al.  (1973)  for  a detailed  treatment  of 
theory  and  practice  of  scale  modeling. 
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In  many  cases,  all  similitude  parameters  for  a given  problem 
cannot  be  satisfied  simultaneously.  This  leads  to  the 
selection  of  alternate  types  of  modeling,  such  as  replica 
versus  nonreplica.  By  definition,  a replica  model  has 
geometrical  similarity  and  is  constructed  of  the  same 
material  as  the  prototype.  The  scaling  is  in  size  alone. 
Replica  models  are  valid  for  a broad  range  of  material  behav- 
ior, including  linear  elastic,  nonlinear  elastic  and  elas- 
tic-plastic; fracture  as  well  may  be  included  when  it  can 
be  defined  by  global  quantities  such  as  maximum  stress  or 
maximum  strain.  However,  in  replica  modeling  certain  impor- 
tant parameters  are  not  scaled  properly;  the  most  signifi- 
cant of  these  relate  to  strain-rate  dependence  and  fracture 
toughness.  Nonreplica  models,  by  definition,  use  materials 
which  are  different  from  those  in  the  prototype.  This 
feature  provides  flexibility  in  fabricating  the  models  and 
in  satisfying  the  required  similitude  parameters. 

If  response  depends  on  fracture  toughness,  a good 
prospect  is  to  seek  materials  which  preserve  the  ratio  of 
the  plastic  zone  size  to  the  plate  thickness;  as  an 
approximation  one  might  use 


where  is  dynamic  fracture  initiation  toughness,  t is 
plate  thickness  and  cr^  is  yield  stress.  A program  of 
experiment  and  analysis  could  be  conducted  to  select  appro- 
priate scaling  parameters  using  this  rule  (or  other  more 
suitable  rules)  and  to  verify  that  model  tests  using  tough- 
ness-scaled materials  will  give  reliable  information  on 
actual  structures  with  flaws . 

The  justification  for  using  small-scale  models  has  been 
greatly  enhanced  by  the  rapid  improvements  in  computer 
coded  analytical  methods  of  predicting  dynamic  response. 

If  the  analysis  can  accurately  predict  the  response  for  the 
model,  there  is  good  reason  to  believe  that  the  prototype 
will  respond  in  a similar  manner,  providing  the  model  cor- 
rectly represents  the  prototype;  depending  on  the  problem, 
it  may  be  necessary  to  account  in  the  model  for  the  effect 
of  imperfections  (flaws  and/or  geometric  distortions) 
and  residual  stresses.  As  discussed  in  earlier  sections. 
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proven  analytical  codes  describing  linear  and  nonlinear 
elastic  behavior  are  numerous,  and  recently  there  has  been 
good  progress  in  developing  analyses  for  elastic-plastic 
response. 

The  capability  to  predict  dynamic  fracture  of  a full- 
scale  fabricated  flawed  structure  solely  by  analytical 
methods  is  not  likely  in  the  near  future.  For  certain 
structures,  nondestructive  testing  (NDT)  procedures, 
applied  in  service,  can  alleviate  this  problem.  In  the 
interim,  it  is  believed  to  be  more  realistic  to  use 
analytical  techniques  to  determine  overall  response  of  a 
structure  and  to  make  a full  thickness  model  of  a smail 
segment,  or  substructure,  for  fracture  toughness  testing. 
In  all  these  situations,  there  is  no  doubt  that  continued 
improvements  in  analysis  directly  enhances  the  use  of 
small-scale  models  and  substructure  models. 
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Metallurgical  Advisory  Board,  through  a series  of  changes  and  expansion  of  scope,  it 
became  the  Materials  Advisory  Board  and,  in  January  1969,  the  National  Materials 
Advisory  Board.  In  consonance  with  the  scope  of  the  two  Academies,  the  general  purpose 
of  the  Board  is  the  advancement  of  materials  science  and  engineering,  in  the  national 
interest.  The  Board  fulfills  its  purpose  by:  providing  advice  and  assistance,  on  request,  to 
government  agencies  and  to  private  organizations  on  matters  of  materials  science  and 
technology  affecting  the  national  interest;  focusing  attention  on  the  materials  a^ects 
of  national  problems  and  opportunities,  both  technical  and  nontechnical  in  nature,  and 
making  appropriate  recommendations  as  to  the  solution  of  such  problems  and  the 
exploitation  of  these  opportunities;  performing  studies  and  critical  analyses  on  mate- 
rials problems  of  a national  scope,  recommending  approaches  to  the  solution  of  these 
problems,  and  providing  continuing  guidance  in  the  implementation  of  resulting 
activities;  identifying  problems  in  the  interactions  of  materials  disciplines  with  other 
technical  functions,  and  defining  approaches  for  the  effective  utilization  of  materials 
technologies;  cooperating  in  the  development  of  advanced  educational  concepts  and 
approaches  in  the  materials  disciplines;  communicating  and  disseminating  information 
on  Board  activities  and  related  national  concerns;  promoting  cooperation  with  and 
among  the  materials-related  professional  societies;  maintaining  an  awareness  of  trends 
and  significant  advances  in  materials  technology,  in  order  to  call  attention  to  opportuni- 
ties and  possible  roadblocks,  and  their  implications  for  other  fields,  and  recognizing  and 
promoting  the  development  and  application  of  advanced  concepts  in  materials  and  ma- 
terials processes. 


i 


